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.     

Lateral flow thresholds for aspen forested hillslopes on the Western Boreal Plain,Alberta, Canada

To predict the long‐term sustainability of water resources on the Boreal Plain region of northern Alberta, it is critical to understand when hillslopes generate runoff and connect with surface waters. The sub‐humid climate (P ≤ ET) and deep glacial sediments of this region result in large available soil storage capacity relative to moisture surpluses or deficits, leading to threshold‐dependent rainfall‐runoff relationships. Rainfall simulation experiments were conducted using large magnitude and high intensity applications to examine the thresholds in precipitation and soil moisture that are necessary to generate lateral flow from hillslope runoff plots representative of Luvisolic soils and an aspen canopy. Two adjacent plots (areas of 2·95 and 3·4 m²) of contrasting antecedent moisture conditions were examined; one had tree root uptake excluded for two months to increase soil moisture content, while the second plot allowed tree uptake over the growing season resulting in drier soils. Vertical flow as drainage and soil moisture storage dominated the water balances of both plots. Greater lateral flow occurred from the plot with higher antecedent moisture content. Results indicate that a minimum of 15–20 mm of rainfall is required to generate lateral flow, and only after the soils have been wetted to a depth of 0·75 m (C‐horizon). The depth and intensity of rainfall events that generated runoff > 1 mm have return periods of 25 years or greater and, when combined with the need for wet antecendent conditions, indicate that lateral flow generation on these hillslopes will occur infrequently. Copyright © 2008 John Wiley & Sons, Ltd.     

The influence of climate and microclimate (aspect) on soil creep efficiency: Cinder cone morphology and evolution along the eastern Mediterranean Golan Heights

Although hillslope evolution has been subject to much investigation for more than a century, the effect of climate on the morphology of soil‐mantled hillslopes remains poorly constrained. In this study, we perform numerical simulations of volcanic cinder cones in the Golan Heights (eastern Mediterranean) to estimate soil transport efficiency across a significant north–south gradient in mean annual precipitation (1100 to 500 mm). We use the initial cinder cone morphology (constrained by stratigraphy), the modern hillslope form (surveyed with sub‐meter accuracy) and the eruption age (based on ⁴⁰Ar–³⁹Ar chronology) to predict the best‐fit value of the soil transport coefficient (‘diffusivity’) based on a nonlinear transport model. Our results indicate that the best‐fit diffusivity (K ) varies from 1 to 6 m² ka^?1 among the five cinder cones in our field area. Diffusivity (K ) values vary systematically with precipitation and hillslope aspect; specifically, K is higher on south‐facing (drier) hillslopes and decreases with mean annual precipitation. We interpret this climate dependency to reflect vegetation‐driven variations in apparent soil cohesion, which increases with root network density, and attenuation of rain splash and overland flow erosion, which increases with vegetative ground cover. To assess how vegetative root mass and ground cover vary with precipitation and aspect, we quantified the spatial distribution of NDVI (normalized difference vegetation index) from ASTER satellite images and observed spatial variations that correlate with our calibrated values of K . Analysis of previously studied cinder cones in the USA can be used to extend our framework to arid domains. This endeavor suggests a humped relationship between K and precipitation with maximum diffusivity at mean annual precipitation of 400–600 mm. Copyright © 2017 John Wiley & Sons, Ltd.     

Slope‐area thresholds of road‐induced gully erosion and consequent hillslope–channel interactions

Pikes Peak Highway is a partially paved road between Cascade, Colorado and the summit of Pikes Peak. Significant gully erosion is occurring on the hillslopes due to the concentration of surface runoff, the rearrangement of drainage pathways along the road surface and adjacent drainage ditches, and the high erodibility of weathered Pikes Peak granite that underlies the area. As a result, large quantities of sediment are transported to surrounding valley networks causing significant damage to water quality and aquatic, wetland, and riparian ecosystems. This study establishes the slope/drainage area threshold for gullying along Pikes Peak Highway and a cesium‐137 based sediment budget highlighting rates of gully erosion and subsequent valley deposition for a small headwater basin. The threshold for gullying along the road is S_cr = 0 · 21A^–0·45 and the road surface reduces the critical slope requirement for gullying compared to natural drainages in the area. Total gully volume for the 20 gullies along the road is estimated at 5974 m³, with an erosion rate of 64 m³ yr^–1 to 101 m³ yr^–1. Net valley deposition is estimated at 162 m³ yr^–1 with 120 m³ yr^–1 unaccounted for by gullying. The hillslope–channel interface is decoupled with minimal downstream sediment transport which results in significant local gully‐derived sedimentation. Copyright © 2013 John Wiley & Sons, Ltd.     

Bioturbation on a south‐east Australian hillslope: estimating contributions to soil flux

The contribution of bioturbation to downslope soil transport is significant in many situations, particularly in the context of soil formation, erosion and creep. This study explored the direct flux of soil caused by Aphaenogaster ant mounding, vertebrate scraping and tree‐throw on a wildfire‐affected hillslope in south‐east Australia. This included the development of methods previously applied to Californian gopher bioturbation, and an evaluation of methods for estimating the volume of soil displaced by tree‐throw events. All three bioturbation types resulted in a net downslope flux, but any influence of hillslope angle on flux rates appeared to be overshadowed by environmental controls over the spatial extent of bioturbation. As a result, the highest flux rates occurred on the footslope and lower slope. The overall contribution of vertebrate scraping (57.0 ± 89.4 g m^?1 yr^?1) exceeded that of ant mounding (36.4 ± 66.0 g m^?1 yr^?1), although mean rates were subject to considerable uncertainty. Tree‐throw events, which individually cause major disturbance, were limited in their importance by their scarcity relative to faunalturbation. However, tree‐throw might be the dominant mechanism of biotic soil flux on the mid‐slope provided that it occurs at a frequency of at least 2–3 events ha^?1 yr^?1. Although direct biotic soil flux appears to be geomorphologically significant on this hillslope, such transport processes are probably subordinate to other impacts of bioturbation at this site such as the enhancement of infiltration following wildfire. Copyright © 2011 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.     

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.     

Effects of gradient,distance, curvature and aspect on steep burned and unburned hillslope soil erosion and deposition

Wildfire denudes vegetation and impacts chemical and physical soil properties, which can alter hillslope erosion rates. Post‐wildfire erosion can also contribute disproportionately to long‐term erosion rates and landscape evolution. Post‐fire hillslope erosion rates remain difficult to predict and document at the hillslope scale. Here we use ²¹⁰Pb_aex (lead‐210 mineral‐adsorbed excess) inventories to describe net sediment erosion on steep, convex hillslopes in three basins (unburned, moderately and severely burned) in mountainous central Idaho. We analyzed nearly 300 soil samples for ²¹⁰Pb_aex content with alpha spectrometry and related net sediment erosion to burn severity, aspect, gradient, curvature and distance from ridgetop. We also tested our data against models for advective, linear and non‐linear diffusive erosion. Statistically lower net soil losses on north‐ versus south‐facing unburned hillslopes suggest that greater vegetative cover and soil cohesion on north‐facing slopes decrease erosion. On burned hillslopes, erosion differences between aspects were less apparent and net erosion was more variable, indicating that vegetation influences erosion magnitude and fire drives erosion variability. We estimated net soil losses throughout the length of unburned hillslopes, including through a footslope transition to concave form. In contrast, on burned hillslopes, the subtle shift from convex to concave form was associated with deposition of a post‐fire erosion pulse. Such overall patterns of erosion and deposition are consistent with predictions from a non‐linear diffusion equation. This finding also suggests that concave sections of overall convex hillslopes affect post‐disturbance soil erosion and deposition. Despite these patterns, no strong relationships were evident between local net soil losses and gradient, curvature, distance from ridgetop, or erosion predicted with advection or diffusion equations. The observed relationship between gradient and erosion is therefore likely more complex or stochastic than often described theoretically, especially over relatively short timescales (60–100 years). Copyright © 2016 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.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Changes in carbon and nutrient fluxes from headwaters to ocean in a mountainous temperate to subtropical basin

Water erosion provides major links in global cycles of carbon (C), nitrogen (N) and phosphorus (P). Although significant research on erosion mechanisms has been done, there is still little knowledge on C, N and P fluxes across landscapes to the ocean and their controlling factors in subtropical climates. A four‐year study quantifying and comparing particulate and dissolved C, N and P from multiple scales (microplot, plot, microcatchment, subcatchment, catchment, sub‐basin and basin) was performed in Thukela basin (≈30 000 km²), South Africa. The basin climate was largely subtropical‐humid [mean annual precipitation (MAP) > 980 mm yr^‐1], but temperate (MAP >2000 mm yr^‐1) on the highlands. Open grassland, cropland and bushland were the major land uses. On average, 65, 24 and 4 g m^‐2 yr^‐1 C, N and P were displaced from original topsoil positions, but only 0.33, 0.005 and 0.002 mg m^‐2 yr^‐1 were, respectively, exported to the ocean. The fluxes decreased by 95, 97 and 84%, respectively, from plot to microcatchment outlet; and decreased further in downstream direction by >99% from microcatchment to basin outlet. The hillslope (microplot to microcatchment) fluxes correlated strongly with rainfall parameters. Particulate contributions dominated hillslope fluxes at 73, 81 and 76% of total annual C, N and P, respectively. Although particulate C dominated in the microcatchment‐catchment reach (55%), N (54%) and P (69%) were dominated by dissolved forms. The lower basin zone was dominated by dissolved flux contributions at 93, 81 and 78% for C, N and P for the sub‐basin outlet. These results suggested spatially varying drivers of C, N and P losses from the landscape to the ocean, via the river network. Deposition was envisaged the dominant hillslope level loss process, which gradually gave way to mineralization and biotic uptake in the river network as flux contributions shifted from being predominantly particulate to dissolved forms. Copyright © 2017 John Wiley & Sons, Ltd.     

The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.     

How deep can forest vegetation cover extend their hydrological reinforcing contribution?

An experimental campaign was set up to quantify the contribution of evapotranspiration fluxes on hillslope hydrology and stability for different forest vegetation cover types. Three adjacent hillslopes, respectively, covered by hardwood, softwood, and grass were instrumented with nine access tubes each to monitor soil water dynamics at the three depths of 30, 60, and 100 cm, using a PR2/6 profile probe (Delta‐T Devices Ltd) for about 6 months including wet periods. Soil was drier under softwood and wetter under grass at all the three depths during most of the monitoring period. Matric suction derived via the soil moisture measurements was more responsive to changes in the atmospheric conditions and also recovered faster at the 30 cm depth. Results showed no significant differences between mean matric suction under hardwood (101.6 kPa) with that under either softwood or grass cover. However, a significant difference was found between mean matric suction under softwood (137.5 kPa) and grass (84.3 kPa). Results revealed that, during the wettest period, the hydrological effects from all three vegetation covers were substantial at the 30 cm depth, whereas the contribution from grass cover at 60 cm (2.0 kPa) and 100 cm (1.1 kPa) depths and from hardwood trees at 100 cm depth (1.2 kPa) was negligible. It is surmised that potential instability would have occurred at these larger depths along hillslopes where shallow hillslope failures are most likely to occur in the region. The hydrological effects from softwood trees, 8.1 and 3.9 kPa, were significant as the corresponding factor of safety values showed stable conditions at both depths of 60 and 100 cm, respectively. Therefore, the considerable hydrological reinforcing effects from softwood trees to the 100 cm depth suggest that a hillslope stability analysis would show that hillslopes with softwood trees will be stable even during the wet season.     

Quantifying erosional equilibrium across a slowly eroding,soil mantled landscape

The textbook concept of an equilibrium landscape, which posits that soil production and erosion are balanced and equal channel incision, is rarely quantified for natural systems. In contrast to mountainous, rapidly eroding terrain, low relief and slow-eroding landscapes are poorly studied despite being widespread and densely inhabited. We use three field sites along a climosequence in South Africa to quantify very slow (2-5 m/My) soil production rates that do not vary across hillslopes or with climate. We show these rates to be indistinguishable from spatially invariant catchment-average erosion rates while soil depth and chemical weathering increase strongly with rainfall across our sites. Our analyses imply landscape-scale equilibrium although the dominant means of denudation varies from physical weathering in dry climates to chemical weathering in wet climates. In the two wetter sites, chemical weathering is so significant that clay translocates both vertically in soil columns and horizontally down hillslope catenas, resulting in particle size variation and the accumulation of buried stone lines at the clay-rich depth. We infer hundred-thousand-year residence times of these stone lines and suggest that bioturbation by termites plays a key role in exhuming sediment into the mobile soil layer from significant depths below the clay layer. Our results suggest how tradeoffs in physical and chemical weathering, potentially modulated by biological processes, shape slowly eroding, equilibrium landscapes. © 2019 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Landscape structure,groundwater dynamics,and soil water content influence soil respiration across riparian–hillslope transitions in the Tenderfoot Creek Experimental Forest,Montana

Variability in soil respiration at various spatial and temporal scales has been the focus of much research over the last decade aimed to improve our understanding and parameterization of physical and environmental controls on this flux. However, few studies have assessed the control of landscape position and groundwater table dynamics on the spatiotemporal variability of soil respiration. We investigated growing season soil respiration in a ～393 ha subalpine watershed in Montana across eight riparian–hillslope transitions that differed in slope, upslope accumulated area (UAA), aspect, and groundwater table dynamics. We collected daily‐to‐weekly measurements of soil water content (SWC), soil temperature, soil CO₂ concentrations, surface CO₂ efflux, and groundwater table depth, as well as soil C and N concentrations at 32 locations from June to August 2005. Instantaneous soil surface CO₂ efflux was not significantly different within or among riparian and hillslope zones at monthly timescales. However, cumulative integration of CO₂ efflux during the 83‐day growing season showed that efflux in the wetter riparian zones was ～25% greater than in the adjacent drier hillslopes. Furthermore, greater cumulative growing season efflux occurred in areas with high UAA and gentle slopes, where groundwater tables were higher and more persistent. Our findings reveal the influence of landscape position and groundwater table dynamics on riparian versus hillslope soil CO₂ efflux and the importance of time integration for assessment of soil CO₂ dynamics, which is critical for landscape‐scale simulation and modelling of soil CO₂ efflux in complex landscapes. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessment of long‐term Holocene soil erosion rates in Polish loess areas using sedimentary archives from closed depressions

Soil erosion is a particularly important problem in the loess areas of Central Europe. Numerous studies of past and present soil erosion based on colluvial sediments have so far been conducted. The main problem is the fact that colluvia usually do not represent the complete sedimentation record. Closed depressions (CDs) collect all colluvial sediments from their catchment, therefore, constitute sediment stores enabling the calculation of soil erosion rates. Colluvial sediments and fossil soils, infilling four CDs in the Polish loess belt, were OSL and C‐14 dated. Human settlements near the studied CDs were analyzed. Phases of soil erosion and colluviation from the Neolithic (5400–2900 bc ), from the Middle Bronze Age to the Early Iron Age (1600–0 bc ), and from the Early Middle Ages to Modern Times (500 AD until today) were documented within the CDs studied. Phases of low soil erosion rate and pedogenesis occurred from the Late Vistulian to the Early Neolithic and from the Iron Age to Early Middle Ages. This study reveals that these phases are not synchronous with the soil erosion phases in Central Europe, as the latter mainly occurred in the Bronze Age, Roman Period and Middle Ages. The obtained soil erosion rates were compared with erosion rates in different areas of Central Europe. This study indicates that in loess regions with long‐term agricultural land use, mean erosion rates (i.e. 3.7–5.9 t ha^‐1 yr^‐1) from the Middle Ages to Modern Times were ten times higher than during the entire prehistoric period (0.39–0.67 t ha^‐1 yr^‐1). The mean soil erosion rates for forested CDs was 0.24–0.74 t ha^‐1 yr^‐1. Soil erosion phases are most probably caused by human activities (i.e. land use change) but the early Holocene erosion phase (7.96 +/‐ 0.67 kyr) could have been induced by a climatic fluctuation (e.g. a 8.2 kyr Bond event). Copyright © 2018 John Wiley & Sons, Ltd.     

Asymmetric hillslope erosion following wildfire in Fourmile Canyon,Colorado

Infrequent, high‐magnitude events cause a disproportionate amount of sediment transport on steep hillslopes, but few quantitative data are available that capture these processes. Here we study the influence of wildfire and hillslope aspect on soil erosion in Fourmile Canyon, Colorado. This region experienced the Fourmile Fire of 2010, strong summer convective storms in 2011 and 2012, and extreme flooding in September 2013. We sampled soils shortly after these events and use fallout radionuclides to trace erosion on polar‐ and equatorial‐facing burned slopes and on a polar‐facing unburned slope. Because these radionuclides are concentrated in the upper decimeter of soil, soil inventories are sensitive to erosion by surface runoff. The polar‐facing burned slope had significantly lower cesium‐137 (¹³⁷Cs) and lead‐210 (²¹⁰Pb) inventories (p < 0.05) than either the polar‐facing unburned slope or equatorial‐facing burned slope. Local slope magnitude does not appear to control the erosional response to wildfire, as relatively gently sloping (~20%) polar‐facing positions were severely eroded in the most intensively burned area. Field evidence and soil profile analyses indicate up to 4 cm of local soil erosion on the polar‐facing burned slope, but radionuclide mass balance indicates that much of this was trapped nearby. Using a ¹³⁷Cs‐based erosion model, we find that the burned polar‐facing slope had a net mean sediment loss of 2 mm (~1 kg m^?2) over a one to three year period, which is one to two orders of magnitude higher than longer‐term erosion rates reported for this region. In this part of the Colorado Front Range, strong hillslope asymmetry controls soil moisture and vegetation; polar‐facing slopes support significantly denser pine and fir stands, which fuels more intense wildfires. We conclude that polar‐facing slopes experience the most severe surface erosion following wildfires in this region, indicating that landscape‐scale aridity can control the geomorphic response of hillslopes to wildfires. Copyright © 2018 John Wiley & Sons, Ltd.     

Exploring geochemical controls on weathering and erosion of convex hillslopes: beyond the empirical regolith production function

Landscape curvature evolves in response to physical, chemical, and biological influences that cannot yet be quantified in models. Nonetheless, the simplest models predict the existence of equilibrium hillslope profiles. Here, we develop a model describing steady‐state regolith production caused by mineral dissolution on hillslopes which have attained an equilibrium parabolic profile. When the hillslope lowers at a constant rate, the rate of chemical weathering is highest at the ridgetop where curvature is highest and the ridge develops the thickest regolith. This result derives from inclusion of all the terms in the mathematical definition of curvature. Including these terms shows that the curvature of a parabolic hillslope profile varies with distance from the ridge. The hillslope model (meter‐scale) is similar to models of weathering rind formation (centimeter‐scale) where curvature‐driven solute transport causes development of the thickest rinds at highly curved clast corners. At the clast scale, models fit observations. Here, we similarly explore model predictions of the effect of curvature at the hillslope scale. The hillslope model shows that when erosion rates are small and vertical porefluid infiltration is moderate, the hill weathers at both ridge and valley in the erosive transport‐limited regime. For this regime, the reacting mineral is weathered away before it reaches the land surface: in other words, the model predicts completely developed element‐depth profiles at both ridge and valley. In contrast, when the erosion rate increases or porefluid velocity decreases, denudation occurs in the weathering‐limited regime. In this regime, the reacting mineral does not weather away before it reaches the land surface and simulations predict incompletely developed profiles at both ridge and valley. These predictions are broadly consistent with observations of completely developed element‐depth profiles along hillslopes denuding under erosive transport‐limitation but incompletely developed profiles along hillslopes denuding under weathering limitation in some field settings. Copyright © 2013 John Wiley & Sons, Ltd.