Journal of Geophysical Research: Earth Surface最新文献

Estimating sediment transport in mountain rivers is challenging because of sediment supply limitation, broad grain size distributions, complex flow hydraulics, and large form drag. Consequently, sediment transport equations are lacking for application in rivers where the bed is coarse and largely immobile, but small fractions of finer, transportable sized material contribute disproportionately to bedload transport. We introduce a framework for estimating sediment transport in mountain rivers that addresses two limitations: estimating the shear stress acting on mobile grains, and accounting for the difference between mobility of size fractions, that is, whether or not a specific grain size can move at a given flow, and transportability, which we define as how much of that size present in the bed will be recruited into transport. We use two bedload data sets to develop equations for predicting incipient motion and transport rates of each grain size fraction present in the bed. We tested the new equations against incipient motion and sediment transport data we collected from streams in the Rocky Mountains, USA, and against published regional sediment yield data. Using this method results in transport estimates where the finer fractions, despite being a small fraction of the bed surface, make up a large part of the total yield. Fractions greater than the median bed grain size are mobile only during peak flood flows, consistent with the existing mountain river bedload data sets. The approach is parsimonious, requiring only data that are often readily available or obtainable: a bed grain size distribution, hydraulic geometry measurements, and discharge.

The motion state of a particle is a crucial aspect of sediment transport problems. In this paper, we conceptualized three states: stillness, “transport”, and “non-transport”, considering that not all the particle motions contribute significantly to the mean sediment transport rate. Starting from a data set of bed-load particle tracks obtained from particle tracking velocimetry, we removed the bias from experimental uncertainty and applied one-dimensional, instantaneous, and non-parametric criteria for distinguishing the states. We described the kinematics of particles in transport and non-transport states, presenting some sample trajectories and the distributions of particle velocity and acceleration. While the transport state presents a clear distinction between stream-wise and transverse particle velocity, the non-transport state is related to isotropic particle jiggling, and does not significantly contribute to the bed-load rate. Vice-versa, the particle motions in the non-transport state are relevant for other summary indicators of the transport process, such as the mean number of moving particles and mean particle velocity. We discuss how applying the proposed non-parametric criterion for state separation is beneficial compared to parameter-dependent alternatives available in the literature. Finally, we provide an outlook on possible applications of our concept for the investigation of other sediment transport processes (incipient motion, solid-fluid interface, creeping flow).

Clastic sediment composition constitutes a key archive of Earth history, controlled by allogenic and autogenic processes that impact weathering, erosion, sediment transfer, and deposition. Deciphering those processes can provide valuable insights into ancient and modern tectonic, geomorphic, climatic, and anthropogenic controls that shape sediment routing systems over a wide range of temporal and spatial scales. However, in order to clearly identify the controls on sediment composition, it is necessary to exclude sources of bias that may mask or diminish the original provenance signal. Such biases may be natural, including mineral fertility, sediment recycling, and grain size, or analytical. This special collection arises from the fifth meeting of the working group on sediment generation held at the University Milano-Bicocca in Milan, Italy, from 28 to 30 June 2022. The collation includes studies that investigate biasing factors affecting all steps of the sediment cascade and all stages of sample collection, preparation, and analysis, as well as case studies that aim to disentangle original provenance signals from geological, environmental, or analytical noise.

Understanding large desert formation/evolution contributing to regional-to-global dust cycles remains a challenge. This study presents the geochemical and Sr-Nd isotope compositions of 51 surface sediment samples collected from the widespread hyper-arid Thar Desert in northwestern India. The major objective is to determine sediment provenance for a better understanding of the formation/evolution mechanism of this Great Indian Desert as well as downwind dust contributions toward the Himalayas. The compositionally immature sandy Thar sediments (CIA ∼50 ± 4, WIP ∼49 ± 12, and Eu_N/Eu* ∼0.80 ± 0.13) are recycled materials derived from the Himalayan orogen and later modified by quartz addition and heavy mineral depletion/sorting processes. The ⁸⁷Sr/⁸⁶Sr (0.7259 ± 0.0012 and ε_Nd (−12.5 ± 2.7) in the bulk of these Thar sediments are different from the earlier published compositions of the eolian sand deposits in northwestern India. The subcategories of Thar materials collected from different dune types exposed over different lithologies (Quaternary alluvium vs. Tertiary and Mesozoic sedimentary formations) are geochemically and isotopically indistinguishable, which indicates their cogenetic sources and/or sediment reworking. Thar sediments collected in this study have a predominant Indus origin along with significant contributions from the upwind Ghaggar-Hakra paleochannels. The Indus sediments are most likely wind-eroded from the shelf region exposed during the low sea stand of LGM and afterward deglaciation. Considering the new and published data sets, the Sr-Nd isotope budget of dust deposited in the Himalayan frontal glaciers indicates that atmospheric mineral dust contribution from the upwind Indo-Gangetic Plain proximal to the Himalayas is at par with dust parcels from distant natural deserts.

Sediment budgets are widely used to measure reach-scale sediment accumulation and evacuation. Such measurements, however, cannot determine when the disturbance is major and the measured sediment mass imbalance is reflective of a river adjusting to a new equilibrium state, as opposed to situations when the disturbance is minor, and the mass imbalance is reflective of a river adjusting within its existing behavioral regime. Sediment sorting among channels and floodplains can have a large effect on how a river responds to a disturbance. Fine sediment may accumulate in the floodplains while coarser sediment erodes from the channel bed. We demonstrate that if a sediment budget does not account for the different behavior and destination of grain sizes, the budget cannot reveal important channel adjustments. In this study, we evaluated how a sand bed river responded to increases in sediment supply by partitioning a sediment budget among silt/clay and five sand fractions. On average, 12 metric tons/meter (downstream)/year of sand was evacuated from the system, but sorting caused channel margins to behave differently from vegetated islands, revealing how a river can slightly narrow while in deficit. Floodplain shaving and bed coarsening evacuated sediment while channel geometry barely changed, consistent with a river adjusting to a minor disturbance within its behavioral regime. This study is an important reminder that sediment mass imbalance does not always lead to channel change. Mechanisms such as floodplain shaving and bed textural change help rivers absorb minor disturbances and resist channel change.

Riverbed elevations play a crucial role in sediment transport and flow resistance, making it essential to understand and quantify their effects. This knowledge is vital for various fields, including river engineering and stream ecology. Previous observations have revealed that fluctuations in the bed surface can exhibit both multifractal and monofractal behaviors. Specifically, the probability distribution function (PDF) of elevation increments may transition from Laplace (two-sided exponential) to Gaussian with increasing scales or consistently remain Gaussian, respectively. These differences at the finest timescale lead to distinct patterns of bedload particle exchange with the bed surface, thereby influencing particle resting times and streamwise transport. In this paper, we utilize the fractional Laplace motion (FLM) model to analyze riverbed elevation series, demonstrating its capability to capture both mono- and multi-fractal behaviors. Our focus is on studying the resting time distribution of bedload particles during downstream transport, with the FLM model primarily parameterized based on the Laplace distribution of increments PDF at the finest timescale. Resting times are extracted from the bed elevation series by identifying pairs of adjacent deposition and entrainment events at the same elevation. We demonstrate that in cases of insufficient data series length, the FLM model robustly estimates the tail exponent of the resting time distribution. Notably, the tail of the exceedance probability distribution of resting times is much heavier for experimental measurements displaying Laplace increments PDF at the finest scale, compared to previous studies observing Gaussian PDF for bed elevation.

Understanding the rates and mechanisms of erosion by subglacial quarrying is a major unsolved problem in geomorphology. Stress enhancement due to load concentration on bedrock ledges between cavities is hypothesized to drive the growth of fractures. Prior work assumed the formation of vertically oriented tensile fractures at the downstream margins of cavities as the controlling process, but did not account for the evolution of the stress field as fractures lengthen, and in particular the dominance of the shearing mode at fracture tips. We used 2D finite element analysis and J-integral methods to analyze stress intensity factors and fracture growth potentials at the tips of preexisting fractures in loaded bedrock steps, taking into account normal and shear components and measured rock strengths. By examining different step heights, step riser angles, rock types, prior fracture locations and orientations, and extents of ice-rock contact zones, we identified some situations favorable for fracture growth, especially in brittle rock types. Typically, however, the growth direction will not be vertically downward but angled up-glacier away from the step riser, a situation unfavorable for quarrying. Moreover, in many situations, the normal stress across fracture planes will be compressive. Non-vertical step risers buttress the bedrock and also suppress fracture growth. In contrast, reducing the sizes of ice-rock contact zones not only increases the loading magnitude, as previously recognized, but also increases intensification of tensile stress at the tips of fractures located just up-glacier. Thus, larger cavities, and hence, fast sliding and low effective pressures, favor quarrying more strongly than previously recognized.

Erosional perturbations from changes in climate or tectonics are recorded in the profiles of bedrock rivers, but these signals can be challenging to unravel in settings with non-uniform lithology. In layered rocks, the surface lithology at a given location varies through time as erosion exposes different layers of rock. Recent modeling studies have used the Stream Power Model (SPM) to highlight complex variations in erosion rates that arise in bedrock rivers incising through layered rocks. However, these studies do not capture the effects of coarse sediment cover on channel evolution. We use the “Stream Power with Alluvium Conservation and Entrainment” (SPACE) model to explore how sediment cover influences landscape evolution and modulates the topographic expression of erodibility contrasts in horizontally layered rocks. We simulate river evolution through alternating layers of hard and soft rock over million-year timescales with a constant and uniform uplift rate. Compared to the SPM, model runs with sediment cover have systematically higher channel steepness values in soft rock layers and lower channel steepness values in hard rock layers. As more sediment accumulates, the contrast in steepness between the two rock types decreases. Effective bedrock erodibilities back-calculated assuming the SPM are strongly influenced by sediment cover. We also find that sediment cover can significantly increase total relief and timescales of adjustment toward landscape-averaged steady-state topography and erosion rates.

Bedload transport can fluctuate considerably over relatively short periods of time and for a given quasi-constant flow rate. What are the implications of replacing the fluctuating signal with a smoothed signal when calculating bedload transport using averaged values, as is common practice? This question was investigated with the BedloadR code, which allows 1D bedload calculation as well as Monte Carlo simulations using a new data set collected in the Severaisse River (French Ecrins massif). Four bedload equations (Camenen & Larson, 2005, https://doi.org/10.1016/j.ecss.2004.10.019; Meyer-Peter & Mueller, 1948; Parker, 1990, https://doi.org/10.1080/00221689009499058; Recking, 2013a, https://doi.org/10.1061/(asce)hy.1943-7900.0000653) were selected for their performance relative to the measured bedload (except for and Meyer-Peter and Mueller) and because each equation has a different mathematical form and degree of nonlinearity. They were used in a Monte Carlo approach, with input probability distributions fitted to the measured river width, slope, bed grain-size distribution, and to the associated (computed) Shields stress. The results show that accounting for natural variability in the calculation reproduces bedload fluctuations well. But overall, when calculating the bedload volume transported by a flow event, accounting for variability systematically leads to higher estimated volumes (of the order of 20%) than those obtained with a deterministic approach using average input parameters. This is a direct consequence of the nonlinearity of the equations.

Permafrost warming has been observed all around the Arctic, however, variations in temperature trends and their drivers remain poorly understood. We present a comprehensive analysis of climatic changes spanning 25 years (1998–2023) at Bayelva (78.92094°N, 11.83333°E) on Spitzbergen, Svalbard. The quality controlled hourly data set includes air temperature, radiation fluxes, snow depth, rainfall, active layer temperature and moisture, and, since 2009, permafrost temperature. Our Bayesian trend analysis reveals an annual air temperature increase of 0.9 ± 0.5°C/decade and strongest warming in September and October. We observed a significant shortening of the snow cover by −14 ± 8 days/decade, coupled with reduced winter snow depth. The active layer simultaneously warmed by 0.6 ± 0.7°C/decade at the top and 0.8 ± 0.5°C/decade at the bottom. While the soil surface got drier, in particular during summer, soil moisture below increased in accordance with the longer unfrozen period and higher winter temperatures. The thawed period prolonged by 10–15 days/decade at different depths. In contrast to earlier top-soil warming, we observed stable temperatures since 2010 and only little permafrost warming (0.14 ± 0.13°C/decade). This is likely due to recently stable winter air temperature and continuously decreasing winter snow depth. This recent development highlights a complex interplay among climate and soil variables. Our distinctive long-term data set underscores (a) the changes in seasonal warming patterns, (b) the influential role of snow cover decline, and (c) that air temperature alone is not a sufficient indicator of change in permafrost environments, thereby highlighting the importance of investigating a wider range of parameters, such as soil moisture and snow characteristics.