River channel width controls blocking by slow-moving landslides in California's Franciscan melange

This is the final version. Available from European Geosciences Union via the DOI in this record. To explore the sensitivity of rivers to blocking from landslide debris, we exploit two similar geomorphic settings in California's Franciscan mélange where slow-moving landslides, often referred to as earthflows, impinge on river channels with drainage areas that differ by a factor of 30. Analysis of valley widths and river long profiles over ∼19 km of Alameda Creek (185 km2 drainage area) and Arroyo Hondo (200 km2 drainage area) in central California shows a very consistent picture in which earthflows that intersect these channels force tens of meters of gravel aggradation for kilometers upstream, leading to apparently long-lived sediment storage and channel burial at these sites. In contrast, over a ∼30 km section of the Eel River (5547 km2 drainage area), there are no knickpoints or aggradation upstream of locations where earthflows impinge on its channel. Hydraulic and hydrologic data from United States Geological Survey (USGS) gages on Arroyo Hondo and the Eel River, combined with measured size distributions of boulders input by landslides for both locations, suggest that landslide derived boulders are not mobile at either site during the largest floods (>2-year recurrence) with field-measured flow depths. We therefore argue that boulder transport capacity is an unlikely explanation for the observed difference in sensitivity to landslide inputs. At the same time, we find that earthflow fluxes per unit channel width are nearly identical for Oak Ridge earthflow on Arroyo Hondo, where evidence for blocking is clear, and for the Boulder Creek earthflow on the Eel River, where evidence for blocking is absent. These observations suggest that boulder supply is also an unlikely explanation for the observed morphological differences along the two rivers. Instead, we argue that the dramatically different sensitivity of the two locations to landslide blocking is related to differences in channel width relative to typical seasonal displacements of earthflows. A synthesis of seasonal earthflow displacements in the Franciscan mélange shows that the channel width of the Eel River is ∼5 times larger than the largest annual seasonal displacement. In contrast, during wet winters, earthflows are capable of crossing the entire channel width of Arroyo Hondo and Alameda Creek. In support of this interpretation, satellite imagery shows that immobile earthflow-derived boulders are generally confined to the edges of the channel on the Eel River. By contrast, immobile earthflow-derived boulders jam the entire channel on Arroyo Hondo. Our results imply that lower drainage area reaches of earthflow-dominated catchments may be particularly prone to blocking. By inhibiting the upstream propagation of base-level signals, valley-blocking earthflows may therefore promote the formation of so-called “relict topography”.National Science Foundatio

Pedogenic pathways and deep weathering controls on soil organic carbon in Pacific Northwest forest soils

Characterizing the distribution and dynamics of organic carbon in soil is critical for quantifying changes in the global carbon cycle. In particular, weathering controls on near-surface and deep (>1 m) soil organic carbon (SOC) dynamics have been proposed but limited data prevents us from predicting SOC over topographically complex landscapes and quantifying how changes in climate and perturbations, such as wildfire or land management, influence SOC stocks. To advance our understanding of how weathering alters soil geochemistry and influences SOC storage, we synthesize previous data with a new analysis of the Siuslaw River soil chronosequence from terraces in the Oregon Coast Range, a region that harbors the richest SOC inventories in the continental US. We analyze how the relationships between soil geochemistry, physical properties, and SOC storage vary with weathering status and pathways across soils that span 0.041 to 990 kyr and vary in depth from 1 m to >10 m. To distinguish the key properties and processes influencing SOC storage at different depths, we break our analysis into three depth intervals: 0–30, 30–100, and >100 cm. Our results suggest that the processes that control SOC stocks vary systematically with time and depth owing to weathering impacts on soil properties and pedogenic development. At 30 kyr we observe a peak in SOC stock in the top 100 cm coincident with a peak in oxalate extractable Al and Fe concentrations, representing secondary poorly crystalline minerals, which is consistent with previous studies. We also observe a decline in shallow SOC stock for >30 kyr soils as poorly crystalline minerals are replaced by more stable crystalline forms and soils become clay dominated. At 120 kyr, SOC below 100 cm starts to contribute significantly to the total SOC profile inventory and by 990 kyr, this fraction composes >40% of the total SOC stock. Taken together, our results indicate that total SOC stock increases with soil age as the increased intensity of bedrock weathering deepens the critical zone, creating accommodation space for deep SOC storage. These findings reveal the intimate link between poorly crystalline minerals and SOC and suggest that systematic analysis of soil development in the critical zone provides a first-order constraint on SOC stocks

Quantifying erosion rates and weathering pathways that maximize soil organic carbon storage

Primary minerals that enter soils through bedrock weathering and atmospheric deposition can generate poorly crystalline minerals (PCM) that preferentially associate with soil organic carbon (SOC). These associations hinder microbial decomposition and the release of CO₂ from soils to the atmosphere, making them a critical geochemical control on terrestrial carbon abundance and persistence. Studies that explore these relationships are typically derived from soil chronosequences that experience negligible erosion and thus do not readily translate to eroding landscapes. Here, we propose a theoretical framework to estimate steady-state PCM density and stocks for hilly and mountainous settings by coupling geochemical and geomorphic mass balance equations that account for soil production from bedrock and dust, soil erosion, PCM formation from weathering, and the transformation of PCMs into crystalline phases. We calculate an optimal erosion rate for maximum PCM abundance that arises because PCMs are limited by insufficient weathering at faster erosion rates and loss via “ripening” into more crystalline forms at slower erosion rates. The optimal erosion rate for modeled hilltop soil is modulated by reaction rate constants that govern the efficiency of primary mineral weathering and PCM ripening. By comparing our analysis with global compilations of erosion and soil production rates derived from cosmogenic nuclides, we show that landscapes with slow-to-moderate erosion rates may be optimal for harboring abundant PCM stocks that can facilitate SOC sequestration and limit turnover. Given the growing array of erosion-topography metrics and the widespread availability of high-resolution topographic data, our framework demonstrates how weathering and critical zone processes can be coupled to inform landscape prioritization for persistent SOC storage potential across a broad range of spatial and temporal scales

A discrete random model describing bedrock profile abrasion

We use a simple, collision-based, discrete, random abrasion model to compute the profiles for the stoss faces in a bedrock abrasion process. The model is the discrete equivalent of the generalized version of a classical, collision based model of abrasion. Three control parameters (which describe the average size of the colliding objects, the expected direction of the impacts and the average volume removed from the body due to one collision) are sufficient for realistic predictions. Our computations show the robust emergence of steady state shapes, both the geometry and the time evolution of which shows good quantitative agreement with laboratory experiments.Comment: 9 pages, 6 figure

A simple and efficient GIS tool for volume calculations of submarine landslides

A numeric tool is presented for calculating volumes of topographic voids such as slump scars of landslides, canyons or craters (negative/concave morphology), or alternatively, bumps and hills (positive/convex morphology) by means of digital elevation models embedded within a geographical information system (GIS). In this study, it has been used to calculate landslide volumes. The basic idea is that a (singular) event (landslide, meteorite impact, volcanic eruption) has disturbed an intact surface such that it is still possible to distinguish between the former (undisturbed) landscape and the disturbance (crater, slide scar, debris avalanche). In such cases, it is possible to reconstruct the paleo-surface and to calculate the volume difference between both surfaces, thereby approximating the volume gain or loss caused by the event. I tested the approach using synthetically generated land surfaces that were created on the basis of Shuttle Radar Topography Mission data. Also, I show the application to two real cases, (1) the calculation of the volume of the Masaya Slide, a submarine landslide on the Pacific continental slope of Nicaragua, and (2) the calculation of the void of a segment of the Fish River Canyon, Namibia. The tool is provided as a script file for the free GIS GRASS. It performs with little effort, and offers a range of interpolation parameters. Testing with different sets of interpolation parameters results in a small range of uncertainty. This tool should prove useful in surface studies not exclusively on earth

Biologically inspired simulation of livor mortis

We present a biologically motivated livor mortis simulation that is capable of modelling the colouration changes in skin caused by blood pooling after death. Our approach consists of a simulation of post mortem blood dynamics and a layered skin shader that is controlled by the haemoglobin and oxygen levels in blood. The object is represented by a layered data structure made of a triangle mesh for the skin and a tetrahedral mesh on which the blood dynamics are simulated. This allows us to simulate the skin discolouration caused by livor mortis, including early patchy appearance, fixation of hypostasis and pressure induced blanching. We demonstrate our approach on two different models and scenarios and compare the results to real world livor mortis photographic examples

Pathways of soil genesis in the Coast Range of Oregon, USA

Background and Aims: Soil chronosequences on marine terraces along the Pacific Coast of California and Oregon show evidence of podzolization, though soils ultimately evolve to Ultisols. It is not clear if this pathway of soil evolution can be extended to the humid, inland Oregon Coast Range. Methods: We analyzed soil properties for a fluvial terrace chronosequence sampled along the Siuslaw River (Oregon, USA) about 50 km from the Pacific coast. The seven terraces ranged in age from <3.5 ky to nearly 1,000 ky. Results: There was no evidence of early podsolization. Instead, evidence was found that andisolization starts early and occurs even in older soils when pedogenic iron accumulation and clay synthesis and illuviation dominate. Soils develop the morphology characteristic of Ultisols sometime between 20 and 70 ky, but high levels of oxalate extractable iron and aluminum satisfy criteria of an andic subgroup. Alfisols are not formed as an intermediary stage. Conclusions: The lack of Spodosols inland is due to the inland shift from udic to ustic or xeric moisture regime, which favors summer drying and ripening of short-range order minerals rather than deep leaching or translocation. Other factors are higher pH, different organic chemistry and faster calcium cycling under the Douglas fir inland when compared to the Sitka spruce of the coastal terraces. © 2013 Springer Science+Business Media Dordrecht

How steady are steady-state landscapes? Using visible-near-infrared soil spectroscopy to quantify erosional variability

Although topographic steady state is often used as a simplifying assumption in sediment yield studies and landscape evolution models, the temporal and spatial scales over which this assumption applies in natural landscapes are poorly defined. We used visible–near-infrared (visNIR) spectroscopy to measure the weathering of hilltop soils and quantify local erosional variability in two watersheds in the Oregon Coast Range (United States). One watershed appears adjusted to base-level lowering driven by rock uplift in the Cascadia forearc, while the other is pinned by a gabbroic dike that locally slows river incision and hillslope erosion. Models for uniformly eroding hillslopes imply uniform soil residence times; instead, we observe significant variability around the mean value of 18.8 k.y. (+31.2/–11.8 k.y.) for our adjusted watershed. The magnitude of erosional variability likely reflects the time scales associated with stochastic processes that drive bedrock weathering, soil production, and soil transport (e.g., tree turnover). The residence time distribution for our pinned watershed has a mean value of 72.9 k.y. (+165.6/–50.6 k.y.) and is highly skewed with a substantial fraction of long residence time soils. We speculate that this pattern results from the lithologic control of base level and lateral divide migration driven by erosional contrasts with neighboring catchments. Our novel and inexpensive methodology enables us to quantify for the first time the magnitude of erosional variability in a natural landscape, and thus provides important geomorphic context for studies characterizing regolith development. More generally, we demonstrate that soils can record catchment-scale landscape dynamics that may arise from lithologic controls or forcing due to climate or tectonics

Soil Particle Transport and Mixing Near a Hillslope Crest: 1. Particle Ages and Residence Times

We examine probabilistic elements of soil particle ages and residence times measured from their entry into the mechanically active soil column, focusing on steady conditions near a hillslope crest, where particles steadily move into and through a soil mantle with fixed thickness in the presence of erosion. Our objective is to clarify consequences of the geometry of the particle trajectories and disturbance-driven mixing in relation to the arrangement of the surfaces through which the particles enter and leave the soil element. We introduce an Eulerian-Lagrangian description of particle motions to characterize effects of a nonuniform (two-dimensional) mean motion with superimposed depth-dependent mixing. With weak mixing, the forms and moments of the probability distributions of particle residence times and ages are distinct and are largely controlled by the geometry of the mean particle motion within the soil element. The average particle age is systematically less than the average residence time. These distributions converge to exponential forms with equal means only in the idealized limit of a “well-mixed” system. The analysis provides the foundation for considering particle weathering and chemical losses in relation to residence times and ages and for describing the spatiotemporal structure of tracers, for example, cosmogenic nuclide concentrations and optically stimulated luminescence particle ages, in relation to what these reveal about particle mixing