Search

Topographic relief and bedrock fractures caused by fault activity can accelerate surface erosion, suggesting that erosion rates increase with decreasing distance to the fault core. To test this hypothesis, we used 10Be cosmogenic nuclides to quantify catchment-wide erosion rates at various distances from a fault in the footwall of the Daqing Shan fault system, a 200-km-long active normal fault in Inner Mongolia, northern China. In addition, in order to examine the influence of topography and bedrock fractures on the surface erosion processes in the fault-weakened zone, their relationships with hillslope erosion, fluvial sediment grain size and transport efficiency, erosion rates, and erosion coefficients were studied. The resulting catchment-wide erosion rates ranged from 0.01 ± 0.00–0.13 ± 0.01 mm yr−1. These erosion rates and coefficients were unaffected by distance from the fault. The analysis of topography, bedrock fractures, and hillslope erosion yielded distinct hillslope erosion processes at different distances from the fault. Rapid mass wasting events, such as landslides, occurred at <15 km from the fault in steep, rocky hillslopes with abundant bedrock fractures. In contrast, at >15 km from the fault, hillslopes gradually become more gentle and soil-mantled. These spatial-slope configuration differences result in a decrease in sediment grain size delivered to river channels with increasing distance from the fault. Once sediments enter river channels, they impact fluvial transport capacity by influencing the initial motion threshold of sediment. Correlation analyses revealed that sediment transport efficiency in channels directly controls the catchment-averaged erosion rate and coefficient. This mechanism explains why the erosion rate remains insensitive to variations in topography and bedrock fractures in fault-weakened zones. • Study of patterns and controls on erosion rate across Daqing Shan fault zone. • Cosmogenic nuclide 10Be used to derive catchment-averaged erosion rate. • Erosion rates do not vary with distance from the fault. • Erosion rate remains insensitive to variations in topography and bedrock fractures. • Sediment transport efficiency in channels controls erosion rate. [ABSTRACT FROM AUTHOR]

• Contributions were quantified of bed erosion, confluence and emission to P transport. • An integrated model covering the flow-sediment-P transport was developed. • Impacts were studied of nonuniform sediment and riverbed coarsening on P transport. Phosphorus (P) transport plays a crucial role in the aquatic ecology of natural rivers. However, our understanding still remains unclear that how P transport is affected in a river-lake connected system downstream of a dam. This system usually undergoes both severe channel degradation and complex exchange of flow-sediment-phosphorus between the mainstem and tributaries. In the current study, a method was proposed firstly to determine the individual contribution of different sources to P recover based on the calculation of phosphorus budget; then an integrated model was developed, covering the modules of flow, nonuniform sediment and phosphorus transport. The application of the proposed method in the 955-km-long Middle Yangtze River (MYR) shows that the type of P transportation was predominantly changed from particulate phosphorus to dissolved phosphorus after the operation of the Three Gorges Project (TGP), but a significant longitudinal recovery of total phosphorus (TP) flux was observed. The TP flux exporting from the MYR was mainly from the Upper Yangtze River (44%), and 12%, 18% and 26% of that were originated from channel erosion, tributary confluence and anthropogenic emission. Moreover, the effects were investigated of nonuniform sediment transport and bed-material coarsening on P transport in the MYR, based on the proposed integrated model. Obtained results show that the TP transport process in the MYR was more reasonable simulated using the nonuniform sediment mode, and it is also confirmed that the process of bed-material coarsening after the TGP operation would lead to the decrease of particulate phosphorus flux in the MYR. [Display omitted] [ABSTRACT FROM AUTHOR]

Climate change is affecting the hydromorphological system. In many places, changes in sediment dynamics are closely correlated to changes in precipitation frequency and magnitude. However, in nivo-pluvial regimes, the hydromorphological response to climate is more challenging to predict as it is not only the amount and occurrence of precipitation that is changing. The changes in precipitation type (i.e., snow or rain), snow accumulation, and snowmelt rates will also have a significant effect on the catchment net precipitation (composed of direct precipitation plus snowmelt contribution), and this may affect overland flow, erosion, stream discharge, sediment transport, and deposition. We investigated the impacts of climate change on hydrology and geomorphology in a small catchment (Emme, 127 km2) located in the Swiss pre-Alps by simulating the difference between the hydromorphological response to net precipitation in the present climate and in three climate scenarios at the end of the century using the CAESAR-Lisflood landscape evolution model. For the most extreme climate scenario (RCP8.5), simulations showed that despite the reduction in net precipitation (by 7 %) and discharge (by 4 %), sediment yield at the outlet of the catchment increased by 6 %. This is not only because precipitation falls more as rain than snow during the cold months, but also because heavy precipitation is expected to intensify. On a seasonal scale, we found that the amount of net precipitation, discharge, and erosion will increase in winter at the end of the century, while it will decrease in spring. In all three climate scenarios, net precipitation is projected to decrease in summer, but sediment yields may both decrease or increase. Autumn is the season with the greatest changes in erosion, while net precipitation remains constant or only slightly increases. Furthermore, we found that erosion and deposition patterns are changing spatially, with more erosion in mid-elevations and more deposition in valleys. Although our results are specific to the study site, we expect similar trends in other catchments of the pre-Alpine region. • Effects of climate change on erosion in a nivo-pluvial catchment were studied. • Erosion rates are foreseen to increase despite a decline in precipitation amounts. • This results from heavy precipitation intensification and a shift from snow to rain. • Erosion rates are found to increase in winter and autumn, and decrease in spring. • Valleys will transition from erosion to deposition, but slope erosion will increase. [ABSTRACT FROM AUTHOR]

Computer-based coupled soilscape-landform evolution models provide an avenue for estimating the erosion potential of a landform over geological timescales. They also provide tools to understand soil-landform relationships in the context of both field and modelled pedogenesis processes. In this study, we used the State Space Soil Production and Assessment Model (SSSPAM) and CAESAR-Lisflood landform evolution model to simulate the evolution of a natural catchment where soil erosion rates have been closely monitored over the last 20 years. Both models predict comparable erosion rates to that of field measurements and produce similar catchment scale geomorphological patterns. Rock content data measured in the field along two transects were compared with SSSPAM simulations at the same representative points. The results showed that SSSPAM produced similar rock content variation along the respective transects. The findings demonstrate 'SSSPAM's ability to predict catchment scale erosion and surface soil distribution. The insights demonstrated here provide confidence in the SSSPAM model for Soilscape predictions for natural, agricultural and post-mining landform assessment. [Display omitted] • SSSPAM and CAESAR-Lisflood models were used to simulate the evolution of a natural catchment. • Both models produced catchment erosion rates comparable with field observations. • SSSPAM reproduced surface rock content variation on 2 transects observed in the field. • SSSPAM also revealed a possible relationship between area, slope and soil grading. • This study demonstrates SSSPAM's ability to predict erosion and surface soil distribution. [ABSTRACT FROM AUTHOR]

Training walls to improve navigability at river mouths alter natural sediment transport regimes, occasioning morphological changes in the nearby littoral zone. In addition, harbor channels and shorelines are susceptible to on-going changes in global circulation conditions and large-scale engineering interventions. Here we present a process-based modeling approach to determine the relationship between spatial patterns of erosion and deposition, seasonal river discharge, and the geometry of coastal defenses. The study uses the Magdalena River delta as an example, where the redistribution of freshwater and sediments after the construction of the "Tajamar" training wall heralded significant morphological changes in the area. The numerical experiments are used to describe how the erosion and deposition patterns within the estuary and delta of the Magdalena River are linked to the seasonal cycle of the Magdalena River and to the geometry of the hard structures. We demonstrate that the construction of the Tajamar and complementary hard structures would have decreased depositional fluxes in the littoral zone, leading to the observed shoreline retreat and disappearance of an extensive coastal lagoon system. It is shown that an aperture in the training walls may help restore the wetlands without compromising navigability. In addition, projected increases in mean sea level are shown to decrease velocities within the lower estuary, potentially causing increased sedimentation within the channel and more complex conditions for the management of safe navigation over the Magdalena River estuary. • Seawall construction on the Magdalena River, Colombia, resulted in the redistribution of freshwater and sediments. • Numerical experiments provide evidence that shoreline retreat was caused by the altered sedimentation regime. • A diversion system in the training walls may help restore sediment supply to the shorelines without risking navigability. • Sea level rise may increase rates of sedimentation within the estuary. [ABSTRACT FROM AUTHOR]

• High-turbidity events are triggered by wind waves with wind speeds of > 7 m s−1. • Sediment erosion and deposition are synchronous during HTEs. • Synchronism causes a near-lakebed turbid water layer and sustains the lutocline. High-turbidity events (HTEs) frequently occur in shallow lakes and can deeply influence lake ecosystems. However, the processes of sediment erosion and deposition are not fully understood. To investigate these processes, two investigative campaigns that integrate bed sediment measurements, sediment settling experiments, in-situ observations, or a hydrodynamics-sediment model were conducted to quantitatively study wind-induced HTEs in Lake Taihu, China. Cohesive sediments were massively entrained from the lakebed by wind waves when the wind speed exceeded its critical value of 7 m s−1, as estimated by water depth, triggering HTEs with turbidity waves of > 100 nephelometric turbidity units. The suspended sediment concentration was the net sediment flux of erosion (average: 5.9 mg m−2 s−1) and deposition (average: 5.5 mg m−2 s−1), which occurred synchronously for most of the HTEs and resulted in high temporal variation in suspended sediment concentration. This synchronism was determined by the critical erosional stress of lakebed sediments, which was spatially heterogeneous (0.017–0.055 N m−2) in a large lake. Equations for erosion and deposition processes were derived following the investigative campaigns and used to develop a cohesive sediment transport model for simulating HTEs in Lake Taihu. Simulations indicate that the synchronism of sediment erosion and deposition led to a near-lakebed turbid water layer during an HTE, which sustained the sediment lutocline in Lake Taihu. Further investigation should be conducted on the role of the near-lakebed turbid water layer during endogenous pollution in large eutrophic lakes. [ABSTRACT FROM AUTHOR]

• CA models for describing runoff are established based on three water distribution rules which significantly impact on confluence. • Dynamic erosion CA models are built by distinguish inter-rill, rill and gully erosions based on critical water depths. • CA models excellently reproduce the runoff and erosion in both the hillslope with rills and the watershed with gullies. • Remain average multi-flow directions algorithm performs the best in runoff and erosion predictions. Rainfall erosion has been calculated and predicted by various soil loss equations, but the progressions of overland flow and sediment transport are less dynamically described. In this study, cellular automata (CA) models based on single flow direction algorithm (SFD), average value multi-flow direction algorithm (AMFD) and remaining average value multi-flow direction algorithm (RAMFD) were constructed. Simultaneously, erosion model was constructed by dividing soil erosion into inter-rill, rill and gully erosions according to critical water depth. The CA models were validated at three scales: a theoretical slope, an in-situ model slope and the natural basin. The results demonstrated the efficient performance of RAMFD in the runoff ascension and recession stages while SFD and AMFD presented unconcentrated runoff distribution, especially in the in-situ slope and basin simulation. Besides, RAMFD model occurred sediment deposition in the basin upstream while SFD and AMFD presented continuous erosion. Furthermore, the runoff Nash-Sutcliffe efficiency (NSE) of three models were 0.85 ∼ 0.95 and 0.65 ∼ 0.94, and the erosion root mean square error (RMSE) were 0.5 ∼ 1.0 kg/min and 0.3 ∼ 0.45 102kg/min in the theoretical slope and natural basin, respectively. Meanwhile, the NSE and RMSE values of RAMFD exhibited the best performance, indicating that this model effectively balanced water distribution while controlling the flow direction to a greater extent. Overall, it is justified to develop an erosion prediction model based on the classification of erosion types, and the rules governing water flow allocation will inevitably result in qualitative and quantitative differences of runoff and erosion. [ABSTRACT FROM AUTHOR]

Submarine canyons are key elements in source-to-sink systems that are commonly developed along continental margins. They act as major conduits transferring sediment and pollutants from continental shelves to deep-water basins, and control the general morphology and evolution of continental margins. This work uses multibeam bathymetric and high-resolution (two- and three-dimensional) seismic data to investigate the main factors controlling the morphology of the Shenhu Canyon System in the northern South China Sea, as well as its detailed morphological character. The Shenhu Canyon System consists of nineteen (19) submarine canyons whose morphologies vary from southwest to northeast along the continental slope. Canyons (C1–C10) in the southwest show greater incision depths, and steeper thalwegs and walls, when compared to their counterparts to the northeast (C11–C17). The southwest canyons are located close to the shelf edge, where the upper continental slope is relatively steep and multiple landslides are imaged. We show that the thalwegs and walls of the southwest canyons were more actively eroded by sediment flows, with respect to the northeast canyons, making them deeper and steeper. Hence, submarine canyons in the southwest, with a more linear geometry, are now directly connected to the Pearl River Canyon. In parallel, seafloor fault scarps act as barriers for sediment transported to the heads of the northeast canyons. This research highlights how sediment supply, sediment pathways, and seafloor scarps can influence submarine canyon morphology along continental slopes. It contributes to a better understanding of the factors controlling canyon morphology worldwide. • The morphology of submarine canyons changes gradually along a continental slope. • High sediment supply and frequent sediment flows promote canyon erosion. • Canyons to the southwest act as preferential pathways for sediment transportation. • Seafloor scarps form physical barriers for sediment transported downslope. [ABSTRACT FROM AUTHOR]

Western Port is a large tide-dominated embayment located 75 km east of Melbourne, Victoria that is flanked by a rapidly growing population, significant industrial infrastructure, grazing lands and protected RAMSAR habitats. Remarkably little is known about the sediment dynamics of Western Port despite concerning rates of erosion of the swamp and floodplain deposits including cohesive clay cliffs that fringe the sunkland. As in other tide-dominated systems, sediment dynamics in Western Port are complicated by strong tidal currents, multiple sediment sources and complex circulation. We applied the coastal sediment cell framework to Western Port, in conjunction with high resolution satellite derived rates of coastal change to elucidate sediment sharing among the sub-embayment beaches. Using end member modelling of 20 sediment assemblages based on elemental and mineralogical composition we robustly identified seven end members that correspond to different point and diffuse provenances. Hierarchical clustering of the 20 beaches dispersed across the bay yielded nine coastal cells, delimited by headlands, deep water channels, estuaries or disruptive infrastructure. Analysis of shoreline change within the cells suggests most cells are neutrally or positively balanced but a northward movement of sediments is leading to rotation. The deep channels provide a pathway for Bass Strait sediments to bypass the embayment headlands and tidal currents enable cross shore exchange with the channels beyond wave derived closure depth in some locations. Understanding the sediment cell boundaries and direction of sediment transport can assist with the management of erosion hotspots; the novel approach applied here can resolve complex connected coastal systems such as estuaries and embayments. • Combined datasets used to backtrack sediment transport processes in a tide dominated bay. • We identified multiple self-contained sediment cells within the embayment. • Transport is in the direction of time-integrated wind balance, causing in-cell beach rotation. • Headland bypassing facilitated through deep channel and strong tidal currents. • Despite local subsidence and sea level rise, shorelines are, on average, prograding. [ABSTRACT FROM AUTHOR]

Subaqueous slopes of many actively prograding deltas are heavily modified by widespread submarine sediment failures, forming highly complex geomorphic features that undergo continuous evolution even over short time scales. However, the effects of wave-seabed interactions (especially liquefaction) on the evolution of seabed sediment failures are rarely studied. Here, we use side-scan sonar, sub-bottom profilers, multibeam bathymetry, and in situ and laboratory geotechnical tests to map the bottom morphology and sediment structure of the Yellow River subaqueous delta. Our results indicate that: 1) The morphology along the Yellow River subaqueous delta is highly complex and can be categorized as a moderately disturbed area and heavily disturbed area; 2) Heavily disturbed area shows notable lines of evidences of liquefication-related sediment failure and four separate morphological zones (Zones 1–4). Zone 1 consists of differently sized sediment highs, which considered to be scarp or monadnock features. Zone 2 undulates with widely distributed tabular sediment blocks, gullies, and small-scale seafloor irregularities. Zone 3 is characterized by irregular, elongate sediment ridges. Finally, Zone 4 is characterized by relatively flat terrain with some localized highs; 3) Sub-bottom profiles used to visualize failure features in Zone 1 show acoustically continuous high-amplitude parallel or subparallel reflection, and the geomorphic features in this area are mainly caused by erosion. In comparison, Zones 2, 3 and 4 lack internal structures and display chaotic/disturbed signatures in corresponding sub-bottom data, which might be an indication of liquefaction; and 4) A combination of preconditioning and triggering factors lead to sediment failures and the complex morphological features in this area; however, silt liquefaction under storm conditions is the most important controlling factor. Based on this evidence, we propose formation mechanism for the liquefaction-related sediment failure features in this area. Sediment liquefaction forms an initial depression and gully when fine particles are resuspended along the arc-shaped liquefaction boundary as sediment vents and re-solidification of the liquefied sediment leads to a stiffer bed. Gullies are then widened and deepened by the predominant currents, and the disturbed depressions develop into elongated highs because of the denser structure and higher strength following liquefaction. While disturbed by engineering activities, the sediment highs become fragmented due to intense localized erosion, having the potential for downslope transport and the formation of sediment ridges. The geometries and distributions of the sediment failure features on the Yellow River subaqueous delta provide insights into process mechanisms and causative factors, especially when compared with other, similar phenomena reported elsewhere. • The geomorphology of the Yellow River subaqueous delta was characterised. • Distinct morphological zones are identified and associated processes explored. • Silt liquefaction under storm conditions dominates. • An evolutionary model for liquefaction-related sediment failure is proposed. [ABSTRACT FROM AUTHOR]

• Effective sediment sizes were determined along a gradient of soil structural stability. • Sediment average size exhibited a unimodal variation with soil structural stability. • Suspension-saltation dominated sediment transport for low and high stable soils. • Suspension-saltation and rolling jointly dominated erosion for medium stable soils. • Soil structural stability determined rainfall intensity effect on sediment transport. Sediment transport patterns (suspension, saltation and rolling) not only depend on flow hydrodynamic characteristics, but also on sediment properties, both of which are largely influenced by soil structural stability. However, how sediment sizes and transport patterns response to original soils over a wide range of structural stability remains unclear. Herein, field plot rainfall simulation experiments on five soils along an increased gradient of structural stability were conducted at two rainfall intensities of 45 and 90 mm h−1 to investigate effective sediment size characteristics and associated transport patterns in sheet erosion. Differing from runoff and sediment yield, sediment sizes were most influenced by soil structural stability (F = 386.9) rather than rainfall intensity (F = 246.7, p < 0.001). Mean weight diameter of sediment (0.16–0.84 mm) and the proportion of coarse particles (>0.25 mm) exhibited a unimodal variation with increased soil structural stability. The effective sediment size distribution shifted from unimodal with the peak at < 0.10 mm for low and high stable soils to bimodal with two peaks at < 0.10 mm and 0.5–1 mm for medium stable soils. Correspondingly, suspension-saltation (>74 %) dominated sediment transport for low and high stable soils, which was determined by fine sized fragments from aggregate breakdown and small flow transport capacity, respectively; while both suspension-saltation (<0.10 mm) and rolling (>0.25 mm) jointly contributed to sediment transport for medium stable soils, which was controlled by the interaction between fragment size and flow transport capacity. Moreover, the positive effect of rainfall intensity on sediment sizes was more remarkable for the medium than for the low and high stable soils. A conceptual model of sediment transport mechanisms in response to soil structural stability is firstly proposed, which merits consideration in the future development of erosion modeling and soil conservation techniques. [ABSTRACT FROM AUTHOR]

• Tool for optimal selection and implementation (TSI) of the eco-remediation measures. • TSI supports implementation of mitigation measures on the most critical source areas. • TSI reduces the size of required area and cost of implementing measures. Erosion processes in the watershed and sediment transport cause hydro-morphological changes, eutrophication, and the loss of storage capacity in reservoirs. This study examines the tool for the optimal selection and implementation (TSI) of eco-remediation measures (ERM) in the river basin area to improve water quality and to reduce storage loss in the reservoir in question. The main purpose of this tool is to support decisions and measures taken to correct defined problems and to improve water quality and storage capacity in the watershed while minimising sediment transport. This tool enables the effective and necessary implementation of these measures to the most critical source areas (CSAs). In order to verify its operability, we selected the Ledava reservoir with a transboundary area of 105.25 km2 in NE Slovenia and SE Austria. With the use of the Soil and Water Assessment Tool (SWAT), critical source areas were determined and the effects of eight different scenarios on sediment yield and load transport were simulated. The results showed that CSAs occupy 12% of the watershed and that sediment inflow into the Ledavsko jezero reservoir could be reduced by up to 30%. After determining the CSAs and which measures would be most effective, the implementation plan could be defined. Within this framework, the TSI enabled the selection of effective measures and contributed to the long-term improvement of the ecological status of surface waters required by Water Framework Directive (2000/60/EC), improving the quality of water bodies of all types to safeguard water ecosystems from harmful consequences. [ABSTRACT FROM AUTHOR]

Abstract The complex interactions between flow hydrodynamics, sediment transport, geometry, and changes in base level control the morphodynamic evolution of tidal channels, inlets, and estuaries. Research on these complex settings have involved decades of field, laboratory, and modeling work. Studies on the one-dimensional morphodynamic of well-mixed, tide-dominated channels demonstrate that the evolution of a horizontal channel subject to tidal fluctuations at the ocean boundary is characterized by erosion in the ocean part of the channel and the formation of a landward migrating shoal that, depending on the channel length, may reach the landward boundary. When the shoal reaches the landward end of the domain, it is reflected and a beach forms as the flow and sediment transport reach conditions of equilibrium. At equilibrium, the elevation of the channel bed profile remains constant over a tidal cycle. Building on previous work, this study focuses on the poorly explored effects of nonuniformity of the sediment size distribution and of the interaction between tidal channel and longshore sediment transport on the equilibrium of tidal channels. The model is verified at laboratory scale, and it is then applied to predict equilibrium conditions with nonuniform sediment and without an imposed input of sediment from the ocean. Model results show that nonuniformity of the grain size distribution of the sediment increases the profile concavity. The equilibrium deposits tend to coarsen from the land to the ocean, and the sediment of the landward deposits presents an upward fining pattern. When the imposed input of sediment from the ocean is relatively small, the equilibrium profile presents a stable shoal at the oceanward part of the channel. In the case of a large input of sediment from the ocean, the channel gradually fills with sediment. The grain size distribution of these channel fill deposits coarsen in the oceanward and upward directions. The comparison between model results and field observations shows that the model can qualitatively capture some of the main features of tidal channels and can thus be of aid in coastal management and restoration with proper site-specific modifications. Highlights • The non-uniformity of the sediment and its caliber play a prime control on estuarine morphodynamics; • Estuarine morphodynamics strongly depends on the interaction between the tidal channel and the longshore sediment transport. • A 1D model of a tidal channel can capture some characteristics of tide dominated estuaries. • The combination of sinusoidal tidal forcing and constant input of sand from the ocean results in cyclic deposition. [ABSTRACT FROM AUTHOR]

Abstract The impact of wind waves and salt marsh erosion on sediment transport dynamics of macro-tidal estuarine systems has been examined using a numerical model. Morphological changes associated with salt marsh erosion facilitate the propagation of waves to the marsh edges, and increase the resuspension and export of sediments from the estuary. Results also highlight the impact of tidally-modulated changes in wave action upon sediment transport dynamics throughout a tidal cycle and from neap to spring tide. In particular, wave action at the landward margin of the estuary is limited by the shallow water depth, and tidally induced increments in water depth only cause a modest increase in wave shear stresses. At the seaward side of the domain, the impact of waves is limited by the deeper water, and tidally induced water level variations have a more significant control on wave shear stresses. The outcomes contribute new insights that are important for the sediment budget of macro-tidal estuaries facing climate change and have implications for the long-term morphological evolution of estuarine wetlands. Highlights • Enhanced sediment export from estuaries caused by wind waves • Enhanced salt marsh erosion caused by wind waves • Salt marsh erosion due to morphological changes of estuaries • Wave impact changes with varying water depth [ABSTRACT FROM AUTHOR]

Abstract In dual wave farms, i.e., arrays of wave energy converters (WECs) with a dual function – generation of renewable power and mitigation of coastal erosion – the spacing between the WECs is a fundamental design parameter. The present research has the objective of establishing how this parameter affects the shoreline evolution behind the array and, on this basis, to propose and apply a method to determine the optimum spacing for coastal protection. The method is demonstrated on a beach subjected to severe erosion. Five case studies are considered: four with different inter-WEC spacings, and one without the wave farm (baseline). A spectral wave propagation model is applied to analyse the variations in significant wave height behind the WEC array. Longshore sediment transport rates are calculated, and a shoreline model is applied. We find that in all the case studies the dry beach area is greater than in the baseline (no farm) case study, which proves the capacity of the dual WEC array to mitigate the erosive trends of the system. Importantly, we obtain that the inter-WEC spacing plays a fundamental role in the evolution of the shoreline and, consequently, in the effectiveness of the WEC array for coastal protection. The case studies with intermediate spacings yield the best performance in terms of dry beach area. More generally, the benefits of dual wave farms in terms of protection of coastal properties and infrastructure, and the ensuing savings in conventional coastal defence measures (coastal structures, beach nourishment, etc.) contribute to the development of wave energy by enhancing its economic viability. The methodology presented in this paper can be used to optimize the design of dual wave farms elsewhere. Graphical Abstract Highlights • Dual wave farms are proposed for carbon-free energy production and coastal defence. • Their coastal defence efficiency is investigated based on the shoreline evolution. • The spacing between wave energy converters plays a fundamental role. • We propose and apply a method to determine the optimum spacing for coastal defence. • Dual wave farms are shown to have potential as a tool for coastal erosion control. [ABSTRACT FROM AUTHOR]

Highlights • Two crucial periods for suspended sediment transport were June and November/December. • Analysis of 236 floods highlighted their intermittency and did not show relationships. • Medium-term continuous time series are necessary to assess range of variations of SSF. Abstract In mountainous catchments, soil erosion and sediment transport are highly variable throughout time and their quantification remains a major challenge for the scientific community. Understanding the temporal patterns and the main controls of sediment yields in these environments requires a long term monitoring of rainfall, runoff and sediment flux. This paper analyses this type of data collected during 7 years (2007–2014), at the outlet of the Galabre River, a 20 km2 watershed, in south eastern France, representative of meso-scale Mediterranean mountainous catchments. This study is based on a hybrid approach using continuous turbidity records and automated total suspended solid sampling to quantify the instantaneous suspended sediment concentrations (SSC), sediment fluxes, event loads and yields. The total suspended sediment yield was 4661 Mg km−2 and was observed during flood events. The two crucial periods for suspended sediment transport at the outlet were June and November/December (63% of the total). The analysis of suspended sediment transport dynamics observed during 236 flood events highlighted their intermittency and did not show any clear relationship between rainfall, discharge and SSC. The most efficient floods were characterised by counter-clockwise hysteresis relationships between SSC and discharges. The floods with complex hysteresis were the more productive in the long term, during this measuring period exceeding a decade. Nevertheless, the current research outlines the need to obtain medium-term (five years) continuous time series to assess the range of variations of suspended sediment fluxes and to outline clearly the seasonality of suspended sediment yields. Results suggest the occurrence of a temporal dis-connectivity in meso-scale catchments over short time-scales between the meteorological forcing and the sediment yields estimated at the outlet. These findings have important methodological impacts for modelling and operational implications for watershed management. [ABSTRACT FROM AUTHOR]

Abstract The lower Amazon basin contains vast floodplains that exchange sediment with the main river. The exchange of sediment between the floodplain and the channel follows a seasonal cycle that is anticorrelated with the hydrological cycle. At low water stages, sediment that has been stored in the floodplain for potentially several thousands of years is eroded and transferred from the floodplain to the mainstem. During high water stages, most sediment transported in the main channel stems directly from the eroding source with little intermittent storage. We apply the meteoric cosmogenic 10Be to stable 9Be ratio (10Be/9Be) as a denudation and weathering proxy to investigate this seasonality in sediment transport. Single meteoric 10Be concentrations ([10Be]) have previously been shown to record floodplain storage; whereas fractions of mobilized 9Be, a trace metal released during weathering, provide degrees of weathering. The resulting 10Be/9Be ratio provides denudation rates of the sediment sources. We compare 10Be/9Be measured on suspended sediment from river depth profiles at two sites (Óbidos/Macapa) and during two hydrological seasons (high/low water stage in 2013). We show that careful construction of river depth profiles in combination with Rouse modeling is necessary to exclude the grain size dependence that together with hydrodynamic sorting introduces a bias into single [10Be] and [9Be]. In contrast, the 10Be/9Be is largely uniform across grain sizes and is thus unaffected by sorting. Our 10Be/9Be ratios, the denudation rates derived from them, and mobilized 9Be fractions do not show variations outside their uncertainties across seasons or sites. Thus, 10Be/9Be show that upstream sediment sources were well mixed by multiple cycles of storage and remobilization in the floodplain. Highlights • Meteoric cosmogenic 10Be/9Be record denudation rates of Amazon sediment sources • No seasonal variations in 10Be/9Be resolvable due to sediment mixing • Depth-integrated single 10Be show small seasonal variations due to sediment storage. • River depth integration required for single 10Be to exclude fluvial sorting bias. [ABSTRACT FROM AUTHOR]

Highlights • Erosion rates in the MTJ region vary from 0.21 to 0.69 mm/yr. • Erosion rates are higher in drainage basins close to the MTJ. • Erosion rates are similar to uplift rates inferred for the past 96–305 ka. • Erosion rates are lower than uplift rates inferred for the past ∼72 ka. • The MTJ region landscape shows a transient response to an increase in uplift rates. Abstract The Mendocino Triple Junction (MTJ) region in northern California is an archetype for studying landscape response to varying rock uplift rates as they increase toward the triple junction. Underpinning these studies lies the assumption that this landscape has reached a dynamic equilibrium in which rock uplift and erosion rates are equal; however, no study has shown that such an equilibrium actually exists around the MTJ region. Here, we report 10Be- and 26Al-derived erosion rates calculated from isotope concentrations in detrital, fluvial sediment from coastal drainage basins; we then compare these rates with uplift rates inferred from marine terraces that were formed and preserved by uplift during the last ∼305 ka. Erosion rates in the more slowly uplifting southern part of the region range from 0.21–0.32 mm/yr and are consistent with rock uplift rates since 305 ka. However, in the northern transition zone, where uplift rates apparently started to increase about ∼100 ka due to northward migration of the MTJ, erosion rates are higher, 0.43 to 0.69 mm/yr. These rates are similar to uplift rates from 96–305 ka, but substantially less than recent uplift rates of ∼3.5–4 mm/yr inferred for the past ∼72 ka. In the central part of the King Range, finite erosion rates cannot be determined due to (i) the presence of excess 10Be that does not appear to have originated from in-situ cosmic-ray production during erosion and transport of the sediment, and (ii) 26Al concentrations below detection limits. The difference between catchment-wide erosion and recent uplift rates in coastal catchments may be related to a lag in response of hillslopes and tributaries to changes in rates of tectonically-driven relative base-level fall or the uncertainties in rock uplift rates inferred from marine terraces. This study suggests that measurements of both erosion and uplift rates are crucial in assessing the equilibrium state of landscapes and in understanding the topographic features made by surface and tectonic processes. [ABSTRACT FROM AUTHOR]

In this paper, we examine the effect of vegetation on soil erosion, runoff generation and sediment transport on saline rangeland hillslopes. Rainfall simulations were conducted at a fixed 114 mm/h intensity on 6 m × 2 m erosion plots with varying degrees of vegetation. Plots were grouped into three categories (L, M and H) based on their canopy cover (L: <5%, M: 5–19%, H: >19%) and selected to limit variations in slope across canopy cover groups. Runoff and sediment samples were combined with three dimensional (3D) reconstruction data used to monitor soil surface microtopographic changes. Runoff initiation was significantly delayed on the L plots but cumulative runoff after 20 min of rainfall simulation indicated a positive effect of vegetation on infiltration processes. Cumulative sediment after 20 min of rainfall was similar across vegetation cover categories. The 3D data suggest that vegetation reduced net sediment delivery from the plots by primarily increasing opportunities for deposition while marginally affecting gross soil erosion. Plots with H vegetation cover experienced lower 3D-estimated erosion volumes but average erosion depth on L plots was lower than that on plots with M and H vegetation covers. Lower runoff volumes on M and H plots may have been compensated by greater runoff erosivity on these plots as runoff was concentrated in a narrower inter-patch space compared to L plots. This study highlights the need for an increased integration between traditional runoff measurement techniques and 3D reconstruction methods. [ABSTRACT FROM AUTHOR]

For wave energy to become a fully-fledged renewable and thus contribute to the much-needed decarbonisation of the energy mix, the effects of wave farms (arrays of wave energy converters) on coastal systems must be addressed. The objective of this work is to investigate the effects of wave farms on the longshore sediment transport and shoreline evolution of a gravel-dominated beach and, in particular, its sensitivity to the longshore position of the farm based on eight scenarios. Nearshore wave propagation patterns are computed by means of a spectral wave propagation model (SWAN), variations in sediment transport rates induced by the farm are calculated, and a one-line model is applied to determine the shoreline position and dry beach area. The significant wave height at breaking is reduced in the lee of the wave farm, dampening sediment transport. We find that changes in the dry beach area induced by the wave farm are highly sensitive to its alongshore position, and may result in: (i) erosion relative to the baseline scenario (without wave farm) in three of the eight scenarios, (ii) accretion in three other scenarios, and (iii) negligible effects in the remaining two. These results prove that the alongshore position of the wave farm controls the response of the beach to the extent that it may shift from accretionary to erosionary, and provide evidence of its effectiveness in countering erosion if appropriately positioned. This effectiveness opens up the possibility of using wave farms not only to generate carbon-free energy but also to manage coastal erosion, thus strengthening the case for the development of wave energy. [ABSTRACT FROM AUTHOR]

Prescribed fire is a common management practice in Mediterranean region to reduce the amount of fuel and to decrease the wildfire risk. The objective of this study is to assess the effect of a prescribed fire on some soil properties, hydrological response and vegetation recovery in experimental plots 5 years after. The results showed that: i) with the exception of electrical conductivity, the prescribed fire did not affect the analyzed soil properties, and ii) overland flow and sediment transport were increased during the first two years, returning to levels pre-fire, 5 years post-fire. The rainfall threshold for overland flow generation was lower during the following months after the prescribed fire owing to the depletion of the vegetation cover. Immediately after the fire, the vegetation cover was of 1.9%, being the three main soil surface components that dominated the hydrological response: charcoal and decayed wood; blackish and greyish ash, and bare soil. However, these areas decreased very rapidly during the second year following the fire, principally, owing to the regrowth of herbaceous plants and shrubs. In addition, the vegetation cover in burned plots was 16.1% higher than that measured in the unburned ones. Overall, the prescribed fire only had an impact on runoff and sediment transport in the two years post-fire, as consequence of vegetation removal. [ABSTRACT FROM AUTHOR]

Velocity-skewed oscillatory sheet flow transport is studied by an approximate model for instantaneous sediment transport resulting from exponential profiles of sediment concentration and velocity over a mobile bed. The phase lag effect, and asymmetry in wave friction factor, oscillatory flow orbital amplitude, roughness height, bed shear stress and wave boundary layer thickness are modelled so that the net boundary layer stream of pure velocity-skewed oscillatory flow can be obtained. The proposed model is found consistent with the previous classical formulas, and qualitatively refines a series of exponents of velocity power function for the approximation of instantaneous sediment transport rate under different flow conditions. The sediment fluxes can be well predicted with accuracy as that of a two-phase numerical model. The net boundary layer stream and sediment transport rates are also well predicted, especially the negative net stream and rates with large phase residual. The phase residual is shown a useful but insufficient factor in the net sediment transport in pure velocity-skewed flows. In addition, the net boundary layer stream resulting from velocity-skewed wave boundary layer development difference between onshore and offshore flow durations is an essential factor as well. [ABSTRACT FROM AUTHOR]

• Demonstration of rise and fall of carbon discharge with time. • Demonstration of enrichment of sand in loose layer above soil. • Demonstration of no carbon enrichment is sediment at the steady state. Soil carbon losses from sheet and interrill erosion areas have an impact on the global carbon cycle. Sediment discharged in many experiments where artificial rainfall has been applied to sloping soil surfaces show greater proportions of carbon in the sediment than found in the cohesive soil source. In many cases, conclusions drawn from these experiments fail to recognise that the enrichment of carbon in the sediment discharge is a result of unsteady conditions that occur before the steady state develops. It is known that particles of sand of different sizes move at different rates when raindrop induced salation occurs with the result that initially the particle size distribution is finer than the source but tends towards that of the source as time progresses. The presence of carbon in aggregates influences the density of aggregates. Here a qualitative model of sediment transport of sand and coal particles by raindrop induced saltation is used to demonstrate time varying enrichment of coal resulting from the effect of carbon on particle settling velocity. The model also demonstrates how the layer of loose particles that develops on top of the cohesive soil surface influences the relative amounts of carbon and sand in the sediment discharges. The steady state is not achieved until the slowest moving particle detached at the farthest distance from the downslope boundary is discharged. That may not be achieved in many of the experiments reported in the literature. [ABSTRACT FROM AUTHOR]

The Sediment Transport and Erosion Prediction (STEP) model is methodically described in this paper. The STEP model was implemented at a 3.47 km × 3.47 km grid resolution in the Kalisapi sub-watershed, Central Java, Indonesia. It was calibrated and validated in the period of 2003–2008 and 2010–2014 respectively. The result shows that the STEP model functions satisfactorily at monthly and annual time window in the calibration and validation periods. In both periods, the STEP model is able to produce fruitful annual sediment transport rates with R2 values of 0.71 and 0.9 in the calibration and validation period respectively. Moreover, the STEP model exhibits more equitable spatial variability than that of the widely known USLE model. Correspondingly, it can be surmised that the STEP model is capable of simulating spatial and temporal variability of sediment transport from a watershed at monthly and annual time scales. • STEP model is developed to simulate soil erosion, deposition and sediment transport processes in time and space. • STEP model has been successfully implemented, calibrated and validated in the Kalisapi sub-watershed, Indonesia. • Calibrating STEP model to fit the observation is straightforward as it only needs to modify two model parameters. • STEP model source code and input data for the study area are freely available in the Bitbucket repository service. [ABSTRACT FROM AUTHOR]

Fluid mud supported by waves can move downslope on a gentle slope under gravity to transport massive sediment across the continental shelf and cause significant geomorphological changes. In order to explore the formation and transportation of fluid mud under extreme conditions, we undertook in-situ measurements and observed fluid-mud-related processes triggered by a typhoon, Danas, in the subaqueous Changjiang Delta. Under the action of waves, firstly, the seabed was eroded slowly for ∼930 min caused by the wave-enhanced bed shear stress; next, rapid erosion occurred when the bed erodibility M e / ρ B doubled the value during no waves. The great change in bed erodibility implied that bed liquefaction occurred due to waves. The fluid mud occurred concurrently with the strongest waves and went through five stages with an interval of 0.5–3.5 h between each stage in which waves weakened continuously. In addition, the five fluid mud stages all underwent three identical phases, i.e. , formation, stabilization, and decay. In the formation and decay phases, the bottom friction force exceeded the buoyancy gravitational force, and the gradient Richardson number (Ri) was lower than the typical threshold of 0.25, indicating that fluid mud was unstable due to high turbulence energy enhanced by Kelvin-Helmholtz instabilities. During the stabilization phase, the buoyancy gravity force was almost balanced with the drag friction force, and Ri was close to or higher than 0.25, suggesting a laminar fluid mud status. Eventually, the fluid mud disappeared when significant wave height decreased to <1.39 m. Our observation results showed that bed erosion and liquefaction both contributed to the formation of the fluid mud, and the fluid mud played an important role in sediment transportation and bed landforms in subaqueous deltas during extreme events. • Wave-induced sediment-bed liquefaction both contributed to sediment resuspension and reduced bed elevation. • Resuspended sediment settled and aggregated near seabed to form fluid mud moving downward along the slope. • Fluid mud was observed to undergo the formation, stabilization, and decay phases. [ABSTRACT FROM AUTHOR]

In this study the numerical model Delft3D is implemented on high-resolution grids to evaluate beach morphological response to groins and submerged breakwaters under sea-breeze conditions. The numerical model is compared with field measurements from three sites on the northern Yucatan coast of Mexico. At these sites, relatively low energy waves with significant wave heights arrive at a steep oblique angle to the coast that drives high sediment transport rates. After calibration of the model, by adjustment of model free parameters, the model results show accurate representations of the diurnal variation of hydrodynamic conditions (wave height and current velocity) driven by the daily sea-breezes. Moreover, the morphodynamic model shows a good representation of the rapid morphological changes, particularly of the accretion pattern caused by a temporary groin and the subsequent recovery following the structure removal. However, erosion patterns at the downdrift areas of structures were not reproduced, owing to the numerical model limitations of simulating swash morphodynamics. Model validation with field observations of the beach response behind submerged breakwaters enabled the model to be used for a parametric study. Following previous work, new equations are developed in the present study that show an improvement for predicting shoreline changes near submerged breakwaters under sea-breeze conditions. Despite the limitations in reproducing downdrift erosion, overall the high-resolution modelling approach satisfactorily simulates the hydrodynamic and the updrift accretion responses to coastal structures on sandy beaches with oblique waves. • Simulations of beach change near groins and submerged breakwaters are performed. • Results are compared to measurements at three sites with high sediment transport rates. • New equations show improved prediction of shoreline change near submerged structures. [ABSTRACT FROM AUTHOR]

[Display omitted] • Temporal evolution of fine sediment deposits on a gravel bar during a two-year period. • Simple physical-based models introduced to better understand fine sediment dynamics. • Flow-related events can yield thick deposits or erosion depending on the discharge-concentration time-series. • Strong wind-related events disperse fine sediments and yield erosion. • Limitation of the image-analysis due to the seasonal variability of the image quality. In this study, the temporal evolution of fine sediment deposits on a gravel bar is analysed during a two-year period. The time-series of fine deposit surface area is obtained using the methodology proposed in Part 1 of this work. By combining the hydro-meteorological data, the main evolutions of the fine deposits observed on the bar surface are found to be triggered by events such as natural floods, dam flushes, and strong wind valley breezes, and are quantified for each of these events. A simple processed-based flow model is proposed to investigate the sediment dynamics during the flood events, which balances the erosional/depositional processes of fine sediments. Variability in the fine deposit evolution after floods is mainly due to the different flow and sediment concentration conditions. For small events for which the gravel bar is not totally submerged, a slight erosion generally occurs. For larger events for which the gravel bar is totally submerged, a large net deposition is often observed. The fine deposit surface area time-series also shows sharp decreases of the surface deposit coverage after some strong valley breeze events, with an evolution of the same order of magnitude than floods. The application of a wind transport model shows that strong wind valley breezes can transport fine particles from the bar surface to the river channel, thus leading to a loss of fine deposits on the bar surface. Also, the deposit initial moisture is a crucial element having impact on the fine sediment transport on the bar surface. Dry deposits are much more likely to be transported than wet deposits. A discussion is also provided on the potential impact of wind and vegetation on the long-term dynamics of fine sediment deposits. [ABSTRACT FROM AUTHOR]

• Physically based 2D hydrodynamic and sediment transport model, CCHE2D. • High resolution mesh, topography photogrammetry, experimental landform evolution. • Watershed hydrology, raindrop soil detachment, displacement, erosion and transport. • Distributed solution, tilled watershed, cornfield tillage ridge and furrow morphology. • Simulated runoff, sediment transport, landform evolution; channel network analysis. Soil and gully erosion is a critical threat to sustainable agricultural land management and productivity. Soil erosion by rain splash, sheet runoff, and concentrated flow are complex problems conditioned by the combined interactions of soil physical properties, hydrology, human activity, landscape topography, and climate. Effective management of erosion processes driven by rainfall and surface runoff requires a combined effort of field observation, physical experimentation, and numerical simulation. A physically-based hydrodynamic numerical model, CCHE2D with updated raindrop erosion effects, is applied to simulate landscape evolution processes as a result of raindrop impact and overland flow. The numerical model simulates a thin layer of runoff and concentrated flow over complex terrains, and predicts sediment transport due to raindrop splash detachment, transport and soil surface erosions. Soil erosion and landform evolution simulations due to rainfall are presented for an experimental landscape and a tilled large scale agricultural field. Photogrammetric field data were collected using an unmanned aerial system (UAS) following planting and subsequent rainstorms. Numerical simulation results agreed well with the observed experimental soil surface topography, sediment yield, and drainage network developments. The simulated soil and gully erosion in the tilled field also agreed well with observations. This research successfully demonstrates how high resolution spatial and temporal topographic measurements in experimental and field landscapes can be simulated and interrogated by 2D physically-based hydrodynamic and morphodynamic models. The accurate predictions of tillage field and gully erosion process in agriculture watersheds with complex tillage conditions also indicated the usefulness of the model for studying soil and landscape degradation. [ABSTRACT FROM AUTHOR]

Coastal dunes provide many ecosystem services including protection of infrastructure from wave overtopping and habitat for native species. Foredunes grow at different rates and assume different forms (i.e., short and wide to tall and narrow) depending on a range of factors including pre-existing beach and dune morphology, wind, wave, and water levels, sediment grain size, and vegetation characteristics, yet the relative importance of these factors on foredune growth is understudied. Here, we quantify foredune evolution (2016–2019) and explore the relative influence of a suite of metocean, sedimentary, and ecological factors for dune growth on three barrier islands in the U.S. North Carolina Outer Banks (Cape Lookout National Seashore). We incorporate observed and hindcast wind, wave, and water level data into the process-based, coupled beach-dune evolution model, Windsurf, to explore the relative contribution of factors likely to influence foredune growth at the annual timescale (2016–2017). Our cross-shore topographic profile observations show varied interannual foredune change rates and characteristics, including horizontal retreat and progradation at the dune toe and vertical erosion and accretion at the dune crest. Model results indicate that, of the factors explored, pre-existing morphology had the greatest influence on the type of foredune growth that occurred (i.e., incipient dune development, widening of the dune, and/or vertical accretion), which dramatically altered the final shape of the dune. Variations in wind and wave climates were associated with the relative contribution of marine- and aeolian-driven bed elevation changes and were particularly influential during storms. In addition, increases in the minimum elevation of vegetation on the dune profile (analogous to the cross-shore distance between the shoreline and established vegetation line, i.e., the vegetation limit, as coined by Duran and Moore, 2013) increased dune crest height and dune volume. Moreover, variations in the minimum vegetation elevation resulted in a larger range of dune crest elevations and dune volumes than differences in the median sand grain size. We suggest that insights into the relative influence of metocean, sedimentary, and ecological factors to dune growth can assist in the development of practical coastal management strategies. • Pre-existing beach-dune morphology had the largest influence on dune growth type. • Exposure to the dominant wave direction led to increased dune volume. • Higher minimum vegetation elevation and finer sand resulted in taller dunes. • Minimum vegetation elevation influenced dune growth more than sand grain size. [ABSTRACT FROM AUTHOR]

• Mudstone deposits behave like noncohesive sediments in terms of erosion resistance. • High erosion resistance of laterites cannot be predicted merely by particle size. • Shields' parameter predicted by water content, clay content, and plasticity index. Soil erosion and sedimentation are major factors in landscape evolution and river morphology. The erosion resistance of sediments plays an important role in landscape evolution processes because it determines the threshold of an acting force at the initiative sediment motion, and the amount of sediment transports. This study examined the sediment properties and their correlations to the erosion resistance of mudstone deposits and laterites, which are two representative natural cohesive soils in hill and terrace topography in Taiwan. Soil physical properties and pore water chemical properties, and the erosion resistance were investigated and determined using laboratory tests and hydraulic flume experiments, respectively. Investigations of sediment properties showed that the mudstone deposits and laterites had significant differences in median particle size, soil texture, level of plasticity, internal structure, and pH value and electrical conductivity of their aqueous slurry. Results from the hydraulic flume experiments showed that the laterites had stronger erosion resistance, indicating by larger values of critical shear stress and the Shields' parameter, than the mudstone deposits. Through the stepwise multiple regression analysis, a quantitative equation was proposed to estimate the Shields' parameter using fine sediment properties including water content, clay content, and the plasticity index expressed in decimal fractions, of the two soils. Despite more experiments of the two natural soils with different ranges of the sediment properties are needed to validate the empirical equation for predicting the Shields' parameter, the proposed methodology provides an essential reference in which erosion resistance is an important parameter. [ABSTRACT FROM AUTHOR]

The effects of climate change on sediment yield and transport dynamics in cold climate regions are not well understood or reported. In this study, the Soil and Water Assessment Tool (SWAT) has been built-up, calibrated, and validated against streamflow and sediment load at several monitoring stations in a cold climate region watershed - the Athabasca River Basin (ARB) in Alberta, Canada. The model was then fed with bias-corrected spatial disaggregated high-resolution (~ 10 km) future climate data from three climate models for two emission scenarios (RCP 4.5 and 8.5), and two periods (mid- and end-century). Results show that channel erosion and deposition are the dominant processes over hill slope erosion in the basin. On average, a predicted warmer and wetter future climate has both synergetic and offsetting effects on sediment yield. Changes are sub-region specific and land-use type dependent, thus reflecting a marked spatial and temporal heterogeneity within the basin. Increases on sediment yield in future periods in the agricultural areas are up to 0.94 t/ha/yr, and are greater than reported soil formation rates in the region. Similarly, while substantial increases (by more than two fold) in the sediment load transport through the river reaches were obtained, the changes show both temporal and spatial variability, and are closely aligned with the trend of stream flows. We believe that availability of such models and knowledge of the effect of future climatic conditions would help water managers formulate appropriate scenarios to manage such basins in a holistic way. However, significant uncertainties in future sediment yield and transport, as a result of variations in climatic forcing of different climate models, need to be considered in any adaptation measures. [ABSTRACT FROM AUTHOR]

Sediment connectivity has been shown in recent years to explain how the watershed configuration controls sediment transport. However, we find no studies develop a watershed erosion modeling framework based on sediment connectivity, and few, if any, studies have quantified sediment connectivity for gently rolling systems. We develop a new predictive sediment connectivity model that relies on the intersecting probabilities for sediment supply, detachment, transport, and buffers to sediment transport, which is integrated in a watershed erosion model framework. The model predicts sediment flux temporally and spatially across a watershed using field reconnaissance results, a high-resolution digital elevation models, a hydrologic model, and shear-based erosion formulae. Model results validate the capability of the model to predict erosion pathways causing sediment connectivity. More notably, disconnectivity dominates the gently rolling watershed across all morphologic levels of the uplands, including, microtopography from low energy undulating surfaces across the landscape, swales and gullies only active in the highest events, karst sinkholes that disconnect drainage areas, and floodplains that de-couple the hillslopes from the stream corridor. Results show that sediment connectivity is predicted for about 2% or more the watershed’s area 37 days of the year, with the remaining days showing very little or no connectivity. Only 12.8 ± 0.7% of the gently rolling watershed shows sediment connectivity on the wettest day of the study year. Results also highlight the importance of urban/suburban sediment pathways in gently rolling watersheds, and dynamic and longitudinal distributions of sediment connectivity might be further investigated in future work. We suggest the method herein provides the modeler with an added tool to account for sediment transport criteria and has the potential to reduce computational costs in watershed erosion modeling. [ABSTRACT FROM AUTHOR]

Underneath rock slopes of rubble mound structures often one or more granular filter layers are present. These filter layers prevent base material washout. In contrast to traditional filters, geometrically open filters allow for some movement of base material. In order to design such open filters the amount of erosion and accretion of base material needs to be predicted. Based on 2D physical model tests Van Gent and Wolters (2015) presented a method to estimate the erosion and accretion of sand underneath rock slopes. Since the developed method was only valid for perpendicular wave attack, a new study has been performed to assess effects of oblique wave attack on the amount of erosion and accretion of sand underneath a rock slope. For that purpose 3D physical model tests were performed in a wave basin with 5 different wave angle between 0° (perpendicular wave attack) and 75°. Based on the test results it was shown how the effects of wave obliquity can be taken into account in the design of open filters. In addition, the effect of the storm duration on open filters has been studied. This led to an extended design method to predict the amount of erosion and accretion of sand underneath rock slopes under wave loading. [ABSTRACT FROM AUTHOR]

Examination of large scale, alluvial estuarine morphology and associated time evolution is of particular importance regarding management of channel navigability, ecosystem, etc. In this work, we analyze morphological evolution and changes of the channel-shoal system in the Changjiang Estuary, a river- and tide-controlled coastal plain estuary, based on bathymetric data between 1958 and 2016. We see that its channel-shoal pattern is featured by meandering and bifurcated channels persisting over decades. In the vertical direction, hypsometry curves show that the sand bars and shoals are continuously accreted while the deep channels are eroded, leading to narrower and deeper estuarine channels. Intensive human activities in terms of reclamation, embankment, and dredging play a profound role in controlling the decadal morphological evolution by stabilizing coastlines and narrowing channels. Even though, the present Changjiang Estuary is still a pretty wide and shallow system with channel width-to-depth ratios >1000, much larger than usual fluvial rivers and small estuaries. In-depth analysis suggests that the Changjiang Estuary as a whole exhibited an overall deposition trend over 59 years, i.e., a net deposition volume of 8.3 × 10 8  m 3 . Spatially, the pan-South Branch was net eroded by 9.7 × 10 8  m 3 whereas the mouth bar zone was net deposited by 18 × 10 8  m 3 , suggesting that the mouth bar zone is a major sediment sink. Over time there is no directional deposition or erosion trend in the interval though riverine sediment supply has decreased by 2/3 since the mid-1980s. We infer that the pan-South Branch is more fluvial-controlled therefore its morphology responds to riverine sediment load reduction fast while the mouth bar zone is more controlled by both river and tides that its morphological response lags to riverine sediment supply changes at a time scale >10 years, which is an issue largely ignored in previous studies. We argue that the time lag effect needs particular consideration in projecting future estuarine morphological changes under a low sediment supply regime and sea-level rise. Overall, the findings in this work can have implications on management of estuarine ecosystem, navigation channel and coastal flooding in general. [ABSTRACT FROM AUTHOR]

Lateral inputs from hillslopes and tributaries exert a variable impact upon the longitudinal connectivity of sediment transfer in river systems with differing drainage network configurations. Network topology influences channel slope and confinement at confluence zones, thereby affecting patterns of sediment storage and the conveyance of sediments through catchments. Rates of disturbance response, patterns of sediment propagation, and the implications for connectivity and recovery were assessed in two neighbouring catchments with differing network configurations on the East Cape of New Zealand. Both catchments were subject to forest clearing in the late 1940s and a major cyclonic storm in 1988. However, reconstruction of landslide runout pathways, and characterization of connectivity using a Tokunaga framework, demonstrates different patterns and rates of sediment transfer and storage in a dendritic network relative to a more elongate, herringbone drainage network. The dendritic network has a higher rate of sediment transfer between storage sites in successive Strahler orders, whereas longitudinal connectivity along the fourth-order mainstem is disrupted by lateral sediment inputs from multiple low-order tributaries in the more elongate, herringbone network. In both cases the most dynamic (‘hotspot’) reaches are associated with a high degree of network side-branching. [ABSTRACT FROM AUTHOR]

The 2011 Great East Japan Tsunami caused severe morphological changes along the affected coastline on the northeast part of Japan. This study investigates the recovery of river mouths and breaching of sandy coast at the place used to be river mouth in Miyagi Prefecture, Japan after severe erosion induced by the tsunami through analyzing aerial photograph and investigating field data. After several years since the tsunami took place, the morphology of some river mouths has mostly recovered while others still have some damage. Results on the analysis of aerial photograph clearly indicate that the total length of sandy coast on both sides of river mouth plays important role in the recovery process. Moreover, the relevant phenomenon occurring during the morphological recovery process is also investigated. The findings demonstrated the intrusion of sand spit into river mouths due to the tsunami having scoured river mouth deeper than the depth of closure. Interestingly, the intrusion distance was found to be positively correlated with the width of river mouth. In addition, there is an important finding on the practical problem which is related to the intrusion upstream of sand spit. That intrusion shifts the boundary of fresh-sea water upstream, therefore external wave forces must be considered when implementing the design of river structures. [ABSTRACT FROM AUTHOR]

Relation among estuarine processes, mixing mechanisms and turbidity maximum zone (TMZ) formation was evaluated in a microtidal tropical estuary of high fluvial discharge (Magdalena River, northwestern South America). Particularly, assessing the effects of vertical stratification over mixing, estuarine circulation and suspended sediment transport processes. Measurements of salinity, temperature, density, currents and suspended sediment concentration (SSC) were performed at the Magdalena river mouth during high and low streamflows to distinguish among two contrasting seasonal conditions (flow conditions, wind patterns, wave energy). Salinity, temperature, density and location of the freshwater-saltwater interface (FSI) exhibited a sharp seasonal variability. The river mouth and delta front were stratified regardless of the hydrologic season, whereas the riverine section was stratified during the low streamflow season, reflecting the balance between fluvial advection and vertical stratification. Reduction of the vertical mixing led to a more clearly defined and extended stratification. Bottom friction is the main mixing mechanism, but as stratification is strengthened along the river mouth, the vertical mixing weakens, promoting the developing of TMZ (SSC >4500 mg l −1 ). Our results indicate that TMZ developing is caused by the convergence of stratified layers near the nodal point and the weakening of bottom turbulence by stratification. Although the TMZ exhibited a strong seasonal variation in its magnitude and extent, it usually formed in the convergence front demarcated by saltwater intrusion. Streamflow is the major influence on the TMZ formation since it allows the development of stratified conditions, which weaken turbulent mixing, and thus promotes particulate matter settling. The Magdalena River estuary can be classified as turbid (4000 mg l −1  < SSC max  < 10,000 mg l −1 ) and extremely turbid during low streamflow seasons (SSC max  > 10,000 mg l −1 ). These values are of the same order as the SSC reported for the Yangtze and Yellow estuaries, classified as having the largest SSC values of the world. [ABSTRACT FROM AUTHOR]

Mega-nourishments, aiming at providing long-term coastal safety, nature qualities and recreational space, have been applied recently at the Holland coast and are considered at various other places in the world. Methods to quickly evaluate the potential and lifetime of these coastal mega nourishments are therefore very much desired, which is the main objective of this research. Two types of mega nourishments can be distinguished: feeder-type mega nourishments may erode freely to feed adjacent coasts for a more natural, dynamic dune growth while permanent mega-nourishments are designed to preserve safety levels and need to maintain their size and shape and thus needs to be nourished themselves. The design and impact assessment studies for both types of mega nourishments require detailed morphological studies to determine the morphological evolution. In this paper 2DH (Delft3D) and 1D (UNIBEST-CL+ and LONGMOR) numerical models were calibrated using data of the Sand Motor mega-nourishment and were then applied to model a series of mega-nourishments with various width over length ratios and volumes in order to derive relations and design graphs for erosion rates, life span and maintenance volumes. These relations and design graphs can be used in project initiation phases and feasibility studies. The magnitude of the modelled wave-driven longshore sediment transport rates in 1D coastline models depend on the representation of wave refraction on the lower shoreface, since a distinction should be made between the non-rotating lower shoreface and active surfzone. It was shown that the life time of nourishments is mainly determined by the dimensions of the nourishment and incoming wave energy. [ABSTRACT FROM AUTHOR]

The process of tree uprooting plays an important role in sediment transport within forested slopes. In this study we determine the intensity of sediment transport by uprooting and its dependencies. We use the Principal Component Analysis (PCA) to determine the most important factors influencing root plate sizes, and we calculate sediment flux by windthrow for the area of the Tatra National Park (211.6 km 2 ), for the year 2013, when an extreme windthrow event occurred. For this purpose we combine field data with the results of windthrow areas mapping based on satellite imagery. We also use data recorded by the Tatra National Park concerning the amount of destroyed wood in the forests of the Tatra National Park to calculate sediment flux by windthrow for the period between 1968 and 2015 (48 years). To calculate the transport distance by the uprooting needed to compute sediment flux, we modified the Gabet et al.'s (2003) Eq. (11), which in our opinion gives results which are twice the actual size. The PCA identified two the most important factors influencing root plate sizes. The first factor shows that all root plate dimensions increase with increasing tree DBH (diameter at breast height), and, to a lesser degree, with an increase in the content of coarse material. The second factor shows that the depth of a root plate increases with the downslope angle of tree fall, decreasing slope inclination, and decreasing the content of coarse material. The sediment flux by windthrow in 2013 in the area of the Tatra National Park was 3.55 × 10 − 4 m 3 m − 1 yr − 1 . Mean annual sediment flux for the period of the last 48 years was 2.76 × 10 − 5 m 3 m − 1 yr − 1 . [ABSTRACT FROM AUTHOR]

Silver nanoparticles (AgNPs) are being used in many products as they have unique antimicrobial-biocidal properties. After disposal of these products AgNPs can reach the soil environment possibly affecting soil organisms and disrupting plants. This work aimed to study the transport of AgNPs by water and sediment during overland flow and soil erosion. This was done in a laboratory setting, using a flume and rainfall simulator. A low concentration of AgNPs (50 μg·kg − 1 ) was applied to two soil-flumes with slope percentages of 20% and 10%. The rainfall was applied in four events of 15 min each with a total amount of rainfall of 15 mm during each event. After applying the rainfall, samples of the non-transported background soil (BS) and the transported sediment (Sf) were collected from the flume surface. Runoff sediment (RS) and water (RW) were collected from the outlet. AgNPs were detected in all samples collected. However, concentration varied according to sample type (soil or water), time of collection (for runoff water and sediment) and the slope of the soil flume. Higher concentrations of AgNPs in soil were detected in the BS than in the Sf likely due to the BS having more fine particles (silt and clay). The AgNPs concentration in the runoff sediments increased with subsequent applied rain events. In addition, increasing the slope of the flume from 10% to 20% increased the total AgNPs transported with the runoff sediment by a factor 1.5. The study confirms that AgNPs can be transported by both overland flow and sediment due to erosion. [ABSTRACT FROM AUTHOR]

Posidonia oceanica plays a significant role in the stabilization and protection of the coast in Gulf of Giens. Unfortunately, its distribution has been declining remarkably due to both anthropogenic interventions and natural factors. The present study focuses on the numerical simulation of the presence of Posidonia as well as the influence of its disappearance on hydrodynamics and sediment transport along Alamanarre beach. The model results indicate that the regression of Posidonia leads to significant changes of the hydrodynamic parameters, sedimentation rates, and causes severe erosion along Almanarre beach. In the annual wind conditions, the mean current speed is increased by 50.38%–80.65%, the mean significant wave height increased by 14.03%–16.58%, and the total load of sediment increased by 252%–426%, in absence of Posidonia. Regarding the seasonal wind variation, the disappearance of Posidonia induces an increase of the mean current speed about 52.86%–58.34%, the mean significant wave height by 8.23%–11.27%, and sediment transport rate by 358.99%–394.53%. Under the impact of extreme events with the disappearance of Posidonia, the mean current speed, the mean significant wave height, and the sediment transport rate increase by 12.53%–58.48%, 23.91%–34.67%, and 65.15%–154.62%, respectively. [ABSTRACT FROM AUTHOR]

Efficient mitigation of environmental loads from agricultural fields to surface waters requires sufficient knowledge of dominating erosion processes and sediment transport pathways. Computational models are often applied in sediment load assessments but their structure inherently includes assumptions and uncertainties, which can hinder their predictive and explanatory capabilities. In this study, a 3D dual-permeability model was applied with field-scale data from a high-latitude site to investigate sediment balances and structural uncertainties in sediment transport components. The two-year data encompassed hourly records of water flow and sediment concentration composite samples from tillage layer runoff and drain discharge in two adjacent clayey fields with different slopes (1% and 5%). Three model structures with different assumptions of sediment transport pathways were built to test their performance against the data. The simulations demonstrated how different model structures can reproduce the data with varying results on sediment balance. The varying results revealed the importance of flow, erosion, and sediment transport observations from the fields to improve the simulations. Structural uncertainty analysis revealed uncertainties which parameter sensitivity analysis could not describe. Concentration data was shown to include more information about erosion and sediment transport processes than solely the load data. The results suggest that a major part (48–69%) of the detached particles remained in the field and that lateral subsurface transport contributed to load generation (10–21% of total loads) especially in the steeper field. The results demonstrated that the majority (84–87%) of sediment loads occurred via subsurface drain discharge and groundwater outflow with the slope gradient of 1–5%, which suggests that load mitigation measures should also be directed to decrease loads via subsurface transport pathways. The simulations demonstrated how transport processes were controlled not only by water flow but also by soil structure. [ABSTRACT FROM AUTHOR]

Reservoirs behind dams act as deposition sites for much of the sediment being transported by rivers. As a result, the downstream river flow can be well below the transport capacity for bed-material. This promotes bed erosion and other geomorphic changes over some length of river located immediately downstream from a dam. These adjustments have been characterized for the Trinity River, TX, downstream of Livingston Dam. Field measurements and results from a 1D numerical model define a 50–60 river kilometer segment of river undergoing bed erosion as the transport capacity for bed material is reestablished. Consequences of this erosion include lowering of the channel bed, reduction in the sediment volume of channel bars, coarsening of sediment on bar tops, steepening of channel banks, and reduction in lateral migration rates of river bends. Repeat surveys of the river long profile reveals that 40 yr of dam closure has produced up to seven meters of channel-bottom incision downstream of the dam, transforming an initially linear profile into a convex-up long profile. The model output matches this observed change, providing confidence that calculated estimates for spatial and temporal changes in bed-material sediment flux can be used to explore the long-term signature of dam influence on the geomorphology of a sand-bed channel. Measurements of channel geometry, profile, lateral migration, and grain size of the lower Trinity River with distance downstream define both the trend and expected variability about the trend associated with the disruption to the bed-material load. [ABSTRACT FROM AUTHOR]

This study is the first to simultaneously compare flow and sediment transport through a blowout and over an adjacent foredune, and the first study of flow within a highly sinuous, slot and cauldron blowout. Flow across the foredune transect is similar to that observed in other studies and is primarily modulated by across-dune vegetation density differences. Flow within the blowout is highly complex and exhibits pronounced accelerations and jet flow. It is characterised by marked helicoidal coherent vortices in the mid-regions, and topographically vertically forced flow out of the cauldron portion of the blowout. Instantaneous sediment transport within the blowout is significant compared to transport onto and/or over the adjacent foredune stoss slope and ridge, with the blowout providing a conduit for suspended sediment to reach the downwind foredune upper stoss slope and crest. Medium term (4 months) aeolian sedimentation data indicates sand is accumulating in the blowout entrance while erosion is taking place throughout the majority of the slot, and deposition is occurring downwind of the cauldron on the foredune ridge. The adjacent lower stoss slope of the foredune is accreting while the upper stoss slope is slightly erosional. Longer term (16 months) pot trap data shows that the majority of foredune upper stoss slope and crest accretion occurs via suspended sediment delivery from the blowout, whereas the majority of the suspended sediment arriving to the well-vegetated foredune stoss slope is deposited on the mid-stoss slope. The results of this study indicate one mechanism of how marked alongshore foredune morphological variability evolves due to the role of blowouts in topographically accelerating flow, and delivering significant aeolian sediment downwind to relatively discrete sections of the foredune. [ABSTRACT FROM AUTHOR]

Digital elevation models (DEMs) built from repeated topographic surveys permit producing DEM of Difference (DoD) that enables assessment of elevation variations and estimation of volumetric changes through time. In the framework of sediment transport studies, DEM differencing enables quantitative and spatially-distributed representation of erosion and deposition within the analyzed time window, at both the channel reach and the catchment scale. In this study, two high-resolution Digital Terrain Models (DTMs) derived from airborne LiDAR data (2 m resolution) acquired in 2005 and 2011 were used to characterize the topographic variations caused by sediment erosion, transport and deposition in two adjacent mountain basins (Gadria and Strimm, Vinschgau - Venosta valley, Eastern Alps, Italy). These catchments were chosen for their contrasting morphology and because they feature different types and intensity of sediment transfer processes. A method based on fuzzy logic, which takes into account spatially variable DTMs uncertainty, was used to derive the DoD of the study area. Volumes of erosion and deposition calculated from the DoD were then compared with post-event field surveys to test the consistency of two independent estimates. Results show an overall agreement between the estimates, with differences due to the intrinsic approximations of the two approaches. The consistency of DoD with post-event estimates encourages the integration of these two methods, whose combined application may permit to overcome the intrinsic limitations of the two estimations. The comparison between 2005 and 2011 DTMs allowed to investigate the relationships between topographic changes and geomorphometric parameters expressing the role of topography on sediment erosion and deposition (i.e., slope and contributing area) and describing the morphology influenced by debris flows and fluvial processes (i.e., curvature). Erosion and deposition relations in the slope-area space display substantial differences between the Gadria and the Strimm basins. While in the former erosion and deposition clusters are reasonably well discriminated, in the latter, characterized by a complex stepped structure, we observe substantial overlapping. Erosion mostly occurred in areas that show persistency of concavity or transformation from convex and flat to concave surfaces, whereas deposition prevailingly took place on convex morphologies. Less expected correspondences between curvature and topographic changes can be explained by the variable sediment transport processes, which are often characterized by alternation of erosion and deposition between different events and even during the same event. [ABSTRACT FROM AUTHOR]

Experimental results involving shoreface nourishment scenarios are presented. The experiments are performed at small scale and the effects of nourishment placement and timing in relation to long-term cross-shore profile development are investigated. Four different nourishment scenarios are tested experimentally, with nourishment positioned both in the trough of a pre-developed profile as well as along the bar of the profile. The results demonstrate that under all scenarios the erosion of the shoreline slows relative to the case without nourishment. The two cases where the nourishment was placed along the bar reduced erosion of the shoreline more compared to the two cases with nourishment in the trough onshore of the bar. Compared to past experiments on nourishment presented in the literature, the present experiments were run for a longer duration, and the concepts of equilibrium profile and the development towards such an equilibrium were investigated. Curiously, the experiments showed a reversal of the sediment transport rate from being primarily offshore directed across the entire profile to onshore directed in the inner surf zone region closest to the shoreline. This reversal of transport rate occurred without changing the incoming wave conditions and resulted in deposition at the shoreline. This either questions the existence of equilibrium beach profiles or at the very least shows that the development towards such an equilibrium will not always be monotonic in time. • Shoreface nourishments are performed at model scale investigating the influence of timing and placement on long term profile development. • Placement of nourishment can disturb the cross-shore profile shape over long time. • After long time running with the same wave climate the transport rate changed from offshore directed to onshore directed in the inner surf zone. • The change in sign was related to cessation of wave-breaking and resulted in sand returning to an otherwise eroding beach. [ABSTRACT FROM AUTHOR]

[Display omitted] • Total volume and peak discharge are the most important drivers of sediment transport. • Antecedent moisture condition is correlated only with the initial discharge of floods. • Hysteresis loops between SSC and Q are analyzed with the hysteresis index (HI). • Negative HI are predominant indicating delayed SS transport and distant SS sources. • In the catchment soil erosion from hillslopes is the dominant process. This study aimed to evaluate the drivers of suspended sediment (SS) transport by analyzing the flood events recorded from July 2010 to June 2011 in the Celone River basin (Apulia, S-E Italy). High-resolution data of rainfall, streamflow (Q), and measurements of suspended sediment concentrations (SSC) were used to carry out the multivariate analysis. Among the initial 18 variables, the Principal Component Analysis revealed that total volume (m3), peak discharge (Qmax, m3/s), and the flood duration were the most important factors influencing SS transport. Total rainfall and total kinetic energy recorded in the gauging station located in the lowland were also important factors, whereas the same rainfall variables recorded in the remote mountainous area of the basin had a minor influence on SS transport. Antecedent moisture conditions were significantly correlated only with the initial discharge of the flood. Hysteresis loops between SSC and Q were analyzed for each flood event by using the hysteresis index (HI). Results showed that this catchment was dominated by anticlockwise loops (negative values of HI) representative of delayed SS transport and of distant SS source areas in terms of the gauging station. These results also suggested that soil erosion from hillslopes was the dominant process compared with channel and bank erosion and that there was no exhaustion of sediment availability in the basin. [ABSTRACT FROM AUTHOR]

[Display omitted] • The SWAT model supports better management of erosion processes. • Computational models can be used with methods of visual recording erosion. • Simple statistical models aligned to SWAT is adequate to analyze erosive processes. • The spatially distributed erosion identifies the susceptible areas to water erosion. Water-driven erosion associated with the transport and deposition of sediments in watersheds may lead to contamination problems causing significant damage to land and infrastructure. The key factors that influence erosive processes are well known and include the precipitation rate, terrain topography, soil type, and vegetation cover. The soil and water assessment tool (SWAT) considers these factors and may be used to identify areas susceptible to soil erosion and to estimate soil loss and sediment transport and deposition. To corroborate this assumption, real data on rill erosion and sediment transport from the Mogi Guaçu River Basin in southeastern Brazil were used to validate SWAT estimates. The study shows that the model is able to identify soil-eroded areas and replicate the scale and variation patterns of sediment production, although not with great precision. Hence, it can successfully be used for prioritizing areas that require action and for evaluating alternative soil management and conservation practices to mitigate erosive processes. [ABSTRACT FROM AUTHOR]

Predicting coastal erosion requires an accurate morphodynamic model. XBeach has been widely adopted to simulate storm-induced coastal erosion. Because of the large number of model input parameters and their uncertainties, sensitivity analysis is a crucial step toward reliable results. This requires running a large number of computationally expensive simulations. Here, we adopt a computationally cost-effective approach based on the Non-Intrusive Polynomial Chaos Expansion method to quantify the sensitivity of XBeach to its input parameters. The method is applied to the coastal erosion event caused by Hurricane Sandy along the barrier islands of New Jersey in the United States. The results show a spatial variation of the model sensitivity with increasing boundary conditions. Parameters' interaction is found to decrease as the magnitude of boundary conditions increases, causing a reduction in the nonlinearity of the model behavior, and consequently leading to a reduction in the uncertainty of the model output. • Sensitivity of XBeach model to its input parameters is investigated. • Sensitivity is quantified using Non-Intrusive Polynomial Chaos Expansions method. • Effect of interactions among input parameters on model behavior is quantified. • Effects of boundary conditions and site characteristics on sensitivity are quantified. [ABSTRACT FROM AUTHOR]