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.     

Climate controls on coupled processes of chemical weathering,bioturbation, and sediment transport across hillslopes

Most hillslope studies examining the interplay between climate and earth surface processes tend to be biased towards eroding parts of landscapes. This limitation makes it difficult to assess how entire upland landscapes, which are mosaics of eroding and depositional areas, evolve physio‐chemically as a function of climate. Here we combine new soil geochemical data and published ¹⁰Be‐derived soil production rates to estimate variations in chemical weathering across two eroding‐to‐depositional hillslopes spanning a climate gradient in southeastern Australia. At the warmer and wetter Nunnock River (NR) site, rates of total soil (–3 to –14 g m^‐2 yr^‐1; negative sign indicates mass loss) and saprolite (–18 to –32 g m^‐2 yr^‐1) chemical weathering are uniform across the hillslope transect. Alternatively, the drier hillslope at Frog's Hollow (FH) is characterized by contrasting weathering patterns in eroding soils (–30 to –53 g m^‐2 yr^‐1) vs. depositional soils (+91 g m^‐2 yr^‐1; positive sign indicates mass addition). This difference partly reflects mineral grain size sorting as a result of upslope bioturbation coupled with water‐driven soil erosion, as well as greater vegetative productivity in moister depositional soils. Both of these processes are magnified in the drier climate. The data reveal the importance of linking the erosion–deposition continuum in hillslope weathering studies in order to fully capture the coupled roles of biota and erosion in driving the physical and chemical evolution of hillslopes. Our findings also highlight the potential limitations of applying current weathering models to landscapes where particle‐sorting erosion processes are active. Copyright © 2018 John Wiley & Sons, Ltd.     

Weathering the escarpment: chemical and physical rates and processes,south‐eastern Australia

Differences in chemical weathering extent and character are expected to exist across topographic escarpments due to spatial gradients of climatic and/or tectonic forcing. The passive margin escarpment of south‐eastern Australia has a debated but generally accepted model of propagation in which it retreated (within 40 Ma) to near its current position following rifting between Australia and New Zealand 85–100 Ma before present. We focus on this escarpment to quantify chemical weathering rates and processes and how they may provide insight into scarp evolution and retreat. We compare chemical weathering extents and rates above and below the escarpment using a mass balance approach coupling major and trace element analyses with previous measurements of denudation rates using cosmogenic nuclides (¹⁰Be and ²⁶Al). We find a slight gradient in saprolite chemical weathering rate as a percentage of total weathering rate across the escarpment. The lowlands area, encompassing the region extending from the base of the escarpment to the coast, experiences a greater extent of chemical weathering than the highland region above the escarpment. Percents of denudation attributable to saprolite weathering average 57 ± 6% and 47 ± 7% at low and high sites respectively. Furthermore, the chemical index of alteration (CIA), a ratio of immobile to mobile oxides in granitic material that increases with weathering extent, have corresponding average values of 73·7 ± 3·9 and 65·5 ± 3·4, indicating lower extents of weathering above the escarpment. Finally, we quantify variations in the rates and extent of chemical weathering at the hillslope scale across the escarpment to suggest new insight into how climate differences and hillslope topography help drive landscape evolution, potentially overprinting longer term tectonic forcing. Copyright © 2009 John Wiley & Sons, Ltd.     

The ‘humped’ soil production function: eroding Arnhem Land,Australia

We report erosion rates and processes, determined from in situ‐produced beryllium‐10 (¹⁰Be) and aluminum‐26 (²⁶Al), across a soil‐mantled landscape of Arnhem Land, northern Australia. Soil production rates peak under a soil thickness of about 35 cm and we observe no soil thicknesses between exposed bedrock and this thickness. These results thus quantify a well‐defined ‘humped’ soil‐production function, in contrast to functions reported for other landscapes. We compare this function to a previously reported exponential decline of soil production rates with increasing soil thickness across the passive margin exposed in the Bega Valley, south‐eastern Australia, and found remarkable similarities in rates. The critical difference in this work was that the Arnhem Land landscapes were either bedrock or mantled with soils greater than about 35 cm deep, with peak soil production rates of about 20 m/Ma under 35–40 cm of soil, thus supporting previous theory and modeling results for a humped soil production function. We also show how coupling point‐specific with catchment‐averaged erosion rate measurements lead to a better understanding of landscape denudation. Specifically, we report a nested sampling scheme where we quantify average erosion rates from the first‐order, upland catchments to the main, sixth‐order channel of Tin Camp Creek. The low (～5 m/Ma) rates from the main channel sediments reflect contributions from the slowly eroding stony highlands, while the channels draining our study area reflect local soil production rates (～10 m/Ma off the rocky ridge; ～20 m/Ma from the soil mantled regions). Quantifying such rates and processes help determine spatial variations of soil thickness as well as helping to predict the sustainability of the Earth's soil resource under different erosional regimes. Copyright © 2009 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.     

The spatial variation of weathering and soil depth across a Triassic sandstone outcrop

In this study, we used an archive of borehole logs from the British Geological Survey to collect information on the spatial structure of weathering that extends from the surface to competent bedrock across the Triassic Sherwood Sandstone Group outcrop (750 km²), in the East Midlands, UK. The borehole logs were used to estimate the thickness of the soil (n = 280) and soil and saprolite (S&S) to competent rock (n = 500). The weathering profile of the sandstone consisted of soil (median thickness ～ 1·5 m) overlying a transition zone of compacted and weakly cemented weathered sandstone saprolite over bedrock. Topographic analysis using a NEXTMAP 5 m × 5 m digital elevation model (DEM) revealed no significant relationships between slope properties (relief, flow length, flow accumulation or slope angle) and soil or S&S thickness. A weak, but statistically significant correlation was found between the thickness of the soil and S&S (r_s = 0·25, p < 0·001, n = 192). The variation in soil thickness may be related to changes in current and historic and land‐use, variation in sandstone properties and the influence of glacial/peri‐glacial processes. The thickness of the saprolite was more variable towards the southern part of the study area, where it increased to a maximum 40 m. We hypothesize and provide evidence that the greater weathering thickness is related to the occurrence of increased faulting in this part of the study region, allowing increased access to meteoric waters. A possible source of increased water supply is meltwater from Quaternary ice sheets; the overburden of ice may have increased sub‐glacial pore water pressure, with the fractures and faults acting as a drainage system for the removal of dissolved weathering products. British Geological Survey © NERC 2010     

Potential of in situ-produced cosmogenic nuclides for quantifying strength reduction of bedrock in soil-mantled hillslopes

This study presents a semi-empirical model for quantifying the reduction in the mechanical strength of bedrock beneath actively eroding soil-mantled hillslopes. The strength reduction of bedrock controls the rate of physical disintegration of saprolite, which supplies fresh minerals that are then exposed to intense chemical weathering in soil sections. To determine the values of parameters employed in the model requires knowledge of the denudation rate of the hillslope, the thickness of the soil and saprolite layers, the strength of fresh bedrock, and the threshold strength for physical erosion at the uppermost face of the saprolite. These parameters can be obtained from cosmogenic nuclide analyses for quartz samples from the soil–saprolite boundary and basic field- and laboratory-based investigations. Further testing of the model within a diverse range of climatic, tectonic, and lithologic environments is likely to provide clues to the mechanisms responsible for local and regional variations in the rates of soil production and chemical weathering upon hillslopes.     

Topographic correlations with soil and regolith thickness from shallow‐seismic refraction constraints across upland hillslopes in the Valles Caldera,New Mexico

How rock is weathered physically and chemically into transportable material is a fundamental question in critical‐zone science. In addition, the distribution of weathered material (soil and intact regolith) across upland landscapes exerts a first‐order control on the hydrology of watersheds. In this paper we present the results of six shallow seismic‐refraction surveys in the Redondo Mountain region of the Valles Caldera, New Mexico. The P‐wave velocities corresponding to soil (≤ 0.6 km s^?1) were inferred from a seventh seismic survey where soil‐thickness data were determined by pit excavation. Using multivariable regression, we quantified the relationships among slope gradient, aspect, and topographic wetness index (TWI) on soil and regolith (soil plus intact regolith) thicknesses. Our results show that both soil and regolith thicknesses vary inversely with TWI in all six survey areas while varying directly with slope aspect (i.e. thicker beneath north‐facing slopes) and inversely with slope gradient (i.e. thinner beneath steep slopes) in the majority of the survey areas. An empirical model based on power‐law relationships between regolith thickness and its correlative variables can fit our inferred thicknesses with R² ‐values up to 0.880 for soil and 0.831 for regolith in areas with significant topographic variations. These results further demonstrate the efficacy of shallow seismic refraction for mapping and determining how soil and regolith variations correlate with topography across upland landscapes. Copyright © 2016 John Wiley & Sons, Ltd.     

Evidence for climatic and hillslope‐aspect controls on vadose zone hydrology and implications for saprolite weathering

Through the delivery of water in snowmelt, climate should govern the rate and extent of saprolite formation in snow‐dominated mountain watersheds, yet the mechanisms by which water flows deeply into regolith are largely unexplored. In this study we link rainfall, snow depth, and water content data from both soil and shallow saprolite to document vadose zone dynamics in two montane catchments over 2 years. Measurements of snow pack thickness and soil moisture reveal strong contrasts between north‐ and south‐facing slopes in both the timing of meltwater delivery and the duration of significant soil wetting in the shallow vadose zone. Despite similar magnitudes of snowmelt recharge, north‐facing slopes have higher sustained soil moisture compared to south‐facing slopes. To help interpret these observations, we use a 2D numerical model of vadose zone dynamics to calculate the expected space–time moisture patterns on an idealized hillslope under two wetting scenarios: a single sustained recharge pulse versus a set of short pulses. The model predicts that the duration of the recharge event exerts a stronger control on the depth and residence time of water in the upper unsaturated zone than the magnitude of the recharge event. Model calculations also imply that water should move more slowly through the subsurface and downward water flux should be substantially reduced when water is applied in several pulses rather than in one sustained event. The results suggest that thicker soil and more deeply weathered rock on north‐facing slopes may reflect greater water supply to the deep subsurface. Copyright © 2015 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.     

Distribution of caesium-137 in soils across a hillslope hollow

The distribution of soil ¹³⁷Cs in relation to selected soil and landform properties was studied across a 16 ha hillslope hollow in the Hunter valley, New South Wales, Australia. The hillslope was used as grazing for cattle. Caesium-137 was not significantly related to the amount of sand, silt, or clay, the bulk density, the organic matter content, the slope angle or the relative distance downslope. However, ¹³⁷Cs was significantly related to the thickness of the soil A horizon. Spatial variations in ¹³⁷Cs were compared with topographic units and a six-element hillslope model, but there was little correspondence. It was thought that the effects of microtopography could have masked potential interrelationships between ¹³⁷Cs and broader scale landform parameters.     

Hillslope erosion in a grassland environment: Calibration and evaluation of the SIBERIA landscape evolution model

Field measurement and modelling of soil erosion provides insights into landscape systems as well as the potential for enhanced landscape management. There are a number of field and numerical methods by which soil erosion and deposition can be quantified. Here we examine the capability of the SIBERIA landscape evolution model to quantify short-term erosion and deposition on a well-managed cattle grazing landscape on the east coast of Australia. The model is calibrated by two methods (1) a geomorphological approach using a site digital elevation model (DEM) and soil data and (2) a laboratory-scale flume. The two calibration processes resulted in similar model input parameters and estimated erosion rates of 3.1 t ha⁻¹ year⁻¹ and 4.4 t ha⁻¹ year⁻¹, respectively. These were found to closely match erosion rates estimated using the environmental tracer ¹³⁷Cs (2.7–4.8 t ha⁻¹ year⁻¹). However, erosion and deposition estimated at individual points along the hillslope was not well correlated with ¹³⁷Cs at the same position due to the temporal averaging of the model and microtopography. Sensitivity analysis showed the model was more sensitive to parameterisation than sub-DEM-scale topography. This places confidence in the model's ability to estimate erosion and deposition across an entire hillslope and catchment on decadal time scales. We also highlight the robustness and flexibility of the calibration methods.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Rilling of hillslopes reclaimed before 1977 surface mining law,Dave Johnston Mine,Wyoming

The erosion of hillslopes reclaimed following the surface-mining of coal is an important geomorphologic concern. However, progress in the design of post-mining topography and its hydrologic response has been hampered by the paucity of research focusing on the propensity of rill formation on these disturbed surfaces. This investigation is intended to partially rectify this situation through the development of regression equations for the prediction of rill frequency and magnitude based upon site characteristics. Rill width can be estimated using soil bulk density, hillslope age, and hillslope length with a standard error of 0.152 log₁₀ units at this location. Rill depth can be estimated using soil shear strength, hillslope length, and hillslope age with a standard error of 0.114 log₁₀ units. Rill length can be estimated using hillslope length, percent silt in the soil, hillslope age, Bouyoucos Clay Ratio, and soil bulk density with a standard error of 3.515 m. Rill frequency can be estimated using soil bulk density, Bouyoucos Clay Ratio, soil compaction, soil reaction (pH), and hillslope age with a standard error of 0.241 rills/m. The cross-product of rill width times depth can be estimated using soil shear strength, hillslope length, hillslope age, and soil bulk density with a standard error of 0.260 log₁₀ units. These results must be tested further under various environmental conditions. Nevertheless, prediction of rill formation seems to be a problem capable of solution.     

How long is a hillslope?

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The role of peri‐glacial active layer development in determining soil‐regolith thickness across a Triassic sandstone outcrop in the UK

This paper examines the weathering processes that have combined to produce the distribution of soil‐regolith (SR) thickness across the Triassic Sherwood Sandstone Group outcrop (750 km²) in Nottinghamshire, UK. Archive borehole logs (n = 282) taken across the outcrop showed that SR thickness had mean and median depths of ~1·8 and 1·5 m, respectively. Cores were taken from a forested site to depths ~3 m for geochemical analysis. At this site the SR thickness was ~1·7 m. Analysis of the loss of elements, compared to bedrock using mass balance calculations (τ) showed that all the calcite and gypsum cement had been removed to depths of >3 m. Thus the major difference between the SR and the underlying saprolite was that the former exists as loose sand as opposed to a semi‐durable rock. Scanning electron microscopy (SEM) analysis of core samples suggested that the non‐durable rock or saprolite had greater cementation of clay particles. We propose that the mechanism through which the clay cement (and other interlocking grain bonds) was eased apart was through freeze–thaw processes associated with the summer ‘active layer development (ALD)’ during the last glacial activity in the UK. We tested this theory by developing a Monte Carlo simulation based on a simplified version of the Stefan equation. Current Arctic datasets of air and ground temperatures were obtained to provide reasonable starting conditions for input variables. These were combined with known data for thermal conductivity, bulk density and moisture content of the Sherwood Sandstone regolith. Model predictions (n = 1000) of the distribution of SR thickness accurately reflect the observed distribution thickness from the borehole logs. This is strong evidence that freeze–thaw and ‘ALD’ processes are major factors in determining the thickness of SR across this outcrop. British Geological Survey © NERC 2012     

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.     

Effects of differential hillslope‐scale water retention characteristics on rainfall–runoff response at the Landscape Evolution Observatory

Hillslopes turn precipitation into runoff and thus exert important controls on various Earth system processes. It remains difficult to collect reliable data necessary for understanding and modeling these Earth system processes in real catchments. To overcome this problem, controlled experiments are being conducted at the Landscape Evolution Observatory at Biosphere 2, The University of Arizona. Previous experiments have revealed differences in hydrological response between 2 landscapes within Landscape Evolution Observatory, even though both landscapes were designed to be identical. In an attempt to discover where the observed differences stem from, we use a fully 3‐dimensional hydrological model (CATchment HYdrology) to show the effect of soil water retention characteristics and saturated hydraulic conductivity on the hydrological response of these 2 hillslopes. We also show that soil water retention characteristics can be derived at hillslope scale from experimental observations of soil moisture and matric potential. It is found that differences in soil packing between the 2 landscapes may be responsible for the observed differences in hydrological response. This modeling study also suggests that soil water retention characteristics and saturated hydraulic conductivity have a profound effect on rainfall–runoff processes at hillslope scale and that parametrization of a single hillslope may be a promising step in modeling rainfall–runoff response in real catchments.     

Hillslope and point based soil erosion – an evaluation of a Landscape Evolution Model

Excessive soil erosion and deposition is recognised as a significant land degradation issue. Quantifying soil erosion and deposition is a non-trivial task. One of these methods has been the mathematical modelling of soil erosion and deposition patterns and the processes that drive them. Here we examine the capability of a landscape evolution model to predict both soil erosion rate and pattern of erosion and deposition. This numerical model (SIBERIA) uses a Digital Elevation Model (DEM) to represent the landscape and calculates erosion and deposition at each grid point in the DEM. To assess field soil redistribution rates (SRR) and patterns the distribution of the environmental tracer ¹³⁷Cs has been analysed. Net hill slope SRR predicted by SIBERIA (a soil loss rate of 1.7 to 4.3 t ha^-1 yr^-1) were found to be in good agreement with ¹³⁷Cs based estimates (2.1 – 3.4 t ha^-1 yr^-1) providing confidence in the predictive ability of the model at the hillslope scale. However some differences in predicted erosion/deposition patterns were noted due to historical changes in landscape form (i.e. the addition of a contour bank) and possible causes discussed, as is the finding that soil erosion rates are an order of magnitude higher than likely soil production rates. The finding that SIBERIA can approximate independently quantified erosion and deposition patterns and rates is encouraging, providing confidence in the employment of DEM based models to quantify hillslope erosion rates and demonstrating the potential to upscale for the prediction of whole catchment erosion and deposition. The findings of this study suggest that LEMs can be a reliable alternative to complex and time consuming methods such as that using environmental tracers for the determination of erosion rates. The model and approach demonstrates a new approach to assessing soil erosion that can be employed elsewhere. © 2018 John Wiley & Sons, Ltd.