Stochastic processes of soil production and transport: erosion rates,topographic variation and cosmogenic nuclides in the Oregon Coast Range

Landscapes in areas of active uplift and erosion can only remain soil‐mantled if the local production of soil equals or exceeds the local erosion rate. The soil production rate varies with soil depth, hence local variation in soil depth may provide clues about spatial variation in erosion rates. If uplift and the consequent erosion rates are sufficiently uniform in space and time, then there will be tendency toward equilibrium landforms shaped by the erosional processes. Soil mantle thickness would adjust such that soil production matched the erosion. Previous work in the Oregon Coast Range suggested that there may be a tendency locally toward equilibrium between hillslope erosion and sediment yield. Here results from a new methodology based on cosmogenic radionuclide accumulation in bedrock minerals at the base of the soil column are reported. We quantify how soil production varies with soil thickness in the southern Oregon Coast Range and explore further the issue of landscape equilibrium. Apparent soil production is determined to be an inverse exponential function of soil depth, with a maximum inferred production rate of 268 m Ma^?1 occurring under zero soil depth. This rate depends, however, on the degree of weathering of the underlying bedrock. The stochastic and large‐scale nature of soil production by biogenic processes leads to large temporal and spatial variations in soil depth; the spatial variation of soil depth neither supports nor rejects equilibrium morphology. Our observed catchment‐averaged erosion rate of 117 m Ma^?1 is, however, similar to that estimated for the region by others, and to soil production rates under thin and intermediate soils typical for the steep ridges. We suggest that portions of the Oregon Coast Range may be eroding at roughly the same rate, but that local competition between drainage networks and episodic erosional events leads to landforms that are out of equilibrium locally and have a spatially varying soil mantle. Copyright © 2001 John Wiley & Sons, Ltd.     

The sensitivity of hillslope bedrock erosion to precipitation

Decoupling the impacts of climate and tectonics on hillslope erosion rates is a challenging problem. Hillslope erosion rates are well known to respond to changes in hillslope boundary conditions (e.g. channel incision rates) through their dependence on soil thickness, and precipitation is an important control on soil formation. Surprisingly though, compilations of hillslope denudation rates suggest little precipitation sensitivity. To isolate the effects of precipitation and boundary condition, we measured rates of soil production from bedrock and described soils on hillslopes along a semi‐arid to hyperarid precipitation gradient in northern Chile. In each climate zone, hillslopes with contrasting boundary conditions (actively incising channels versus non‐eroding landforms) were studied. Channel incision rates, which ultimately drive hillslope erosion, varied with precipitation rather than tectonic setting throughout the study area. These precipitation‐dependent incision rates are mirrored on the hillslopes, where erosion shifts from relatively fast and biologically‐driven to extremely slow and salt‐driven as precipitation decreases. Contrary to studies in humid regions, bedrock erosion rates increase with precipitation following a power law, from ～1 m Ma^?1 in the hyperarid region to ～40 m Ma^?1 in the semi‐arid region. The effect of boundary condition on soil thickness was observed in all climate zones (thicker soils on hillslopes with stable boundaries compared to hillslopes bounded by active channels), but the difference in bedrock erosion rates between the hillslopes within a climate region (slower erosion rates on hillslopes with stable boundaries) decreased as precipitation decreased. The biotic‐abiotic threshold also marks the precipitation rate below which bedrock erosion rates are no longer a function of soil thickness. Our work shows that hillslope processes become sensitive to precipitation as life disappears and the ability of the landscape to respond to tectonics decreases. Copyright © 2010 John Wiley & Sons, Ltd.     

Very slow erosion rates and landscape preservation across the southwestern slope of the Ladakh Range,India

Erosion rates are key to quantifying the timescales over which different topographic and geomorphic domains develop in mountain landscapes. Geomorphic and terrestrial cosmogenic nuclide (TCN) methods were used to determine erosion rates of the arid, tectonically quiescent Ladakh Range, northern India. Five different geomorphic domains are identified and erosion rates are determined for three of the domains using TCN ¹⁰Be concentrations. Along the range divide between 5600 and 5700 m above sea level (asl), bedrock tors in the periglacial domain are eroding at 5.0 ± 0.5 to 13.1 ± 1.2 meters per million years (m/m.y.)., principally by frost shattering. At lower elevation in the unglaciated domain, erosion rates for tributary catchments vary between 0.8 ± 0.1 and 2.0 ± 0.3 m/m.y. Bedrock along interfluvial ridge crests between 3900 and 5100 m asl that separate these tributary catchments yield erosion rates <0.7 ± 0.1 m/m.y. and the dominant form of bedrock erosion is chemical weathering and grusification. Erosion rates are fastest where glaciers conditioned hillslopes above 5100 m asl by over‐steepening slopes and glacial debris is being evacuated by the fluvial network. For range divide tors, the long‐term duration of the erosion rate is considered to be 40–120 ky. By evaluating measured ¹⁰Be concentrations in tors along a model ¹⁰Be production curve, an average of ~24 cm is lost instantaneously every ~40 ky. Small (<4 km²) unglaciated tributary catchments and their interfluve bedrock have received very little precipitation since ~300 ka and the long‐term duration of their erosion rates is 300–750 ky and >850 ky, respectively. These results highlight the persistence of very slow erosion in different geomorphic domains across the southwestern slope of the Ladakh Range, which on the scale of the orogen records spatial changes in the locus of deformation and the development of an orogenic rain shadow north of the Greater Himalaya. Copyright © 2014 John Wiley & Sons, Ltd.     

Detection of transience in eroding landscapes

Past variations in climate and tectonics have led to spatially and temporally varying erosion rates across many landscapes. In this contribution I examine methods for detecting and quantifying the nature and timing of transience in eroding landscapes. At a single location, cosmogenic nuclides can detect the instantaneous removal of material or acceleration of erosion rates over millennial timescales using paired nuclides. Detection is possible only if one of the nuclides has a significantly shorter half‐life than the other. Currently, the only practical way of doing this is to use cosmogenic in situ carbon‐14 (¹⁴C) alongside a longer lived nuclide, such as beryllium‐10 (¹⁰Be). Hillslope information can complement or be used in lieu of cosmogenic information: in soil mantled landscapes, increased erosion rates can be detected for millennia after the increase by comparing relief and ridgetop curvature. This technique will work as long as the final erosion rate is greater than twice the initial rate. On a landscape scale, transience may be detected based upon disequilibria in channel profiles or ridgetops, but transience can be sensitive to the nature of transient forcing. Where forcing is periodic, landscapes display differing behavior if forcing is driven by changes in base level lowering rates versus changes in the efficiency of either channel or hillslope erosion (e.g. driven by climate change). Oscillations in base level lowering lead to basin averaged erosion rates that reflect a long term average erosion rate despite strong spatial heterogeneity in local erosion rates. This averaging is reflected in ¹⁰Be concentrations in stream sediments. Changes in hillslope sediment transport coefficients can lead to large fluctuations in basin averaged erosion rates, which again are reflected in ¹⁰Be concentrations. The variability of erosion rates in landscapes where both the sediment transport and channel erodibility coefficients vary is dominated by changes to the hillslope transport coefficient. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantifying hillslope erosion rates and processes for a coastal California landscape over varying timescales

The Earth's surface erodes by processes that occur over different spatial and temporal scales. Both continuous, low‐magnitude processes as well as infrequent, high‐magnitude events drive erosion of hilly soil‐mantled landscapes. To determine the potential variability of erosion rates we applied three independent, field‐based methods to a well‐studied catchment in the Marin Headlands of northern California. We present short‐term, basin‐wide erosion rates determined by measuring pond sediment volume (40 years) and measured activities of the fallout nuclides ¹³⁷Cs and ²¹⁰Pb (40–50 years) for comparison with long‐term (>10 ka) rates previously determined from in situ‐produced cosmogenic ¹⁰Be and ²⁶Al analyses. In addition to determining basin‐averaged rates, ¹³⁷Cs and ²¹⁰Pb enable us to calculate point‐specific erosion rates and use these rates to infer dominant erosion processes across the landscape. When examined in the context of established geomorphic transport laws, the correlations between point rates of soil loss from ¹³⁷Cs and ²¹⁰Pb inventories and landscape morphometry (i.e. topographic curvature and upslope drainage area) demonstrate that slope‐driven processes dominate on convex areas while overland flow processes dominate in concave hollows and channels. We show a good agreement in erosion rates determined by three independent methods: equivalent denudation rates of 143 ± 41 m Ma^?1 from pond sediment volume, 136 ± 36 m Ma^?1 from the combination of ¹³⁷Cs and ²¹⁰Pb, and 102 ± 25 m Ma^?1 from ¹⁰Be and ²⁶Al. Such agreement suggests that erosion of this landscape is not dominated by extreme events; rather, the rates and processes observed today are indicative of those operating for at least the past 10 000 years. Copyright © 2006 John Wiley & Sons, Ltd.     

The critical role of climate and saprolite weathering in landscape evolution

Landscapes evolve in response to external forces, such as tectonics and climate, that influence surface processes of erosion and weathering. Internal feedbacks between erosion and weathering also play an integral role in regulating the landscapes response. Our understanding of these internal and external feedbacks is limited to a handful of field‐based studies, only a few of which have explicitly examined saprolite weathering. Here, we report rates of erosion and weathering in saprolite and soil to quantify how climate influences denudation, by focusing on an elevation transect in the western Sierra Nevada Mountains, California. We use an adapted mass balance approach and couple soil‐production rates from the cosmogenic radionuclide (CRN) ¹⁰Be with zirconium concentrations in rock, saprolite and soil. Our approach includes deep saprolite weathering and suggests that previous studies may have underestimated denudation rates across similar landscapes. Along the studied climate gradient, chemical weathering rates peak at middle elevations (1200–2000 m), averaging 112·3 ± 9·7 t km^–2 y^–1 compared to high and low elevation sites (46·8 ± 5·2 t km^?2 y^?1). Measured weathering rates follow similar patterns with climate as those of predicted silica fluxes, modeled using an Arrhenius temperature relationship and a linear relationship between flux and precipitation. Furthermore, chemical weathering and erosion are tightly correlated across our sites, and physical erosion rates increase with both saprolite weathering rates and intensity. Unexpectedly, saprolite and soil weathering intensities are inversely related, such that more weathered saprolites are overlain by weakly weathered soils. These data quantify exciting links between climate, weathering and erosion, and together suggest that climate controls chemical weathering via temperature and moisture control on chemical reaction rates. Our results also suggest that saprolite weathering reduces bedrock coherence, leading to faster rates of soil transport that, in turn, decrease material residence times in the soil column and limit soil weathering. Copyright © 2009 John Wiley & Sons, Ltd.     

Erosion in northwest Tibet from in‐situ‐produced cosmogenic 10Be and 26Al in bedrock

Concentrations of in‐situ‐produced cosmogenic nuclides ¹⁰Be and ²⁶Al in quartz were measured by accelerator mass spectrometry for bedrock basalts and sandstones located in northwest Tibet. The effective exposure ages range between 23 and 134 ka (¹⁰Be) and erosion rates between 4·0 and 24 mm ka^?1. The erosion rates are significantly higher than those in similarly arid Antarctica and Australia, ranging between 0·1 and 1 mm ka^?1, suggesting that precipitation is not the major control of erosion of landforms. Comparison of erosion rates in arid regions with contrasting tectonic activities suggests that tectonic activity plays a more important role in controlling long‐term erosion rates. The obtained erosion rates are, however, significantly lower than the denudation rate of 3000–6000 mm ka^?1 beginning at c. 5‐3 Ma in the nearby Godwin Austen (K2) determined by apatite fission‐track thermochronology. It appears that the difference in erosion rates within different time intervals is indicative of increased tectonic activity at c. 5–3 Ma in northwest Tibet. We explain the low erosion rates determined in this study as reflecting reduced tectonic activity in the last million years. A model of localized thinning of the mantle beneath northwest Tibet may account for the sudden increased tectonic activity at c. 5–3 Ma and the later decrease. Copyright © 2006 John Wiley & Sons, Ltd.     

Coupling chemical weathering with soil production across soil‐mantled landscapes

Soil‐covered upland landscapes constitute a critical part of the habitable world. Our understanding of how they evolve as a function of different climatic, tectonic and geological regimes is important across a wide range of disciplines and depends, in part, on understanding the links between chemical and physical weathering processes. Extensive previous work has shown that soil production rates decrease with increasing soil column thickness, but chemical weathering rates were not measured. Here we examine a granitic, soil‐mantled hillslope at Point Reyes, California, where soil production rates were determined using in situ produced cosmogenic nuclides (¹⁰Be and ²⁶Al), and we quantify the extent as well as the rates of chemical weathering of the saprolite from beneath soil from across the landscape. We collected saprolite samples from the base of soil pits and analysed them for abrasion pH as well as for major and trace elements by X‐ray fluorescence spectroscopy, and for clay mineralogy by X‐ray diffraction spectroscopy. Our results show for the first time that chemical weathering rates decrease with increasing soil thickness and account for 13 to 51 per cent of total denudation. We also show that spatial variation in chemical weathering appears to be topographically controlled: weathering rate decreases with slope across the divergent ridge and increases with upslope contributing area in the convergent swale. Furthermore, to determine the best measure for the extent of saprolite weathering, we compared four different chemical weathering indices – the Vogt ratio, the chemical index of alteration (CIA), Parker's index, and the silicon–aluminium ratio – with saprolite pH. Measurements of the CIA were the most closely correlated with saprolite pH, showing that weathering intensity decreases linearly with an increase in saprolite pH from 4·7 to almost 7. Data presented here are among the first to couple directly rates of soil production and chemical weathering with how topography is likely to control weathering at a hillslope scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Erosion of a high‐altitude,low‐relief area on the Korean Peninsula: implications for its development processes and evolution

The processes involved in the development of high‐altitude, low‐relief areas (HLAs) are still poorly understood. Although cosmogenic nuclides have provided insights into the evolution of HLAs interpreted as paleo‐surfaces, most studies focus on estimating how slowly they erode and thereby their relative stability. To understand actual development processes of HLAs, we applied several techniques of cosmogenic nuclides in the Daegwanryeong Plateau, a well‐known HLA in the Korean Peninsula. Our denudation data from strath terraces, riverine sediments, soils, and tors provide the following conclusions: (1) bedrock incision rate in the plateau (~127 m Myr^?1) is controlled by the incision rate of the western part of the Korean Peninsula, and is similar to the catchment‐wide denudation rate of the plateau (~93 m Myr^?1); (2) the soil production function we observed shows weak depth dependency that may result from highly weathered bedrock coupled with frequent frost action driven by alpine climate; (3) a discrepancy between the soil production and catchment‐wide denudation rates implies morphological disequilibrium in the plateau; (4) the tors once regarded as fossil landforms of the Tertiary do not reflect Tertiary processes; and (5) when compared with those of global paleo‐surfaces (<20 m Myr^?1), our rapid denudation rates suggest that the plateau cannot have maintained its probable initial paleo landscape, and thus is not a paleo‐surface. Our data contribute to understanding the surface processes of actively eroding upland landscapes as well as call into question conventional interpretations of supposed paleo‐surfaces around the world. Copyright © 2015 John Wiley & Sons, Ltd.     

Orders of magnitude increase in soil erosion associated with land use change from native to cultivated vegetation in a Brazilian savannah environment

The Brazilian savanna (cerrado) is a large and important economic and environmental region that is experiencing significant loss of its natural landscapes due to pressures of food and energy production, which in turn has caused large increases in soil erosion. However the magnitude of the soil erosion increases in this region is not well understood, in part because scientific studies of surface runoff and soil erosion are scarce or nonexistent in the cerrado as well as in other savannahs of the world. To understand the effects of deforestation we assessed natural rainfall‐driven rates of runoff and soil erosion on an undisturbed tropical woodland classified as ‘cerrado sensu stricto denso’ and bare soil. Results were evaluated and quantified in the context of the cover and management factor (C‐factor) of the Universal Soil Loss Equation (USLE). Replicated data on precipitation, runoff, and soil loss on plots (5 × 20 m) under undisturbed cerrado and bare soil were collected for 77 erosive storms that occurred over 3 years (2012 through 2014). C‐factor was computed annually using values of rainfall erosivity and soil loss rate. We found an average runoff coefficient of ~20% for the plots under bare soil and less than 1% under undisturbed cerrado. The mean annual soil losses in the plots under bare soil and cerrado were 12.4 t ha^‐1 yr^‐1 and 0.1 t ha^‐1 yr^‐1, respectively. The erosivity‐weighted C‐factor for the undisturbed cerrado was 0.013. Surface runoff, soil loss and C‐factor were greatest in the summer and fall. Our results suggest that shifts in land use from the native to cultivated vegetation result in orders of magnitude increases in soil loss rates. These results provide benchmark values that will be useful to evaluate past and future land use changes using soil erosion models and have significance for undisturbed savanna regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Landscape form and millennial erosion rates in the San Gabriel Mountains,CA

It has been long hypothesized that topography, as well as climate and rock strength, exert first order controls on erosion rates. Here we use detrital cosmogenic ¹⁰Be from 50 basins, ranging in size from 1 to 150 km², to measure millennial erosion rates across the San Gabriel Mountains in southern California, where a strong E–W gradient in relief compared to weak variation in precipitation and lithology allow us to isolate the relationship between topographic form and erosion rate. Our erosion rates range from 35 to 1100 m/Ma, and generally agree with both decadal sediment fluxes and long term exhumation rates inferred from low temperature thermochronometry. Catchment-mean hillslope angle increases with erosion rate until ～ 300 m/Ma, at which point slopes become invariant with erosion rate. Although this sort of relation has been offered as support for non-linear models of soil transport, we use 1-D analytical hillslope profiles derived from existing soil transport laws to show that a model with soil flux linear in slope, but including a slope stability threshold, is indistinguishable from a non-linear law within the scatter of our data. Catchment-mean normalized channel steepness index increases monotonically, though non-linearly, with erosion rate throughout the San Gabriel Mountains, even where catchment-mean hillslope angles have reached a threshold. This non-linearity can be mostly accounted for by a stochastic threshold incision model, though additional factors likely contribute to the observed relationship between channel steepness and erosion rate. These findings substantiate the claim that the normalized channel steepness index is an important topographic metric in active ranges.     

Erosion rates in the southern oregon coast range: Evidence for an equilibrium between hillslope erosion and sediment yield

The relationship of hillslope erosion rates and sediment yield is often poorly defined because of short periods of measurement and inherent spatial and temporal variability in erosion processes. In landscapes containing hillslopes crenulated by alternating topographic noses and hollows, estimates of local hillslope erosion rates averaged over long time periods can be obtained by analysing colluvial deposits in the hollows. Hollows act as local traps for a portion of the colluvium transported down hillslopes, and erosion rates can be calculated using the age and size of the deposits and the size of the contributing source area. Analysis of colluvial deposits in nine Oregon Coast Range hollows has yielded average colluvial transport rates into the hollows of about 35cm³cm^?1yr^?1 and average bedrock lowering rates of about 0.07 mm yr^?1 for the last 4000 to 15000 yr. These rates are consistent with maximum bedrock exfoliation rates of about 0.09 mm yr^?1 calculated from six of the hollows, supporting the interpretation that exfoliation rates limit erosion rates on these slopes. Sediment yield measurements from nine Coast Range streams provide similar basin-wide denudation rates of between 0.05 and 0.08mm yr^?1, suggesting an approximate steady-state between sediment production on hillslopes and sediment yield. In addition, modern sediment yields are similar in basins varying in size from 1 to 1500 km², suggesting that erosion rates are spatially uniform and providing additional evidence for an approximate equilibrium in the landscape.     

The use of ablation‐dominated medial moraines as samplers for 10Be‐derived erosion rates of glacier valley walls,Kichatna Mountains,AK

We use cosmogenic ¹⁰Be concentrations in amalgamated rock samples from active, ice‐cored medial moraines to constrain glacial valley sidewall backwearing rates in the Kichatna Mountains, Alaska Range, Alaska. This dramatic landscape is carved into a small ～65 Ma granitic pluton about 100 km west of Denali, where kilometer‐tall rock walls and ‘cathedral’ spires tower over a radial array of over a dozen valley glaciers. These supraglacial landforms erode primarily by rockfall, but erosion rates are difficult to determine. We use cosmogenic ¹⁰Be to measure rockwall backwearing rates on timescales of 10³–10⁴ years, with a straightforward sampling strategy that exploits ablation‐dominated medial moraines. A medial moraine and its associated englacial debris serve as a conveyor system, bringing supraglacial rockfall debris from accumulation‐zone valley walls to the moraine crest in the ablation zone. We discuss quantitatively several factors that complicate interpretation of cosmogenic concentrations in this material, including the complex scaling of production rates in very steep terrain, the stochastic nature of the rockfall erosion process, the unmixed nature of the moraine sediment, and additional cosmogenic accumulation during transport of the sediment. We sampled medial moraines on each of three glaciers of different sizes and topographic aspects. All three moraines are sourced in areas with identical rock and similar sidewall relief of ～1 km. Each sample was amalgamated from 25 to 35 clasts collected over a 1‐km longitudinal transect of each moraine. Two of the glaciers yield similar ¹⁰Be concentrations (～1·6–2·2 × 10⁴ at/g) and minimum sidewall slope‐normal erosion rates (～0·5–0·7 mm/yr). The lowest ¹⁰Be concentrations (8 × 10³ at/g) and the highest erosion rates (1·3 mm/yr) come from the largest glacier in the range with the lowest late‐summer snowline. These rates are reasonable in an alpine glacial setting, and are much faster than long‐term exhumation rates of the western Alaska Range as determined by thermochronometric studies. Copyright © 2010 John Wiley & Sons, Ltd.     

Geomorphic and cosmogenic nuclide constraints on escarpment evolution in an intraplate setting,Darling Escarpment,Western Australia

The ～900 km long Darling Scarp in Western Australia is one of the most prominent linear topographic features on Earth. Despite the presence of over‐steepened reaches in all westerly flowing streams crossing the scarp, and significant seismic activity within 100 km of the scarp, there is no historical seismicity and no reported evidence for Quaternary tectonic displacements on the underlying Darling Fault. Consequently, it is unclear whether the scarp is a rapidly evolving landform responding to recent tectonic and/or climatic forcing or a more slowly evolving landform. In order to quantify late Quaternary rates of erosion and scarp relief processes, we obtained measurements of the cosmic‐ray produced nuclide beryllium‐10 (¹⁰Be) from outcropping bedrock surfaces along the scarp summit and face, in valley floors, and at stream knickpoints. Erosion rates of bedrock outcrops along the scarp summit surface range from 0·5 to 4·0 m Myr^?1. These are in the same range as erosion rates of 2·1 to 3·6 m Myr^?1 on the scarp face and similar to river incision rates of 2·6 to 11·0 m Myr^?1 from valley floor bedrock straths, indicating that the Darling Scarp is a slowly eroding ‘steady state’ landform, without any significant contemporary relief production over the last several 100 kyr and possibly several million years. Knickpoint retreat rates determined from ¹⁰Be concentrations at the bases of two knickpoints on small streams incised into the scarp are 36 and 46 m Myr^?1. If these erosion rates were sustained over longer timescales, then associated knickpoints may have initiated in the mid‐Tertiary to early Neogene, consistent with early‐mid Tertiary marginal uplift. Ongoing maintenance of stream disequilibrium longitudinal profiles is consistent with slow, regional base level lowering associated with recently proposed continental‐scale tilting, as opposed to differential uplift along discrete faults. Cosmogenic ¹⁰Be analysis provides a useful tool for interpreting the palaeoseismic history of intraplate near‐fault landforms over 10⁵ to 10⁶ years. Copyright © 2010 John Wiley & Sons, Ltd.     

Subglacial abrasion rates at Goldbergkees,Hohe Tauern,Austria, determined from cosmogenic 10Be and 36Cl concentrations

We report concentrations of cosmogenic ¹⁰Be and ³⁶Cl used to determine erosion depths in the recently deglaciated bedrock at Goldbergkees in the Eastern Alps. The glacier covered the sampling sites during the Little Ice Age (LIA) until c. 1940. The youngest ages calculated from these concentrations match the known exposure time after the post‐LIA exposure of <100 years. The apparent age (no cover, no erosion) of most samples, however, is significantly older. We show that the measured nuclide concentrations represent subglacial erosion depths, rather than exposure times. In particular, erosion depths calculated using ¹⁰Be and ³⁶Cl concentrations of individual samples match well, whereas apparent ³⁶Cl ages are consistently older than ¹⁰Be ages. The bedrock at the ‘youngest’ surfaces was deeply eroded (≥ 297 cm) by the Goldbergkees during the late Holocene. In contrast, bedrock at the margin of the LIA ice extent was eroded ≤35 cm. These values convert to subglacial erosion rates on the order of 0.1 mm/a to >5 mm/a. While modeled erosion rates depend on the duration of glacial cover and erosion intrinsic to the different exposure scenarios used for calculation (700–3300 years), modeled total erosion depths are insensitive (5–20% change). Analysis of erosion depths on the transverse valley profile shows a general trend of greatest erosion part way up the valley side and less erosion under thin ice at the lateral margin. A second profile along the valley axis indicates depth of erosion is greatest where the ice abuts the foot of the investigated bedrock riegel and at its lee side just beyond the crest. Copyright © 2016 John Wiley & Sons, Ltd.     

Tracking denudation and sediment production and transport with cosmogenic 10Be in arid,high-altitude Himalayan half-grabens,Zanskar, northern India

Understanding the extent to which local factors, including bedrock and structure, govern catchment denudation in mountainous environments as opposed to broader climate or tectonic patterns provides insight into how landscapes evolve as sediment is generated and transported through them, and whether they have approached steady-state equilibrium. We measured beryllium-10 (¹⁰Be) concentrations in 21 sediment samples from glaciated footwall and hanging wall catchments, including a set of nested catchments, and 12 bedrock samples in the Puga and Tso Morari half-grabens located in the high-elevation, arid Zanskar region of northern India. In the Puga half-graben where catchments are underlain by quartzo-feldspathic gneissic bedrock, bedrock along catchment divides is eroding very slowly, about 5 m/Ma, due to extreme aridity and ¹⁰Be concentrations in catchment sediments are the highest (~60–90 × 10⁵ atoms/g SiO₂) as colluvium accumulates on hillslopes, decoupled from their ephemeral streams. At Puga, ¹⁰Be concentrations and the average erosion rates of a set of six nested catchments demonstrate that catchment denudation is transport-limited as sediment stagnates on lower slopes before reaching the catchment outlet. In the Tso Morari half-graben, gneissic bedrock is also eroding very slowly but ¹⁰Be concentrations in sediments in catchments underlain by low grade meta-sedimentary rocks, are significantly lower (~10–35 × 10⁵ atoms/g SiO₂). In these arid, high-elevation environments, ¹⁰Be concentrations in catchment sediments have more to do with bedrock weathering and transport times than steady-state denudation rates. © 2020 John Wiley & Sons, Ltd.     

Bistable plant–soil dynamics and biogenic controls on the soil production function

Soil formation results from opposite processes of bedrock weathering and erosion, whose balance may be altered by natural events and human activities, resulting in reduced soil depth and function. The impacts of vegetation on soil production and erosion and the feedbacks between soil formation and vegetation growth are only beginning to be explored quantitatively. Since plants require suitable soil environments, disturbed soil states may support less vegetation, leading to a downward spiral of increased erosion and decline in ecosystem function. We explore these feedbacks with a minimal model of the soil–plant system described by two coupled nonlinear differential equations, which include key feedbacks, such as plant‐driven soil production and erosion inhibition. We show that sufficiently strong positive plant–soil feedback can lead to a ‘humped’ soil production function, a necessary condition for soil depth bistability when erosion is assumed to vary monotonically with vegetation biomass. In bistable plant–soil systems, the sustainable soil condition engineered by plants is only accessible above a threshold vegetation biomass and occurs in environments where the high potential rate of erosion exerts a strong control on soil production and erosion. Vegetation removal for agriculture reduces the stabilizing effect of vegetation and lowers the system resilience, thereby increasing the likelihood of transition to a degraded soil state. Copyright © 2016 John Wiley & Sons, Ltd.     

Soil formation and weathering in a permafrost environment of the Swiss Alps: a multi‐parameter and non‐steady‐state approach

Spatially discontinuous permafrost conditions frequently occur in the European Alps. How soils under such conditions have evolved and how they may react to climate warming is largely unknown. This study focuses on the comparison of nearby soils that are characterised by the presence or absence of permafrost (active‐layer thickness: 2–3 m) in the alpine (tundra) and subalpine (forest) range of the Eastern Swiss Alps using a multi‐method (geochemical and mineralogical) approach. Moreover, a new non‐steady‐state concept was applied to determine rates of chemical weathering, soil erosion, soil formation, soil denudation, and soil production. Long‐term chemical weathering rates, soil formation and erosion rates were assessed by using immobile elements, fine‐earth stocks and meteoric ¹⁰Be. In addition, the weathering index (K + Ca)/Ti, the amount of Fe‐ and Al‐oxyhydroxides and clay minerals characteristics were considered. All methods indicated that the differences between permafrost‐affected and non‐permafrost‐affected soils were small. Furthermore, the soils did not uniformly differ in their weathering behaviour. A tendency towards less intense weathering in soils that were affected by permafrost was noted: at most sites, weathering rates, the proportion of oxyhydroxides and the weathering stage of clay minerals were lower in permafrost soils. In part, erosion rates were higher at the permafrost sites and accounted for 79–97% of the denudation rates. In general, soil formation rates (8.8–86.7 t/km²/yr) were in the expected range for Alpine soils. Independent of permafrost conditions, it seems that the local microenvironment (particularly vegetation and subsequently soil organic matter) has strongly influenced denudation rates. As the climate has varied since the beginning of soil evolution, the conditions for soil formation and weathering were not stable over time. Soil evolution in high Alpine settings is complex owing to, among others, spatio‐temporal variations of permafrost conditions and thus climate. This makes predictions of future behaviour very difficult. Copyright © 2016 John Wiley & Sons, Ltd.     

Bedrock erosion and relief production in the northern Flinders Ranges,Australia

Cosmogenic ¹⁰Be concentrations in exposed bedrock surfaces and alluvial sediment in the northern Flinders Ranges reveal surprisingly high erosion rates for a supposedly ancient and stable landscape. Bedrock erosion rates increase with decreasing elevation in the Yudnamutana Catchment, from summit surfaces (13·96 ± 1·29 and 14·38 ± 1·40 m Myr^?1), to hillslopes (17·61 ± 2·21 to 29·24 ± 4·38 m Myr^?1), to valley bottoms (53·19 ± 7·26 to 227·95 ± 21·39 m Myr^?1), indicating late Quaternary increases to topographic relief. Minimum cliff retreat rates (9·30 ± 3·60 to 24·54 ± 8·53 m Myr^?1) indicate that even the most resistant parts of cliff faces have undergone significant late Quaternary erosion. However, erosion rates from visibly weathered and varnished tors protruding from steep bedrock hillslopes (4·17 ± 0·42 to 14·00 ± 1·97 m Myr^?1) indicate that bedrock may locally weather at rates equivalent to, or even slower than, summit surfaces. ¹⁰Be concentrations in contemporary alluvial sediment indicate catchment‐averaged erosion at a rate dominated by more rapid erosion (22·79 ± 2·78 m Myr^?1), consistent with an average rate from individual hillslope point measurements. Late Cenozoic relief production in the Yudnamutana Catchment resulted from (1) tectonic uplift at rates of 30–160 m Myr^?1 due to range‐front reverse faulting, which maintained steep river gradients and uplifted summit surfaces, and (2) climate change, which episodically increased both in situ bedrock weathering rates and frequency–magnitude distributions of large magnitude floods, leading to increased incision rates. These results provide quantitative evidence that the Australian landscape is, in places, considerably more dynamic than commonly perceived. Copyright © 2006 John Wiley & Sons, Ltd.     

Hillslope response to tectonic forcing in threshold landscapes

Hillslopes are thought to poorly record tectonic signals in threshold landscapes. Numerous previous studies of steep landscapes suggest that large changes in long‐term erosion rate lead to little change in mean hillslope angle, measured at coarse resolution. New LiDAR‐derived topography data enables a finer examination of threshold hillslopes. Here we quantify hillslope response to tectonic forcing in a threshold landscape. To do so, we use an extensive cosmogenic beryllium‐10 (¹⁰Be)‐based dataset of catchment‐averaged erosion rates combined with a 500 km² LiDAR‐derived 1 m digital elevation model to exploit a gradient of tectonic forcing and topographic relief in the San Gabriel Mountains, California. We also calibrate a new method of quantifying rock exposure from LiDAR‐derived slope measurements using high‐resolution panoramic photographs. Two distinct trends in hillslope behavior emerge: below catchment‐mean slopes of 30°, modal slopes increase with mean slopes, slope distribution skewness decreases with increasing mean slope, and bedrock exposure is limited; above mean slopes of 30°, our rock exposure index increases strongly with mean slope, and the prevalence of angle‐of‐repose debris wedges keeps modal slopes near 37°, resulting in a positive relationship between slope distribution skewness and mean slope. We find that both mean slopes and rock exposure increase with erosion rate up to 1 mm/a, in contrast to previous work based on coarser topographic data. We also find that as erosion rates increase, the extent of the fluvial network decreases, while colluvial channels extend downstream, keeping the total drainage density similar across the range. Our results reveal important textural details lost in 10 or 30 m resolution digital elevation models of steep landscapes, and highlight the need for process‐based studies of threshold hillslopes and colluvial channels. Copyright © 2012 John Wiley & Sons, Ltd.