The width of fluvial valley floors is a key parameter to quantifying the morphology of mountain regions. Valley floor width is relevant to diverse fields including sedimentology, fluvial geomorphology, and archaeology. The width of valleys has been argued to depend on climatic and tectonic conditions, on the hydraulics and hydrology of the river channel that forms the valley, and on sediment supply from valley walls. Here, we derive a physically based model that can be used to predict valley width and test it against three different datasets. The model applies to valleys that are carved by a river migrating laterally across the valley floor. We conceptualize river migration as a Poisson process, in which the river changes its direction stochastically at a mean rate determined by hydraulic boundary conditions. This approach yields a characteristic timescale for the river to cross the valley floor from one wall to the other. The valley width can then be determined by integrating the speed of migration over this timescale. For a laterally unconfined river that is not uplifting, the model predicts that the channel-belt width scales with river flow depth. Channel-belt width corresponds to the maximum width of a fluvial valley. We expand the model to include the effects of uplift and lateral sediment supply from valley walls. Both of these effects lead to a decrease in valley width in comparison to the maximum width. We identify a dimensionless number, termed the mobility–uplift number, which is the ratio between the lateral mobility of the river channel and uplift rate. The model predicts two limits: at high values of the mobility–uplift number, the valley evolves to the channel-belt width, whereas it corresponds to the channel width at low values. Between these limits, valley width is linked to the mobility–uplift number by a logarithmic function. As a consequence of the model, valley width increases with increasing drainage area, with a scaling exponent that typically has a value between 0.4 and 0.5, but can also be lower or higher. We compare the model to three independent datasets of valleys in experimental and natural uplifting landscapes and show that it closely predicts the first-order relationship between valley width and the mobility–uplift number.
This study implements the least-squares inversion method for solving the exhumation history from the thermochronologic age–elevation relationship (AER) based on the linear equation among exhumation rate, age and total exhumation from the closure depth to the Earth surface. Modeling experiments suggest significant and systematic influence of initial geothermal model, the a priori exhumation rate and the time interval length on the a posteriori exhumation history. Lessons learned from the experiments include that (i) the modern geothermal gradient can be used for constraining the initial geothermal model, (ii) a relatively high a priori exhumation rate would lead to systematically lower a posteriori exhumation and vice versa, (iii) the variance of the a priori exhumation rate controls the variation in the inverted exhumation history, and (iv) the choice of time interval length should be optimized for resolving the potential temporal changes in exhumation. To mitigate the dependence of inverted erosion history on these initial parameters, we implemented a new stepwise inverse modeling method for optimizing the model parameters by comparing the observed and predicted thermochronologic data and modern geothermal gradients. Finally, method demonstration was performed using four synthetic datasets and three natural examples of different exhumation rates and histories. It is shown that the inverted rock exhumation histories from the synthetic datasets match the whole picture of the “truth”, although the temporal changes in the magnitude of exhumation are underestimated. Modeling of the datasets from natural samples produces geologically reasonable exhumation histories. The code and data used in this work are available on Zenodo (https://doi.org/10.5281/zenodo.10839275).
Understanding the response of coastal barriers to future changes in rates of sea level rise, sediment availability, and storm intensity/frequency is essential for coastal planning, including socioeconomic and ecological management. Identifying drivers of past changes in barrier morphology, as well as barrier sensitivity to these forces, is necessary to accomplish this. Using remote sensing, field, and laboratory analyses, we reconstruct the mesoscale (decades–centuries) evolution of central Fire Island, a portion of a 50 km barrier island fronting Long Island, New York, USA. We find that the configuration of the modern beach and foredune at Fire Island is radically different from the system's relict morphostratigraphy. Central Fire Island is comprised of at least three formerly inlet-divided rotational barriers with distinct subaerial beach and dune–ridge systems that were active prior to the mid-19th century. Varying morphologic states reflected in the relict barriers (e.g., progradational and transgressive) contrast with the modern barrier, which is dominated by a tall and nearly continuous foredune and is relatively static, except for erosion and drowning of its fringing marsh. We suggest that this state shift indicates a transition from a regime dominated by inlet-mediated gradients in alongshore sediment availability to one where human impacts exerted greater influence on island evolution from the late 19th century onward. The retention of some geomorphic capital in Fire Island's relict subaerial features combined with its static nature renders the barrier increasingly susceptible to narrowing and passive submergence. This may lead to an abrupt geomorphic state shift in the future, a veiled vulnerability that may also exist in other stabilized barriers.
The lateral movement of drainage divides is co-influenced by tectonics, lithology, and climate and therefore archives a wealth of geologic and climatic information. It also has wide-ranging implications for topography, the sedimentary record, and biological evolution and thus has drawn much attention in recent years. Several methods have been proposed to determine drainage divides' migration state (direction and rate), including geochronological approaches (e.g., 10Be) and topography-based approaches (e.g., χ−1, whereas the migration rates at the Yingwang Shan are approximately zero. This study indicates that the drainage-divide stability can be determined more accurately using high-resolution topographic data. Furthermore, this study takes the cross-divide differences in the uplift rate of channel heads into account in the measurement of drainage-divide migration rate for the first time.
We show a long-term erosion monitoring of several geomorphologically active gully systems on Little Ice Age lateral moraines in the European Central–Eastern Alps, covering a total time period from 1953 to 2019 and including several survey periods in order to identify corresponding morphodynamic trends. For the implementation, DEM (digital elevation model) of Differences (DoDs) were calculated, based on multitemporal high-resolution digital elevation models from historical aerial images (generated by structure from motion photogrammetry with multi-view stereo) and light detection and ranging from airborne platforms. Two approaches were implemented to achieve the corresponding objectives. First, by calculating linear regression models using the accumulated sediment yield and the corresponding catchment area (on a log–log scale), the range of the variability in the spatial distribution of erosion values within the sites. Second, we use volume calculations to determine the total and the mean sediment yield (as well as erosion rates) of the entire sites. Subsequently, both the sites and the different time periods of both approaches are compared. Based on the slopes of the calculated regression lines, it can be shown that the highest variability in the sediment yield at the sites occurs in the first time period (mainly 1950s to 1970s). This can be attributed to the fact that within some sites the sediment yield per square metre increases clearly more strongly (regression lines with slopes up to 1.5). In contrast, in the later time periods (1970s to mid-2000s and mid-2000s to 2017/2019), there is generally a decrease in 10 out of 12 cases (regression lines with slopes around 1). However, even at sites with an increase in the variability in the sediment yield over time, the earlier high variabilities are no longer reached. This means that the spatial pattern of erosion in the gully heads changes over time as it becomes more uniform. Furthermore, using sediment volume calculations and corresponding erosion rates, we show a generally decreasing trend in geomorphic activity (amount of sediment yield) between the different time periods in 10 out of 12 sites, while 2 sites show an opposite trend, where morphodynamics increase and remain at the same level. Finally, we summarise the results of long-term changes in the morphodynamics of geomorphologically active areas on lateral moraines by presenting the “sediment activity concept”, which, in contrast to theoretical models, is based on actually calculated erosion. The level of geomorphic activity depends strongly on the characteristics of the sites, such as size, slope length, and slope gradient, some of which are associated with deeply incised gullies. It is noticeable that especially areas with influence of dead ice over decades in the lower slope area show high geomorphic activity. Furthermore, we
The size distributions and the shapes of detrital rock clasts can shed light on the environmental history of the clast assemblages and the processes responsible for clast comminution. For example, mechanical fracture due to the stresses imposed on a basal rock surface by a body of flowing glacial ice releases initial “parent” shapes of large blocks of rock from an outcrop, which then are modified by the mechanics of abrasion and fracture during subglacial transport. The latter processes produce subsequent generations of shapes, possibly distinct in form from the parent blocks. A complete understanding of both the processes responsible for block shape changes and the trends in shape adjustment with time and distance away from the source outcrop is lacking. Field data on edge rounding and shape changes of Shap granite blocks (dispersed by Devensian ice eastwards from the outcrop) are used herein to explore the systematic changes in block form with distance from the outcrop.
The degree of edge rounding for individual blocks increases in a punctuated fashion with the distance from the outcrop as blocks fracture repeatedly to introduce new fresh unrounded edges. In contrast, block shape is conservative, with parent blocks fracturing to produce self-similar “child” shapes with distance. Measured block shapes evolve in accord with two well-known models for block fracture mechanics – (1) stochastic and (2) silver ratio models – towards one or the other of these two attractor states. Progressive reduction in block size, in accord with fracture mechanics, reflects the fact that most blocks were transported at the sole of the ice mass and were subject to the compressive and tensile forces of the ice acting on the stoss surfaces of blocks lying against a bedrock or till surface. The interpretations might apply to a range of homogeneous hard rock lithologies.The Huanghe (Yellow River), one of the largest turbid river systems in the world, has long been recognized as a major contributor of suspended particulate matter (SPM) to the ocean. However, over the last few decades, the SPM export flux of the Huanghe has decreased over 90 % due to the high management, impacting the global export of particulate organic carbon (POC). To better constrain sources and modes of transport of POC beyond the previously investigated transportation of POC near the channel surface, SPM samples were for the first time collected over a whole channel cross-section in the lower Huanghe. Riverine SPM samples were analyzed for particle size and major element contents, as well as for POC content and dual carbon isotopes (13C14C). Clear vertical and lateral heterogeneities of the physical and chemical properties of SPM are observed within the river cross-section. For instance, finer SPM carry more POC in general with higher 14C± 500 radiocarbon years) and dominated by organic carbon (OC) from the biosphere, while the lithospheric fraction is ca. 12 %. The mobilization of aged and refractory OC, including radiocarbon-dead biospheric OC, from deeper soil horizons of the loess–paleosol sequence through erosion in the Chinese Loess Plateau is an important mechanism contributing to fluvial POC in the Huanghe drainage basin. Altogether, anthropogenic activities can drastically change the compositions and transport dynamics of fluvial POC, consequentially altering the feedback of the source-to-sink trajectory of a river system to regional and global carbon cycles.
The grain size 2 mm is the conventional border between sand and gravel. This size is used extensively, and generally without much physical justification, to discriminate between such features as sedimentary deposit type (clast-supported versus matrix-supported), river type (gravel bed versus sand bed), and sediment transport relation (gravel versus sand). Here we inquire as to whether this 2 mm boundary is simply a social construct upon which the research community has decided to agree or whether there is some underlying physics. We use dimensionless arguments to show the following for typical conditions on Earth, i.e., natural clasts (e.g., granitic or limestone) in 20 ∘C water. As grain size ranges from 1 to 5 mm (a narrow band including 2 mm), sediment suspension becomes vanishingly small at normal flood conditions in alluvial rivers. We refer to this range as pea gravel. We further show that bedload movement of a clast in the pea gravel range with, for example, a size of 4 mm moving over a bed of 0.4 mm particles has an enhanced relative mobility compared to a clast with a size of 40 mm moving over a bed of the same 4 mm particles. With this in mind, we use 2 mm here as shorthand for the narrow pea gravel range of 1–5 mm over which transport behavior is distinct from both coarser and finer material. The use of viscosity allows the delineation of a generalized dimensionless bed grain size discriminator between “sand-like” and “gravel-like” rivers. The discriminator is applicable to sediment transport on Titan (ice clasts in flowing methane/ethane liquid at reduced gravity) and Mars (mafic clasts in flowing water at reduced gravity), as well as Earth.
The difficulties of measuring bedload transport in gravel-bed rivers have given rise to the morphological method wherein sediment transport can be inferred from changes in riverbed elevation and estimates of the distance traveled by sediment: its path length. Because current methods for estimating path length are time- and labor-intensive, we present a method to estimate a characteristic path length from repeat digital elevation models (DEMs of difference, i.e., DoDs). We propose an automated method to extract the spacing between erosional and depositional sites on the DoD by the application of variational mode decomposition (VMD), a signal processing method, to quantify the spacing as a proxy for path length. We developed this method using flume experiments where bed topography and sediment flux were measured and then applied it to published field data with physical path length measured from tracer measurements. Our sediment transport estimates were not significantly different than the measured sediment flux at lower discharges in the lab. However, we observed an underestimation of sediment flux at the higher discharges in the flume study. We interpret this as a limit of the method in confined settings, where sediment transport becomes decoupled from morphological changes. We also explore how the time between survey acquisitions, the morphological active width relative to the channel width, and DoD thresholding techniques affect the proposed method and the potential issues they pose for the morphological method in general.
The interplay between climate-forced sea level change, erosional and depositional processes, and flexural isostasy in deep time on passive margin deltas remains poorly understood. We performed a series of conceptual simulations to investigate flexural isostatic responses to high-frequency fluctuations in water and sediment load associated with climatically driven sea level changes. We model a large drainage basin that discharges to a continental margin and produces a large deltaic depocenter, then prescribe synthetic and climatic-driven sea level curves of different frequencies to assess flexural response. Results show that flexural isostatic responses are bidirectional over 100–1000 kyr timescales and are in sync with the magnitude, frequency, and direction of sea level fluctuations and that isostatic adjustments play an important role in driving along-strike and cross-shelf river mouth migration and sediment accumulation. Our findings demonstrate that climate-forced sea level changes produce a feedback mechanism that results in self-sustaining creation of accommodation into which sediment is deposited and plays a major role in delta morphology and stratigraphic architecture.
Chemical weathering of sedimentary rocks can release carbon dioxide (CO2) and consume oxygen (O2) via the oxidation of petrogenic organic carbon and sulfide minerals. These pathways govern Earth's surface system and climate over geological timescales, but the present-day weathering fluxes and their environmental controls are only partly constrained due to a lack of in situ measurements. Here, we investigate the gaseous exchange of CO2O2CO2CO2CO2CO2CO2. However, CO2
Massive sediment pulses in catchments are a key alpine multi-risk component. Substantial sediment redistribution in alpine catchments frequently causes flooding, river erosion, and landsliding and affects infrastructure such as dam reservoirs as well as aquatic ecosystems and water quality. While systematic rock slope failure inventories have been collected in several countries, the subsequent cascading sediment redistribution is virtually unaccessed. For the first time, this contribution reports the massive sediment redistribution triggered by the multi-stage failure of more than 130 000 m3> 600 rockfall events (> 1 m3) from four rock slope catchments with different slope aspects and their volume estimates. In a sediment cascade approach, we evaluated erosion, transport, and deposition from the rock face to the upper channelized erosive debris flow channel, then to the widened dispersive debris flow channel, and finally to the outlet into the braided sediment-supercharged Jochbach river. We observe the decadal flux of more than 400 000 m32–103 mm a−1) to massive erosive regimes (102 mm a−1) within single years, and the massive sediment redistribution (iv) shows limited dependency on rainfall frequency and intensity. This study provides generic information on spatial and temporal patterns of massive sediment pulses in highly sediment-charged alpine catchments.
Wave–current ripples that develop on seabeds of mixed non-cohesive sand and cohesive clay are commonplace in coastal and estuarine environments. While laboratory research on ripples forming in these types of mixed-bed environments is relatively limited, it has identified deep cleaning, the removal of clay below the ripple troughs, as an important factor controlling ripple development. New large-scale flume experiments seek to address this sparsity in data by considering two wave–current conditions with initial clay content, C0, ranging from 0 % to 18.3 %. The experiments record ripple development and pre- and post-experiment bed clay contents to quantify clay winnowing. The present experiments are combined with previous wave-only, wave–current, and current-only experiments to produce a consistent picture of larger and smaller flatter ripples over a range of wave–current conditions and C0. Specifically, the results reveal a sudden decrease in the ripple steepness for C0 > 10.6 %, likely associated with a decrease in hydraulic conductivity of 3 orders of magnitude. Accompanying the sudden change in steepness is a gradual linear decrease in wavelength with C0C0 > 7.4 %. Ultimately, for the highest values of C0, the bed remains flat, but clay winnowing still takes place, albeit at a rate 2 orders of magnitude lower than for rippled beds. For a given flow, the initiation time, when ripples first appear on a flat bed, increases with increasing C0. This, together with the fact that the bed remains flat for the highest values of C0, demonstrates that the threshold of motion increases with C0. The inferred threshold enhancement, and the occurrence of large and small ripples, is used to construct a new three-dimensional phase diagram of bed characteristics involving the wave and current Shields parameters and C0, which has important implications for morphodynamic modelling.
Rock avalanches reach considerably greater runout lengths than predicted by Coulomb friction. While it has been known for a long time that runout length increases with volume, explaining the increase qualitatively is still a challenge. In this study, the widely used Voellmy rheology is reinterpreted and modified. Instead of adding a Coulomb friction term and a velocity-dependent term, the modified rheology assigns the two terms to different regimes of velocity. While assuming a transition between Coulomb friction and flow at a given velocity is the simplest approach, a reinterpretation of an existing model for the kinetic energy of random particle motion predicts a dependence of the crossover velocity on the thickness of the rock avalanche. Analytical solutions for a lumped mass on a simple 1D topography reveal the existence of a slope-dominated and a height-dominated regime within the regime of flow. In the slope-dominated regime, the kinetic energy at the foot of the slope depends mainly on the slope angle, while the absolute height relative to the valley floor has little effect, and vice versa. Both regimes can be distinguished by the ratio of a length scale derived from the rheology and the length scale of the topography. Long runout occurs in the height-dominated regime. In combination with empirical relations between volume, thickness, and height, the approach based on the random kinetic energy model reproduces the scaling of runout length with volume observed in nature very well.
Reliable descriptions of erosion events are foundational to effective frameworks relevant to the fate of tidal landscape evolution. Besides the rhythmic, predictable action of tidal currents, erosion in shallow tidal environments is strongly influenced by the stochastic wave-induced bottom shear stress (BSS), mainly responsible for sediment resuspension on tidal flats. However, the absence of sufficiently long, measured time series of BSS prevents a direct analysis of the combined tide- and wave-driven erosion dynamics and its proper representation in long-term morphodynamic models. Here we test the hypothesis of describing erosion dynamics in shallow tidal environments as a Poisson process by analysing, with the peak-over-threshold theory, the BSS time series computed using a fully coupled, bi-dimensional numerical model. We perform this analysis on the Venice Lagoon, Italy, taking advantage of several historical surveys done in the last 4 centuries, which allow us to investigate the effects of morphological modifications on spatial and temporal erosion patterns. Our analysis suggests that erosion events on intertidal flats can effectively be modelled as a marked Poisson process in different morphological configurations because the interarrival times, durations, and intensities of the over-threshold exceedances are always well described by exponentially distributed random variables. The resulting statistical characterization allows a straightforward computation of morphological indicators, such as the erosion work, and paves the way for a novel synthetic, yet reliable, approach for the long-term morphodynamic modelling of tidal environments.
A proper understanding of sediment resuspension and transport processes is key to the morphodynamics of shallow tidal environments. However, a complete spatial and temporal coverage of suspended sediment concentration (SSC) to describe these processes is hardly available, preventing the effective representation of depositional dynamics in long-term modelling approaches. Through aiming to couple erosion and deposition dynamics in a unique synthetic theoretical framework, we here investigate SSC dynamics, following a similar approach to that adopted for erosion (D'Alpaos et al., 2024). The analysis with the peak-over-threshold theory of SSC time series computed using a fully coupled, bi-dimensional model allows us to identify interarrival times, intensities, and durations of over-threshold events and test the hypothesis of modelling SSC dynamics as a Poisson process. The effects of morphological modifications on spatial and temporal SSC patterns are investigated in the Venice Lagoon, for which several historical configurations in the last 4 centuries are available. Our results show that, similar to erosion events, SSC can be modelled as a marked Poisson process in the intertidal flats for all the analysed morphological lagoon configurations because exponentially distributed random variables describe over-threshold events well. Although erosion and resuspension are intimately intertwined, erosion alone does not suffice to describe also SSC because of the non-local dynamics due to advection and dispersion processes. The statistical characterization of SSC events completes the framework introduced for erosion mechanics, and together, they represent a promising tool to generate synthetic, yet realistic, time series of shear stress and SSC for the long-term modelling of tidal environments.
The Wei River catchment in the southern part of the Chinese Loess Plateau (CLP) is one of the centers of the agricultural revolution in China. The area has experienced intense land use changes since ∼6000 BCE, which makes it an ideal place to study the response of fluvial systems to past anthropogenic land cover change (ALCC). We apply a numerical landscape evolution model that combines the Landlab landscape evolution model with an evapotranspiration model to investigate the direct and indirect effects of ALCC on hydrological and morphological processes in the Wei River catchment since the mid-Holocene. The results show that ALCC has not only led to changes in discharge and sediment load in the catchment but also affected their sensitivity to climate change. When the proportion of agricultural land area exceeded 50 % (around 1000 BCE), the sensitivity of discharge and sediment yield to climate change increased abruptly indicating a regime change in the fluvial catchment. This was associated with a large sediment pulse in the lower reaches. The model simulation results also show a link between human settlement, ALCC and floodplain development: changes in agricultural land use led to downstream sediment accumulation and floodplain development, which in turn resulted in further spatial expansion of agriculture and human settlement.
Large earthquakes can contribute to mountain growth by building topography but also contribute to mass removal from mountain ranges through widespread mass wasting. On annual to decadal or centennial timescales, large earthquakes also have the potential to significantly alter fluvial sediment dynamics if a significant volume of the sediment generated reaches the fluvial network. In this contribution, we focus on the Melamchi–Indrawati and Bhote Koshi rivers in central Nepal, which have both experienced widespread landsliding associated with the 2015 Gorkha (Nepal) earthquake. Using a time series of high-resolution satellite imagery, we have mapped exposed sediment along the rivers from 2012–2021 to identify zones of active channel deposition and document changes over time. Counter to expectations, we show negligible increases in coarse-sediment accumulation along both river corridors since the Gorkha earthquake. However, an extremely high-concentration flow event on 15 June 2021 caused an approximately 4-fold increase in exposed sediment along a 30 km reach of the channel with up to 12 m of channel aggradation in the Melamchi–Indrawati rivers; this event was localised and did not impact the neighbouring Bhote Koshi catchment. Based on published reports, new helicopter-based photography, and satellite data, we demonstrate that this event was sourced from a localised rainfall event between 4500 and 4800 m and that a significant fraction of the sediment was supplied from sources that were unrelated to the landslides generated by the Gorkha earthquake.
Recent studies have shown that hydrophone sensors can monitor bedload flux in rivers by measuring the self-generated noise (SGN) emitted by bedload particles when they impact the riverbed. However, experimental and theoretical studies have shown that the measured SGN depends not only on bedload flux intensity but also the propagation environment, which differs between rivers. Moreover, the SGN can propagate far from the acoustic source and be well measured at distant river positions without bedload transport. It has been shown that this dependency of the measured SGN data on the propagation environment can significantly affect the performance of monitoring bedload flux by hydrophone techniques. In this article, we propose an inversion model to solve the problem of the SGN propagation and integration effect. In this model, we assume that the riverbed acts as SGN source areas with intensity proportional to the local bedload flux. The inversion model locates the SGN sources and calculates their corresponding acoustic power by solving a system of linear algebraic equations, accounting for the actual measured cross-sectional acoustic power (acoustic mapping) and attenuation properties. We tested the model using data from measured bedload SGN profiles (acoustic mapping with a drift boat) and bedload flux profiles (direct sampling with an Elwha sampler) acquired during two field campaigns conducted in 2018 and 2021 on the Giffre river in the French Alps. Results confirm that the bedload flux measured at different verticals on the river cross-section correlates more with the inversed acoustic power than measured acoustic power. Moreover, it was possible to fit data from the two field campaigns with a common curve after inversion, which was not possible with the measured acoustic data. The results of the inversion model, compared to measured data, show the importance of considering the propagation effect when using the hydrophone technique and offer new perspectives for the calibration of bedload flux with SGN in rivers.
Low-intensity but high-frequency coastal flooding, also known as nuisance flooding, can negatively affect low-lying coastal communities with potentially large socioeconomic effects. Partially driven by wave runup, this type of flooding is difficult to predict due to the complexity of the processes involved. Here, we present the results of a probabilistic analysis of flooding events measured on an eroded beach at the Texas coast. A high-resolution time series of the flooded area was obtained from pictures using convolutional neural network (CNN)-based semantic segmentation methods, as described in the first part of this contribution. After defining flooding events using a peak-over-threshold method, we found that their size follows an exponential distribution. Furthermore, consecutive flooding events were uncorrelated at daily timescales but correlated at hourly timescales, as expected from tidal and day–night cycles. Our measurements confirm the broader findings of a recent multi-site investigation of the probabilistic structure of high-water events that used a semi-empirical formulation for wave runup. Indeed, we found a relatively good statistical agreement between our CNN-based empirical flooding data and predictions using total-water-level estimations. As a consequence, our work supports the validity of a relatively simple probabilistic model of high-frequency coastal flooding driven by wave runup that can be used in coastal risk management and landscape evolution models.