Post-deforestation soil loss from steepland hillslopes in Taranaki,New Zealand

Soil erosion on steepland hillslopes in Taranaki, New Zealand, where landsliding is the dominant erosion form, was investigated by comparing mean regolith depths between first-order basins that have had their forest cover removed for different periods of time. Regolith depth and slope angle data were collected along 19 profile lines and 30 profile lines from steepland basins that had been deforested for 10 and 85 years, respectively. These profile lines were subdivided into a total of 236 profile segments of relatively linear slope angle and uniform regolith depth, that averaged 17·5 m in length. The depth of pre-existing regolith on post-deforestation landslide sites is estimated from a regression of regolith depth on slope angle for undisturbed (non-landslide) profile segments. Regolith depletion on landslide sites is in turn estimated by subtracting the depth of regolith on landslide sites from the estimate of pre-existing regolith depth. Regolith depletion by post-deforestation landslides, averaged over the entire length of profile lines, gives an estimate of average surface lowering. For the area deforested for 85 years, average surface lowering by post-deforestation landslides is 0·15 ± 0·04 m, and is the same as the difference in mean depth of 0·15 ± 0·11 m between this area and the area deforested for 10 years. Erosion of regolith from hillslopes by processes other than landsliding appears to be minimal. The 0·15 m average surface lowering represents a regolith depletion rate of 1·8 ± 0±5 mm yr^?1. For hillslopes steeper than 28°, where all post-deforestation landslides occur, average surface lowering is 0·20 ± 0·05 m, and the regolith depletion rate is 2±4 · 0±6 mm yr^?1. Average surface lowering is greatest at 0·23 ± 0·07 m on hillslopes steeper than 32° where most post-deforestation landslides occur. Here, the regolith depletion rate is 2·7 ± 0·8 mm yr^?1. A large-magnitude, low-frequency storm in March 1990, produced an average surface lowering of 0·041 m. There were proportionately more landslides in the area deforested for 10 years, illustrating the importance of previous erosion history of hillslopes on the spatial distribution of landslides. There were also comparatively few landslides on steeper hillslopes because previous lower magnitude storms had already removed much of the deeper regolith.     

Critical notch depths for failure of coastal limestone cliffs: case study at Kuro‐shima Island,Okinawa, Japan

Development of a notch at the base of a cliff reduces cliff stability and often induces a collapse. Pleistocene limestone coastal cliffs of elevation 5?m in Kuro‐shima, Ryukyu Islands, have a prominent notch with a depth of 3–4?m at their bases. Around these coastal cliffs, collapses different from previous studies of cliff collapses in the Ryukyu Islands were found; collapses in Kuro‐shima have a horizontal failure surface. The horizontal failure surface, situated at the height of the failure surface corresponding to the retreat point of the notch, is bounded by vertical joints cutting the whole cliff and the reef flat in front of the cliff. Two types of horizontal failure surface were found, triangular and quadrangular; the distinction appears to depend on the angle between the vertical joints and the front face of the cliff. Prior to collapse, these cliffs appear to have been separated from the adjacent cliffs by the development of vertical joints. Consequently, a cliff that will collapse can be identified in advance; cliff instability is strongly dependent on the development of a notch. To study the effect of notch development on cliff collapse, the notch depth at which collapse occurs was calculated using stability analysis. Instability of a cliff increases with notch depth; collapse occurs at the horizontal failure surface when the ratio of the notch depth to the seaward length of the cliff is approximately 0·5–0·7 for a triangular failure surface, and 0·7–0·9 for a quadrangular failure surface. Copyright © 2010 John Wiley & Sons, Ltd.     

Bank‐toe processes in incised channels: the role of apparent cohesion in the entrainment of failed bank materials

Numerous processes may instigate bank retreat and the consequent collection of failed cohesive materials at the bank toe. Cohesion between the failed material and the substrate can provide additional strength to resist direct fluvial entrainment. Failed, cohesive material can act as a form of natural bank‐toe protection by consuming and diverting flow energy that may otherwise be used to further scour the basal zone of incising channels. Investigations in Goodwin Creek, Mississippi, have revealed the existence of apparent cohesion between failed, cohesive blocks and their underlying surface. The method used to assess this cohesion involved a pulley system mounted on a tripod and supporting a load cell. Mean and maximum apparent‐cohesion values of 1·08 kPa and 2·65 kPa, respectively, were measured in this way, identifying a source that bonds blocks to the underlying surface. Cohesion values and types vary spatially and temporally. Tensiometric tests beneath blocks suggest that cohesion resulting from matric suction alone may be as much as 3·5 kPa in summer and 1·8 kPa in winter. Apparent cohesion is believed to have been sufficient to help prevent removal of the largest blocks by a peak flow of 66·4m³/s on 23 September 1997. Maximum excess shear stress required to entrain a D₇₅ block can be augmented by as much as 97% by the presence of apparent cohesion at the block–substrate interface when compared with a condition with zero apparent cohesion at the block underside. Given these findings, it is no longer sufficient to estimate block entrainment in the basal area from block size or bed roughness alone, as in a Shields‐type approach. Copyright © 2001 John Wiley & Sons, Ltd.     

Degradation of glacial deposits quantified with cosmogenic nuclides,Quartermain Mountains,Antarctica

Many glacial deposits in the Quartermain Mountains, Antarctica present two apparent contradictions regarding the degradation of unconsolidated deposits. The glacial deposits are up to millions of years old, yet they have maintained their meter‐scale morphology despite the fact that bedrock and regolith erosion rates in the Quartermain Mountains have been measured at 0·1–4·0 m Ma^?1. Additionally, ground ice persists in some Miocene‐aged soils in the Quartermain Mountains even though modeled and measured sublimation rates of ice in Antarctic soils suggest that without any recharge mechanisms ground ice should sublimate in the upper few meters of soil on the order of 10³ to 10⁵ years. This paper presents results from using the concentration of cosmogenic nuclides beryllium‐10 (¹⁰Be) and aluminum‐26 (²⁶Al) in bulk sediment samples from depth profiles of three glacial deposits in the Quartermain Mountains. The measured nuclide concentrations are lower than expected for the known ages of the deposits, erosion alone does not always explain these concentrations, and deflation of the tills by the sublimation of ice coupled with erosion of the overlying till can explain some of the nuclide concentration profiles. The degradation rates that best match the data range 0·7–12 m Ma^?1 for sublimation of ice with initial debris concentrations ranging 12–45% and erosion of the overlying till at rates of 0·4–1·2 m Ma^?1. Overturning of the tills by cryoturbation, vertical mixing, or soil creep is not indicated by the cosmogenic nuclide profiles, and degradation appears to be limited to within a few centimeters of the surface. Erosion of these tills without vertical mixing may partially explain how some glacial deposits in the Quartermain Mountains maintain their morphology and contain ground ice close to the surface for millions of years. Copyright © 2010 John Wiley & Sons, Ltd.     

Parameterizing a physically based shallow landslide model in a data poor region

Shallow landslides and consequent debris flows are an increasing concern in the Western Ghats of Kerala, India. Their increased frequency has been associated with deforestation and unfavourable land‐use practices in cultivated areas. In order to evaluate the influence of vegetation on shallow slope failures a physically based, dynamic and distributed hydrological model (STARWARS) coupled with a probabilistic slope stability model (PROBSTAB) was applied to the upper Tikovil River basin (55·6 km²). It was tuned with the limited evidence of groundwater conditions during the monsoon season of 2005 and validated against observed landslide activity in the hydrological year 2001–2002. Given the data poor conditions in the region some modifications to the original model were in order, including the estimation of parameters on the basis of generalized information from secondary sources, pedo‐transfer functions, empirical equations and satellite remote sensing data. Despite the poor input, the model captured the general temporal and spatial pattern of instability in the area. Sensitivity analysis proved root cohesion, soil depth and angle of internal friction as the most dominant parameters influencing slope stability. The results indicate the importance of root cohesion in maintaining stability and the critical role of the management of rubber plantations in this. Interception and evapotranspiration showed little influence on the development of failure conditions. The study also highlights the importance of high resolution digital terrain models for the accurate mechanistic prediction of shallow landslide initiation. Copyright © 2009 John Wiley & Sons, Ltd.     

Modified MMF (Morgan–Morgan–Finney) model for evaluating effects of crops and vegetation cover on soil erosion

Modifications are made to the revised Morgan–Morgan–Finney erosion prediction model to enable the effects of vegetation cover to be expressed through measurable plant parameters. Given the potential role of vegetation in controlling water pollution by trapping clay particles in the landscape, changes are also made to the way the model deals with sediment deposition and to allow the model to incorporate particle‐size selectivity in the processes of erosion, transport and deposition. Vegetation effects are described in relation to percentage canopy cover, percentage ground cover, plant height, effective hydrological depth, density of plant stems and stem diameter. Deposition is modelled through a particle fall number, which takes account of particle settling velocity, flow velocity, flow depth and slope length. The detachment, transport and deposition of soil particles are simulated separately for clay, silt and sand. Average linear sensitivity analysis shows that the revised model behaves rationally. For bare soil conditions soil loss predictions are most sensitive to changes in rainfall and soil parameters, but with a vegetation cover plant parameters become more important than soil parameters. Tests with the model using field measurements under a range of slope, soil and crop covers from Bedfordshire and Cambridgeshire, UK, give good predictions of mean annual soil loss. Regression analysis of predicted against observed values yields an intercept value close to zero and a line slope close to 1·0, with a coefficient of efficiency of 0·81 over a range of values from zero to 38·6 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Stratigraphic evidence of episodic talus accumulation on the isle of skye,Scotland

Evidence of episodic development of alluvial talus is given. Three phases of accumulation are identified within the last 3000 years. Average rates of accumulation range from 0·2–1·7 mm yr^?1 during each period. Conventional models relating climate to slope activity are inappropriate due to the local conditions.     

Long term remeasurements of micro-erosion meter rates,Aldabra Atoll,Indian Ocean

Micro-erosion meter sites, installed on Aldabra Atoll, Indian Ocean in 1969 to measure erosion rates in coastal and subaerial sites, were remeasured in 1971 and again in 1982. The orders of magnitude of calculated erosion rates are comparable but the precise figures and patterns differ. For subaerial sites the short term (1969–71) range of 0·11–0·58 mm a^?1 compares with a long term (1971–82) range of 0·06–0·75 mm a^?1; for coastal sites short term ranges of 0·002–7·5 mm a^?1 compare with long term rates of 0·09–2·7 mm a^?1. Interpretation and extrapolation of short term data should only be made with caution.     

Wood distribution in neotropical forested headwater streams of La Selva,Costa Rica

Surveys of wood along 30 forested headwater stream reaches in La Selva Biological Station in north‐eastern Costa Rica represent the first systematic data reported on wood loads in neotropical streams. For streams with drainage areas of 0·1–8·5 km² and gradients of 0·2–8%, wood load ranged from 3 to 34·7 m³ wood/100 m channel and 41–612 m³ wood/ha channel. These values are within the range reported for temperate streams. The variables wood diameter/flow depth, stream power, the presence of backflooding, and channel width/depth are consistently selected as significant predictors by statistical models for wood load. These variables explain half to two‐thirds of the variability in wood load. These results, along with the spatial distribution of wood with respect to the thalweg, suggest that transport processes exert a greater influence on wood loads than recruitment processes. Wood appears to be more geomorphically effective in altering bed elevations in gravel‐bed reaches than in reaches with coarser or finer substrate. Copyright © 2009 John Wiley & Sons, Ltd.     

Alluvial architecture in headwater streams with special emphasis on step–pool topography

Alluvial mountain streams exhibit a range of channel forms: pool–riffle, plane bed, step–pool and cascades. Previous work suggested that these forms exist within discrete, and progressively steeper slope classes. Measurements conducted at over 100 sites in west‐central and central Idaho confirm that slope steepens progressively as one moves from pool–riffle, to plane bed, to step–pool, and finally to cascades. Median slope for pool–riffle topography is 0·0060, for plane beds 0·013, for step–pools 0·044, and for cascades 0·068. There is substantial overlap in the slopes associated with these channel forms. Pool–riffle topography was found at slopes between 0·0010 and 0·015, plane beds between 0·0010 and 0·035, step–pools between 0·015 and 0·134, and cascades between 0·050 and 0·12. Step–pools are particularly striking features in headwater streams. They are characterized by alternating steep and gentle channel segments. The steep segments (step risers) are transverse accumulations of boulder and cobbles, while the gentle segments (pools) contain finer material. Step wavelength is best correlated to step height which is in turn best correlated to the median particle size found on step risers. This result differs from past studies that have reported channel slope to be the dominant control on step wavelength. The presumed geometry and Froude number associated with the features under formative conditions are consistent with the existence field for antidunes and by extension with the hypothesis that step–pools are formed by antidunes. Copyright © 2000 John Wiley & Sons, Ltd.     

A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated downwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Upwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. This article replaces a previously published version (Hydrological Processes, 19 (17), 2915–2929 (2005) [ DOI:10.1002/hyp.5790 ]. See also retraction notice DOI:10.1002/hyp.6350 Copyright © 2006 John Wiley & Sons, Ltd.     

Process–form relationships in Southern Italian badlands: erosion rates and implications for landform evolution

Characteristic badlands are incised into Plio‐Pleistocene clays in Basilicata, southern Italy, creating steep, scarp slopes with knife‐edge ridges (calanchi) and small dome‐shaped forms (biancane). Erosion pin data for the period 1997–2003 give mean annual erosion rates for dome‐shaped biancane in the range 9–19 mm a^?1, while rates for the calanchi scarps are lower, at 7–10 mm a^?1. The erosion pin data also show a non‐linear relationship with slope angle. Maximum erosion rates coincide with a slope angle of 35°, within an envelope defined by combining the theoretical effects of both rainsplash and surface weathering. Monitoring of surface changes and erosion rates for two 0·5 m² cleared swathes on biancane forms reveals a complex relationship between weathering and erosion. Characteristic forms can develop from large blocks of intact clay bedrock over a time period of less than 30 a. The implications of the measured erosion rates for the landform association of mountain front/pediment/domed inselberg are explored. Copyright © 2005 John Wiley & Sons, Ltd.     

Slope stability and valley formation in glacial outwash deposits,North Norfolk

Dry valleys and gullies in fluvioglacial sands and gravels in North Norfolk are fossil forms. Interpretation of their origin demands detailed morphological analysis, but the inevitable ambiguities of form require additional information on the sedimentology of the deposits and associated regolith. Soil mechanical data and stability analysis predict the observed modal slope angles, assuming the water table to be at or near the surface. This suggests that valley formation involved impeded drainage, and slope processes such as shallow slab slides and solifluction over permafrost. The debris was subsequently evacuated from the valleys by meltwater. The stability analysis also explains the different modal slope angles observed at two locations. Slope angle frequency distributions are composite, however, and include data from varying locations up-valley. Up-valley variations in slope angle are seen to reflect variations in the nature of the regolith and consequent changes in soil mechanical properties, as well as changes in valley relief. The up-valley trends in the nature of the regolith are commensurate with a model of headward valley erosion, with a ‘younger’, less weathered regolith at the present valley head.     

Influence of seepage undercutting on the stability of root‐reinforced streambanks

Several mechanisms contribute to streambank failure including fluvial toe undercutting, reduced soil shear strength by increased soil pore‐water pressure, and seepage erosion. Recent research has suggested that seepage erosion of noncohesive soil layers undercutting the banks may play an equivalent role in streambank failure to increased soil pore‐water pressure. However, this past research has primarily been limited to laboratory studies of non‐vegetated banks. The objective of this research was to utilize the Bank Stability and Toe Erosion Model (BSTEM) in order to determine the importance of seepage undercutting relative to bank shear strength, bank angle, soil pore‐water pressure, and root reinforcement. The BSTEM simulated two streambanks: Little Topashaw Creek and Goodwin Creek in northern Mississippi. Simulations included three bank angles (70° to 90°), four pore‐water pressure distributions (unsaturated, two partially saturated cases, and fully saturated), six distances of undercutting (0 to 40 cm), and 13 different vegetation conditions (root cohesions from 0·0 to 15·0 kPa). A relative sensitivity analysis suggested that BSTEM was approximately three to four times more sensitive to water table position than root cohesion or depth of seepage undercutting. Seepage undercutting becomes a prominent bank failure mechanism on unsaturated to partially saturated streambanks with root reinforcement, even with undercutting distances as small as 20 cm. Consideration of seepage undercutting is less important under conditions of partially to fully saturated soil pore‐water conditions. The distance at which instability by undercutting became equivalent to instability by increased soil pore‐water pressure decreased as root reinforcement increased, with values typically ranging between 20 and 40 cm at Little Topashaw Creek and between 20 and 55 cm at Goodwin Creek. This research depicts the baseline conditions at which seepage undercutting of vegetated streambanks needs to be considered for bank stability analyses. Copyright © 2008 John Wiley & Sons, Ltd.     

Retracted and replaced: A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

This article has been retracted and replaced. See Retraction and Replacement Notice DOI: 10.1002/hyp.6350 Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated upwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Downwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. Copyright © 2005 John Wiley & Sons, Ltd.     

Impact–entrainment relationship in a saltating cloud

The problem of impact–entrainment relationship is one of the central issues in understanding saltation, a primary aeolian transport mode. By using particle dynamic analyser measurement technology the movement of saltating particles at the very near‐surface level (1 mm above the bed) was detected. The impacting and entrained particles in the same impact–entrainment process were identified and the speeds, angle with respect to the horizontal, and energy of the impacting and entrained sand cloud were analysed. It was revealed that both the speed and angle of impacting and entrained particles vary widely. The probability distribution of the speed of impacting and entrained particles in the saltating cloud is best described by a Weibull distribution function. The mean impact speed is generally greater than the mean lift‐off speed except for the 0·1–0·2 mm sand whose entrainment is significantly influenced by air drag. Both the impact and lift‐off angles range from 0° to 180°. The mean lift‐off angles range from 39° to 94° while the mean impact angles range from 40° to 78°, much greater than those previously reported. The greater mean lift‐off and especially the mean impact angles are attributed to mid‐air collisions at the very low height, which are difficult to detect by conventional high‐speed photography and are generally ignored in the existing theoretical simulation models. The proportion of backward‐impacting particles also evidences the mid‐air collisions. The impact energy is generally greater than the entrainment energy except for the 0·1–0·2 mm sand. There exists a reasonably good correlation of the mean speed, angle and energy between the impacting and entrained cloud in the impact–entrainment process. The results presented in this paper deserve to be considered in modelling saltation. Copyright © 2002 John Wiley & Sons, Ltd.     

Seismic active earth pressure of cohesive-frictional soil on retaining wall based on a slice analysis method

The M–O (Mononobe–Okabe) theory is used as a standard method to determine the seismic earth pressure. However, the M–O theory does not consider the influence of soil cohesion, and it cannot determine the nonlinear distribution of the seismic earth pressure. This paper presents a general solution for the nonlinear distribution of the seismic active earth pressure of cohesive-frictional soil using the slice analysis method. A new method is proposed to determine the critical failure angle of the backfill wedge under complex conditions, and an iterative calculation method is presented to determine the tension crack depth of the seismic active earth pressure. The considered parameters in the proposed method include the horizontal and vertical seismic coefficients, wall inclination angle, backfill inclination angle, soil friction angle, wall friction angle, soil cohesion, wall adhesion and uniform surcharge. The classical methods of the M–O and Rankine theories can be regarded as special cases of the proposed method. Furthermore, the proposed method is compared with the test results and previously existing solutions to validate the correctness of the results. Additionally, the parameters׳ effect on the critical failure angle, the resultant force, the application-point position, the tension crack depth and the nonlinear distribution of seismic active earth pressure are studied in graphical form.     

Test of statistical means for the extrapolation of soil depth point information using overlays of spatial environmental data and bootstrapping techniques

Hydrological modelling depends highly on the accuracy and uncertainty of model input parameters such as soil properties. Since most of these data are field surveyed, geostatistical techniques such as kriging, classification and regression trees or more sophisticated soil‐landscape models need to be applied to interpolate point information to the area. Most of the existing interpolation techniques require a random or regular distribution of points within the study area but are not adequate to satisfactorily interpolate soil catena or transect data. The soil landscape model presented in this study is predicting soil information from transect or catena point data using a statistical mean (arithmetic, geometric and harmonic mean) to calculate the soil information based on class means of merged spatial explanatory variables. A data set of 226 soil depth measurements covering a range of 0–6·5 m was used to test the model. The point data were sampled along four transects in the Stubbetorp catchment, SE‐Sweden. We overlaid a geomorphology map (8 classes) with digital elevation model‐derived topographic index maps (2–9 classes) to estimate the range of error the model produces with changing sample size and input maps. The accuracy of the soil depth predictions was estimated with the root mean square error (RMSE) based on a testing and training data set. RMSE ranged generally between 0·73 and 0·83 m ± 0·013 m depending on the amount of classes the merged layers had, but were smallest for a map combination with a low number of classes predicted with the harmonic mean (RMSE = 0·46 m). The results show that the prediction accuracy of this method depends on the number of point values in the sample, the value range of the measured attribute and the initial correlations between point values and explanatory variables, but suggests that the model approach is in general scale invariant. Copyright © 2009 John Wiley & Sons, Ltd.     

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.