Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

Abstract. A proposed null hypothesis for fluvial terrace formation is that internally generated or autogenic processes, such as lateral migration and river-bend cutoff, produce variabilities in channel incision that lead to the abandonment of floodplain segments as terraces. Alternatively, fluvial terraces have the potential to record past environmental changes from external forcings that include temporal changes in sea level and hydroclimate. Terraces in the Trinity River valley have been previously characterized as Deweyville groups and interpreted to record episodic cut and fill during late Pleistocene sea level variations. Our study uses high-resolution topography of a bare-earth digital elevation model derived from airborne lidar surveys along ∼ 88 linear kilometers of the modern river valley. We measure both differences in terrace elevations and widths of paleo-channels preserved on these terraces in order to have two independent constraints on terrace formation mechanisms. For 52 distinct terraces, we quantify whether terrace elevations fit distinct planes – expected for allogenic terrace formation tied to punctuated sea level and/or hydroclimate change – by comparing variability in a grouped set of Deweyville terrace elevations against variability associated with randomly selected terrace sets. Results show Deweyville groups record an initial valley floor abandoning driven by allogenic forcing, which transitions into autogenic forcing for the formation of younger terraces. For these different terrace sets, the slope amongst different terraces stays constant. For 79 paleo-channel segments preserved on these terraces, we connected observed changes in paleo-channel widths to estimates for river paleo-hydrology over time. Our measurements suggest the discharge of the Trinity River increased systematically by a factor of ∼ 2 during the late Pleistocene. Despite this evidence of increased discharge, the similar down-valley slopes between terrace sets indicate that there were likely no increases in sediment-to-water discharge ratios that could be linked to allogenic terrace formation. This is consistent with our elevation clustering analysis that suggests younger terraces are indistinguishable in their elevation variance from autogenic terrace formation mechanisms, even if the changing paleo-channel dimensions might, viewed in isolation, provide a mechanism for allogenic terrace formation. Methods introduced here combine river-reach-scale observations of terrace sets and paleo-hydrology with local observations of terraces and paleo-channels to show how interpretations of allogenic versus autogenic terrace formation can be evaluated within a single river system.

Abstract. We examine the influence of incision thresholds on topographic and scaling properties of landscapes that follow a landscape evolution model (LEM) with terms for stream-power incision, linear diffusion, and uniform uplift. Our analysis uses three main tools. First, we examine the graphical behavior of theoretical relationships between curvature and the steepness index (which depends on drainage area and slope). These relationships plot as straight lines for the case of steady-state landscapes that follow the LEM. These lines have slopes and intercepts that provide estimates of landscape characteristic scales. Such lines can be viewed as counterparts of slope–area relationships, which follow power laws in detachment-limited landscapes but not in landscapes with diffusion. We illustrate the response of these curvature–steepness index lines to changes in the values of parameters. Second, we define a Péclet number that quantifies the competition between incision and diffusion, while taking the incision threshold into account. We examine how this Péclet number captures the influence of the incision threshold on the degree of landscape dissection. Third, we characterize the influence of the incision threshold using a ratio between it and the steepness index. This ratio is a dimensionless number in the case of the LEM that we use and reflects the fraction by which the incision rate is reduced due to the incision threshold; in this way, it quantifies the relative influence of the incision threshold across a landscape. These three tools can be used together to graphically illustrate how topography and process competition respond to incision thresholds.

Abstract. Boulder bars are a common form of riverbed morphology that could be affected by landslide dams. However, few studies have focused on the formation and geometry characteristics of boulder bars due to outburst floods triggered by landslide dam failure. In such a way, eight group landslide dam failure experiments with a movable bed length of 4 to 7 times the dam length with 25 boulder bars were carried out. In addition, 38 boulder bars formed in the field triggered by four landslide dam failures were investigated. The aim of this paper is to study the formation and geometry characteristics of boulder bars along the riverbeds. The results show that boulder bars are formed after peak discharge of outburst flow. The number of boulder bars is 0.4 to 1.0 times the ratio of riverbed length to dam bottom length. Besides, boulder bars have the characteristic of lengthening upstream during the failure process. A boulder bar's upstream edge has a more extensive development than a boulder bar's downstream edge. The length of a boulder bar along the channel changes faster than the boulder bar's width and height. After the dam failure, the boulder bar's length is about 8 to 14 times its width. The relationship between the ratio of boulder bar length to width and the boulder bar's dimensionless length could be described with a hyperbolic equation. The dimensionless area of the boulder bar increases linearly with the dimensionless area of the river section, and the linear ratio is about 0.5. With the field data, this demonstrates that the formation and geometry characteristics of boulder bars in tests are consistent with the field boulder bars. Therefore, the results in this paper are credible and can be applied to the riverbed's geomorphological characteristics analysis triggered by overtopped landslide dam failure. The plentiful experimental and field data could contribute to the community boulder bar research.

Abstract. The detachment of rock fragments from fractured bedrock on hillslopes creates sediment with an initial size distribution that sets the upper limits on particle size for all subsequent stages in the evolution of sediment in landscapes. We hypothesize that the initial size distribution should depend on the size distribution of latent sediment (i.e., fracture-bound blocks in unweathered bedrock) and weathering of blocks both before and during detachment (e.g., disintegration along crystal grain boundaries). However, the initial size distribution is difficult to measure because the interface across which sediment is produced is often shielded from view by overlying soil. Here we overcome this limitation by comparing fracture spacings measured from exposed bedrock on cliff faces with particle size distributions in adjacent talus deposits at 15 talus–cliff pairs spanning a wide range of climates and lithologies in California. Median fracture spacing and particle size vary by more than 10-fold and correlate strongly with lithology. Fracture spacing and talus size distributions are also closely correlated in central tendency, spread, and shape, with b-axis diameters showing the closest correspondence with fracture spacing at most sites. This suggests that weathering has not modified latent sediment either before or during detachment from the cliff face. In addition, talus at our sites has not undergone much weathering after deposition and is slightly coarser than the latent sizes because it contains unexploited fractures inherited from bedrock. We introduce a new conceptual framework for understanding the relative importance of latent size and weathering in setting initial sediment size distributions in mountain landscapes. In this framework, hillslopes exist on a spectrum defined by the ratio of two characteristic timescales: the residence time in saprolite and weathered bedrock and the time required to detach a particle of a characteristic size. At one end of the spectrum, where weathering residence times are negligible, the latent size distribution can be used to predict the initial size distribution. At the other end of the spectrum, where weathering residence times are long, the latent size distribution can be erased by weathering in the critical zone.

Abstract. In mountain ranges, earthquakes can trigger widespread landsliding and mobilize large amounts of organic carbon by eroding soil and vegetation from hillslopes. Following a major earthquake, the landslide-mobilized organic carbon can be exported from river catchments by physical sediment transport processes or stored within the landscape where it may be degraded by heterotrophic respiration. The competition between these physical and biogeochemical processes governs a net transfer of carbon between the atmosphere and sedimentary organic matter, yet their relative importance following a large landslide-triggering earthquake remains poorly constrained. Here, we propose a model framework to quantify the post-seismic redistribution of soil-derived organic carbon. The approach combines predictions based on empirical observations of co-seismic sediment mobilization with a description of the physical and biogeochemical processes involved after an earthquake. Earthquake-triggered landslide populations are generated by randomly sampling a landslide area distribution, a proportion of which is initially connected to the fluvial network. Initially disconnected landslide deposits are transported downslope and connected to rivers at a constant velocity in the post-seismic period. Disconnected landslide deposits lose organic carbon by heterotrophic oxidation, while connected deposits lose organic carbon synchronously by both oxidation and river export. The modeling approach is numerically efficient and allows us to explore a large range of parameter values that exert a control on the fate of organic carbon in the upland erosional system. We explore the role of the climatic context (in terms of mean annual runoff and runoff variability) and rates of organic matter degradation using single pool and multi-pool models. Our results highlight the fact that the redistribution of organic carbon is strongly controlled by the annual runoff and the extent of landslide connection, but less so by the choice of organic matter degradation model. In the context of mountain ranges typical of the southwestern Pacific region, we find that model configurations allow more than 90 % of the landslide-mobilized carbon to be exported from mountain catchments. A simulation of earthquake cycles suggests efficient transfer of organic carbon out of a mountain range during the first decade of the post-seismic period. Pulsed erosion of organic matter by earthquake-triggered landslides is therefore an effective process to promote carbon sequestration in sedimentary deposits over thousands of years.

Abstract. Climate has been viewed as a primary control on the rates and patterns of glacial erosion, yet our understanding of the mechanisms by which climate influences glacial erosion is limited. We hypothesize that climate controls the patterns of glacial erosion by altering the basal thermal regime of glaciers. The basal thermal regime is a first-order control on the spatial patterns of glacial erosion. Polythermal glaciers contain both cold-based portions that protect bedrock from erosion and warm-based portions that actively erode bedrock. In this study, we model the impact of various climatic conditions on glacier basal thermal regimes and patterns of glacial erosion in mountainous regions. We couple a sliding-dependent glacial erosion model with the Parallel Ice Sheet Model (PISM) to simulate the evolution of the glacier basal thermal regime and glacial erosion in a synthetic landscape. We find that both basal thermal regimes and glacial erosion patterns are sensitive to climatic conditions, and glacial erosion patterns follow the patterns of the basal thermal regime. Cold temperature leads to limited glacial erosion at high elevations due to cold-based conditions. Increasing precipitation can overcome the impact of cold temperature on the basal thermal regime by accumulating thick ice and lowering the melting point of ice at the base of glaciers. High precipitation rates, therefore, tend to cause warm-based conditions at high elevations, resulting in intensive erosion near the peak of the mountain range. Previous studies often assessed the impact of climate on the spatial patterns of glacial erosion by integrating climatic conditions into the equilibrium line altitudes (ELAs) of glaciers, and glacial erosion is suggested to be maximal around the ELA. However, our results show that different climatic conditions produce glaciers with similar ELAs but different patterns of basal thermal regime and glacial erosion, suggesting that there might not be any direct correlation between ELAs and glacial erosion patterns.

Abstract. Surface flow on rilled hillslopes tends to produce sediment yields that scale nonlinearly with total hillslope length. The widespread observation lacks a single unifying theory for such a nonlinear relationship. We explore the contribution of rill network geometry to the observed yield–length scaling relationship. Relying on an idealized network geometry, we formally develop probability functions for geometric variables of contributing area and rill length. In doing so, we contribute towards a complete probabilistic foundation for the Hack distribution. Using deterministic and empirical functions, we then extend the probability theory to the hydraulic variables that are related to sediment detachment and transport. A Monte Carlo simulation samples hydraulic variables from hillslopes of different lengths to provide estimates of sediment yield. The results of this analysis demonstrate a nonlinear yield–length relationship as a result of the rill network geometry. Theory is supported by numerical modeling, wherein surface flow is routed over an idealized numerical surface and a natural surface from northern Arizona. Numerical flow routing demonstrates probability functions that resemble the theoretical ones. This work provides a unique application of the Scheidegger network to hillslope settings which, because of their finite lengths, result in unique probability functions. We have addressed sediment yields on rilled slopes and have contributed towards understanding Hack's law from a probabilistic reasoning.

Abstract. Wave-transported boulders represent important records of storm and tsunami impact over geological timescales. Their use for hazard assessment requires chronological information on their displacement that in many cases cannot be achieved by established dating approaches. To fill this gap, this study investigated, for the first time, the potential of optically stimulated luminescence rock surface exposure dating (OSL-RSED) for estimating cliff-detachment ages of wave-transported coastal boulders. The approach was tested on calcarenite clasts at the Rabat coast, Morocco. Calibration of the OSL-RSED model was based on samples with rock surfaces exposed to sunlight for ∼ 2 years, and OSL exposure ages were evaluated against age control deduced from satellite images. Our results show that the dating precision is limited for all targeted boulders due to the local source rock lithology which has low amounts of quartz and feldspar. The dating accuracy may be affected by erosion rates on boulder surfaces of 0.02–0.18 mm yr−1. Nevertheless, we propose a robust relative chronology for boulders that are not affected by significant post-depositional erosion and that share surface angles of inclination with the calibration samples. The relative chronology indicates that (i) most boulders were detached from the cliff by storm waves; (ii) these storms lifted boulders with masses of up to ∼ 24 t; and (iii) the role of storms in the formation of boulder deposits along the Rabat coast is more significant than previously assumed. Although OSL-RSED cannot provide reliable absolute exposure ages for the coastal boulders in this study, the approach has large potential for boulder deposits composed of rocks with larger amounts of quartz or feldspar and less susceptibility to erosion.

Abstract. The presence of bare patches within otherwise vegetated coastal marshes is sometimes considered to be a symptom of marsh dieback and the subsequent loss of important ecosystem services. Here we studied the topographical conditions determining the presence and revegetation of bare patches in three marsh sites with contrasting tidal range, sediment supply, and plant species: the Scheldt estuary (the Netherlands), Venice lagoon (Italy), and Blackwater marshes (Maryland, USA). Based on GIS (geographic information system) analyses of aerial photos and lidar imagery of high resolution (≤2×2 m pixels), we analyzed the topographic conditions under which bare patches occur, including their surface elevation, size, distance from channels, and whether they are connected or not to channels. Our results demonstrate that, for the different marsh sites, bare patches can be connected or unconnected to the channel network and that there is a positive relationship between the width of the connecting channels and the size of the bare patches, in each of the three marsh sites. Further, pixels located in bare patches connected to channels occur most frequently at the lowest elevations and farthest distance from the channels. Pixels in bare patches disconnected from channels occur most frequently at intermediate elevations and distances from channels, and vegetated marshes dominate at highest elevations and shortest distances from channels. In line with previous studies, revegetation in bare patches is observed in only one site with the highest tidal range and highest sediment availability, and it preferentially occurs from the edges of small unconnected bare patches at intermediate elevations and intermediate distances from channels. Although our study is only for three different marsh sites with large variations in local conditions, such as tidal range, sediment availability, and plant species, it suggests that similar topographic conditions determine the occurrence of bare patches. Such insights may inform decision makers on coastal marsh management on where to focus monitoring of early signatures of marsh degradation.