Columbia Mountain landslide: late-glacial emplacement and indications of future failure, Northwestern Montana, USA

The well preserved and undissected Columbia Mountain landslide, which is undergoing suburban development, was studied to estimate the timing and processes of emplacement. The landslide moved westward from a bedrock interfluve of the northern Swan Range in Montana, USA onto the deglaciated floor of the Flathead Valley. The landslide covers an area of about 2 km², has a toe-to-crown height of 1100 m, a total length of 3430 m, a thickness of between 3 and 75 m, and an approximate volume of 40 million m³. Deposits and landforms define three portions of the landslide; from the toe to the head they are: (i) clast-rich diamictons made up of gravel-sized angular rock fragments with arcuate transverse ridges at the surface; (ii) silty and sandy deposits resting on diamictons in an internally drained depression behind the ridges; and (iii) diamictons containing angular and subangular pebble-to block-sized clasts (some of which are glacially striated) in an area of lumpy topography between the depression and the head of the landslide. Drilling data suggest the diamictons cover block-to-slab-sized bedrock clasts that resulted from an initial stage of the failure.The landslide moved along a surface that developed at a high angle to the NE-dipping, thinly bedded metasediments of the Proterozoic Belt Supergroup. The exposed slope of the main scarp dips 30–37°W. A hypothetical initial rotational failure of the lower part of a bedrock interfluve may have transported bedrock clasts into the valley. The morphology and deposits at the surface of the landslide indicate deposition by a rock avalanche (sturzstrom) derived from a second stage of failure along the upper part of the scarp.The toe of the Columbia Mountain landslide is convex-west in planview, except where it was deflected around areas now occupied by glacial kettles on the north and south margins. Landsliding, therefore, occurred during deglaciation of the valley while ice still filled the present-day kettles. Available chronostratigraphy suggests that the ˜1-km thick glacier in the region melted before 12,000 ¹⁴C years BP—within 3000 years of the last glacial maximum. Deglaciation and hillslope failure are likely causally linked. Failure of the faceted interfluve was likely due tensile fracturing of bedrock along a bedding-normal joint set shortly after glacial retreat from the hillslope.Open surficial tension fractures and grabens in the Swan Range are limited to an area above the crown of the landslide. Movement across these features suggests that extensional flow of bedrock (sackung) is occurring in what remains of the ridge that failed in the Columbia Mountain landslide. The fractures and grabens likely were initiated during failure, but their morphologies suggest active extension across some grabens. Continued movement of bedrock above the crown may result in future mass movements from above the previous landslide scarp. Landslides sourced from bedrock above the scarp of the late-glacial Columbia Mountain landslide, which could potentially be triggered by earthquakes, are geologic hazards in the region.     

Evolution of undercut slopes on abandoned incised meanders in the Eastern Highland Rim of Tennessee, USA

Field surveys, location-for-time reasoning, and computer modeling were used to study the evolution of slopes on valley walls of abandoned bedrock meanders on the Eastern Highland Rim, Tennessee. Hillslopes on the undercut slopes of cutoff incised meanders were ordered as to relative age by the height of their meander floors above the modern stream level. The assumption is that the undercut slope is actively eroded by the stream until abandonment of the meander, at which time the slope begins to evolve to a different form. More-advanced stages of evolution occur on walls of higher meanders that were abandoned earlier. The most rapid change in this initial form is the elimination of a free face, which occurs soon after the meander is abandoned. In addition, the hillslopes associated with even the lowest (youngest) cutoff meanders show somewhat gentler overall gradients than the actively undercut slopes. Hillslopes associated with meanders 3 to about 20 m above modern stream level maintain straight segments with angles showing only a slight decrease from the 36–38° associated with the lowest cutoffs; overall angles decrease, however, as the straight segment becomes shorter. The oldest slopes, those on cutoffs 30 m or greater above modern stream level, have developed into convex–concave slopes with maximum slopes of 15°.A hillslope evolution model based on previously published algorithms was used to simulate the transition of actively undercut hillslopes into hillslopes on abandoned meanders. Hillslope modeling is particularly useful in this setting. If the valley incision rate is known, an age can be estimated for the cutoff and hence for the hillslope. Alternatively, if hillslope process rates are known, a model age obtained for the hillslope can be used to estimate an incision rate. Even where both incision rates and hillslope process rates are poorly constrained, as in the present setting, modeling allows assumptions about specific rates to be evaluated by determining their implications for other rates. For example, for three cutoff meanders along one stream, best-fit criteria were used to select process rates for the model. Model ages of hillslopes were then obtained and compared with those calculated from a valley-incision rate measured elsewhere in the same physiographic province. For two of the hillslopes, model ages were found to be much younger than those calculated from the incision rate. In order to make the two ages agree, unreasonably low process rates had to be used in the model, thus implying that the incision rate probably underestimates the actual incision rate in this valley.Experimentation with heights of initial profiles, again using best-fit criteria, suggests that since abandonment of the highest cutoff, the plateau has been downwasting at a rate about one-fourth that of the valley incision rate, a finding in agreement with published rates of chemical denudation in the area.     

Impact of hillslope-derived sediment supply on drainage basin development in small watersheds at the northern border of the central Alps of Switzerland

This paper explores the effects of hillslope mobility on the evolution of a 10-km² drainage basin located at the northern border of the Swiss Alps. It uses geomorphologic maps and the results of numerical models that are based on the shear stress formulation for fluvial erosion and linear diffusion for hillslope processes. The geomorphic data suggest the presence of landscapes with specific cross-sectional geometries reflecting variations in the relationships between processes in channels and on hillslopes. In the headwaters, the landscape displays parabolic cross-sectional geometries indicating that mass delivered to channels by hillslope processes is efficiently removed. In the trunk stream portion, the landscape is (i) V-shaped if the downslope flux of mass is balanced by erosion in channels (i.e. if mass delivered to channels by hillslope processes is efficiently removed) and (ii) U-shaped if in-channel accumulation of hillslope-derived material occurs. This latter situation indicates a non-balanced mass flux between processes in channels and on hillslopes.Information about the spatial pattern of the postglacial depth of erosion allows comparative estimates to be made about the erosional efficiency for the various landscapes that were mapped in the study area. The data suggest that the erosional potential and sediment discharge are reduced for the situation of a non-balanced mass flux between processes in channels and on hillslopes. These findings are also supported by the numerical model. Indeed, the model results show that high hillslope mobility tends to reduce the hillslope relief and to inhibit dissection and formation of channels. In contrast, stable hillslopes tend to promote fluvial incision, and the hillslope relief increases. The model results also show that very low erosional resistance of bedrock promotes backward erosion and steepening of channel profiles in headwaters. Beyond that, the model reveals that sediment discharge generally increases with decreasing erosional resistance of bedrock, but that this increase decays exponentially with increasing magnitudes of fluvial and hillslope mobilities. Very high hillslope diffusivities even tend to reduce the erosional potential of the whole watershed. It appears that besides rates of base-level lowering, factors limiting sediment discharge might be the nonlinear relationships between processes in channels and on hillslopes.     

Late Pleistocene tectonic-geomorphological development within a passive margin — The Cariatá trough, northeastern Brazil

Studies on the geomorphological evolution of the South American passive margin have been based on the pediplanation model, which predicts that its morphology is a response to regional uniform uplift and concomitant development of erosion surfaces. We combined remote sensing, geological mapping, lithostratigraphic and facies analyses, and luminescence dating in the Cariatá trough, northeastern Brazil, in order to determine how brittle tectonics and climate influenced colluviation and the shaping of local landforms in the Quaternary. Our work indicates that Cariatá is an asymmetrical trough  40 km long,  25 km wide, 250–550 m deep, and delimited by ENE–WSW-trending faults to the north and south. We recognized an ENE–WSW-oriented compression related to a strike-slip faulting regime, which corresponds to the present-day stress field in the region. This faulting event led to the deposition of colluvial fans, shed from adjacent uplifted crustal blocks, in a tectonically controlled depression. The colluvial succession is  45 m thick and presents two facies assemblages that filled the southern and, in particular, the northern borders of the trough: non-cohesive debrisflow and mudflow deposits. Optically stimulated luminescence dates of the sedimentary infill yielded ages at 224–128 ka and 45–28 ka, dominated by debrisflow and mudflow deposits, respectively. These ages may be over-estimated due to poor bleaching of colluvium, but they and our field data suggest that the margins of the trough were tectonically uplifted and eroded twice in the Late Pleistocene. The spasmodic colluvial accretion reflects the occurrence of high-magnitude, low-recurrence episodes probably associated with climate shifts in a semi-arid hillslope system. It follows that the present-day low-lying piedmont in which the Cariatá trough occurs is a juxtaposition of surfaces of various ages. This trough may have counterparts across the region. These conclusions do not corroborate the application of the cyclical pediplanation hypothesis in the area.     

Initiation conditions for debris flows generated by runoff at Chalk Cliffs, central Colorado

We have monitored initiation conditions for six debris flows between May 2004 and July 2006 in a 0.3 km² drainage basin at Chalk Cliffs; a band of hydrothermally-altered quartz monzonite in central Colorado. Debris flows were initiated by water runoff from colluvium and bedrock that entrained sediment from rills and channels with slopes ranging from about 14° to 45°. The availability of channel material is essentially unlimited because of thick channel fill and refilling following debris flows by rock fall and dry ravel processes. Rainfall exceeding I = 6.61(D)^− 0.77, where I is rainfall intensity (mm/h), and D is duration (h), was required for the initiation of debris flows in the drainage basin. The approximate minimum runoff discharge from the surface of bedrock required to initiate debris flows in the channels was 0.15 m³/s. Colluvium in the basin was unsaturated immediately prior to (antecedent) and during debris flows. Antecedent, volumetric moisture levels in colluvium at depths of 1 cm and 29 cm ranged from 4–9%, and 4–7%, respectively. During debris flows, peak moisture levels in colluvium at depths of 1 cm and 29 cm ranged from 10–20%, and 4–12%, respectively. Channel sediment at a depth of 45 cm was unsaturated before and during debris flows; antecedent moisture ranged from 20–22%, and peak moisture ranged from 24–38%. Although we have no measurements from shallow rill or channel sediment, we infer that it was unsaturated before debris flows, and saturated by surface-water runoff during debris flows.Our results allow us to make the following general statements with regard to debris flows generated by runoff in semi-arid to arid mountainous regions: 1) high antecedent moisture levels in hillslope and channel sediment are not required for the initiation of debris flows by runoff, 2) locations of entrainment of sediment by successive runoff events can vary within a basin as a function of variations in the thickness of existing channel fill and the rate of replenishment of channel fill by rock fall and dry ravel processes following debris flows, and 3) rainfall and simulated surface-water discharge thresholds can be useful in understanding and predicting debris flows generated by runoff and sediment entrainment.     

Clast Fabric in Relict Periglacial Colluvium, Salamanca Re-Entrant, Southwestern New York, USA

Samples of macrofabric data obtained from colluvial deposits near the late-Wisconsinan glacial margin display widely divergent characteristics. Those from highly weathered and jointed sandstone plot as girdle distributions. Samples from red colluvium of early-Wisconsinan age form moderately strong clusters. Fabrics from shallow depth in medium-textured brown colluvium of late-Wisconsinan age generally have modes aligned with the local slope but are weak, and appear to have been modified by frost heaving. Clasts in the basal layers of an overlying loess unit were emplaced by frost heaving from the brown colluvium and have very weak or uniform fabrics. Particle shape exerts an influence on fabric strength in the brown colluvium, with samples composed of more elongated clasts displaying stronger fabrics. Fabrics from near-surface deposits show generally weaker shape and strength characteristics and higher inclinations than those at depth. Calculated seasonal frost penetration is consistent with the depth to which weaker, more steeply inclined fabrics are found. Sampling in relict periglacial deposits should be performed only below the level of disturbance by seasonal frost.     

Post-wildfire erosion response in two geologic terrains in the western USA

Volumes of eroded sediment after wildfires vary substantially throughout different geologic terrains across the western United States. These volumes are difficult to compare because they represent the response to rainstorms and runoff with different characteristics. However, by measuring the erosion response as the erodibility efficiency of water to detach and transport sediment on hillslopes and in channels, the erosion response from different geologic terrains can be compared. Specifically, the erodibility efficiency is the percentage of the total available stream power expended to detach, remobilize, or transport a mass of sediment. Erodibility efficiencies were calculated for the (i) initial detachment, and for the (ii) remobilization and transport of sediment on the hillslopes and in the channels after wildfire in two different geological terrains.The initial detachment efficiencies for the main channel and tributary channel in the granitic terrain were 10 ± 9% and 5 ± 4% and were similar to those for the volcanic terrain, which were 5 ± 5% and 1 ± 1%. No initial detachment efficiency could be measured for the hillslopes in the granitic terrain because hillslope measurements were started after the first major rainstorm. The initial detachment efficiency in the volcanic terrain was 1.3 ± 0.41%. The average remobilization and transport efficiencies associated with flash floods in the channels also were similar in the granitic (0.18 ± 0.57%) and volcanic (0.11 ± 0.41%) terrains. On the hillslope the remobilization and transport efficiency was greater in the volcanic terrain (2.4%) than in the granitic terrain (0.65%). However, this may reflect the reduced sediment availability after the first major rainstorm (30-min maximum rainfall intensity  90 mm h^− 1) in the granitic terrain, while easily erodible fine colluvium remained on the hillslope after the first rainstorm (30-min maximum rainfall intensity = 7.2 mm h^− 1) in the volcanic terrain. The erosion response in channels and on hillslopes of the granitic and volcanic terrains was similar when compared using erodibility efficiencies.     

Geomorphic inferences from regolith thickness, chemical denudation and CRN erosion rates near the glacial limit, Boulder Creek catchment and vicinity, Colorado

Granitic regolith, developed in the Boulder Creek catchment and adjacent areas, records a history of deep weathering, some of which may predate Quaternary time. Field and well-log measurements of weathering, chemical denudation and rates of erosion derived from ¹⁰Be cosmogenic radionuclide (CRN) data help to quantify rates of landscape change in the post-orogenic Rocky Mountains. The density of oxidized, fractured bedrock ranges from 2.7 to about 2.2 g cm^− 3, saprolite and grus have densities between 2.0 and 1.8 g cm^− 3, and 30 soil samples averaged 1.6 ± 0.2 g cm^− 3. Highly weathered regolith in 540 wells averages 3.3 m thick, mean depth to bedrock in 1661 wells is 7 m, and the weathered thickness exceeds 10 m in relatively large local areas east of the late Pleistocene glacial limit. Thickness of regolith shows no simple relationship to rock type or structure, local slope, or distance from channels. Catchments in the vicinity of the Boulder Creek have an average CRN erosion rate of 2.2 ± 0.7 cm kyr^− 1 for the past 10,000 to 40,000 yr. Annual losses of cations and SiO₂ vary from about 2 to 5 g m^− 2 over a runoff range of 10 to nearly 160 cm.Using measured rates in simple box models shows that if a substantial fraction of void space is created by volume expansion in the weathering rock materials, 7 m of weathered rock materials could form in as little as 230 kyr. If density loss results mainly from chemical denudation and some volume expansion, however, the same weathering profile would take > 1340 kyr to form. Rates of erosion measured by CRN could be balanced by the rate of soil formation from saprolite if the annual solute loss from soil is 2.0 g m^− 2 and 70% of the density decrease from saprolite to grus and soil results from strain. Saprolite, however, forms from oxidized bedrock at a far slower rate and rates of saprolite formation cannot balance soil and grus losses to erosion. The zone of thick weathered regolith is likely an eroding relict landscape. The undulating surface marked by relatively low relief and tors is not literally a topographic surface of Eocene, Oligocene or Miocene age unless it was covered with deposits that were removed in Pliocene or Quaternary time.     

Landscape disequilibrium on 1000–10,000 year scales Marsyandi River, Nepal, central Himalaya

In an actively deforming orogen, maintenance of a topographic steady state requires that hillslope erosion, river incision, and rock uplift rates are balanced over timescales of 10⁵–10⁷ years. Over shorter times, <10⁵ years, hillslope erosion and bedrock river incision rates fluctuate with changes in climate. On 10⁴-year timescales, the Marsyandi River in the central Nepal Himalaya has oscillated between bedrock incision and valley alluviation in response to changes in monsoon intensity and sediment flux. Stratigraphy and ¹⁴C ages of fill terrace deposits reveal a major alluviation, coincident with a monsoonal maximum, ca. 50–35 ky BP. Cosmogenic ¹⁰Be and ²⁶Al exposure ages define an alluviation and reincision event ca. 9–6 ky BP, also at a time of strong South Asian monsoons. The terrace deposits that line the Lesser Himalayan channel are largely composed of debris flows which originate in the Greater Himalayan rocks up to 40 km away. The terrace sequences contain many cubic kilometers of sediment, but probably represent only 2–8% of the sediments which flushed through the Marsyandi during the accumulation period. At 10⁴-year timescales, maximum bedrock incision rates are 7 mm/year in the Greater Himalaya and 1.5 mm/year in the Lesser Himalayan Mahabarat Range. We propose a model in which river channel erosion is temporally out-of-phase with hillslope erosion. Increased monsoonal precipitation causes an increase in hillslope-derived sediment that overwhelms the transport capacity of the river. The resulting aggradation protects the bedrock channel from erosion, allowing the river gradient to steepen as rock uplift continues. When the alluvium is later removed and the bedrock channel re-exposed, bedrock incision rates probably accelerate beyond the long-term mean as the river gradient adjusts downward toward a more “equilibrium” profile. Efforts to document dynamic equilibrium in active orogens require quantification of rates over time intervals significantly exceeding the scale of these millennial fluctuations in rate.     

Spatial patterns of old, deep-seated landslides: A case-study in the northern Ethiopian highlands

During the last decade, slope failures were reported in a 500 km² study area in the Geba–Werei catchment, northern Ethiopia, a region where landslides were not considered an important hazard before. Field observations, however, revealed that many of the failures were actually reactivations of old deep-seated landslides after land use changes. Therefore, this study was conducted (1) to explore the importance of environmental factors controlling landslide occurrence and (2) to estimate future landslide susceptibility. A landslide inventory map of the study area derived from aerial photograph interpretation and field checks shows the location of 57 landslides and six zones with multiple landslides, mainly complex slides and debris flows. In total 14.8% of the area is affected by an old landslide. For the landslide susceptibility modelling, weights of evidence (WofE), was applied and five different models were produced. After comparison of the models and spatial validation using Receiver Operating Characteristic curves and Kappa values, a model combining data on elevation, hillslope gradient, aspect, geology and distance to faults was selected. This model confirmed our hypothesis that deep-seated landslides are located on hillslopes with a moderate slope gradient (i.e. 5°–13°). The depletion areas are expected on and along the border of plateaus where weathered basalts rich in smectite clays are found, and the landslide debris is expected to accumulate on the Amba Aradam sandstone and upper Antalo limestone. As future landslides are believed to occur on inherently unstable hillslopes similar to those where deep-seated landslides occurred, the classified landslide susceptibility map allows delineating zones where human interventions decreasing slope stability might cause slope failures. The results obtained demonstrate that the applied methodology could be used in similar areas where information on the location of landslides is essential for present-day hazard analysis.     

Rock weathering processes and landform development in the Sør Rondane Mountains, Antarctica

Field observations of weathering processes and the related landforms, combined with laboratory analyses of weathering products, permit a synthetic evaluation of Late Cenozoic weathering environments in the Sør Rondane Mountains, Antarctica, an arid upland characterized by low temperatures and strong winds. Rates and character of weathering depend mainly on moisture availability and the bedrock geology. Under the humid weathering regime that occurs only locally around the margin of the present sheet, frequent diurnal freeze-thaw cycles in summer cause relatively rapid rock fragmentation. Most of the mountains are situated in the arid weathering regime, under which rock breakdown is very slow unless the rock contains plenty of salts. Salt weathering becomes more intensive and extensive with exposure age, as a result of salt accumulation in rock, eventually producing soils as small as fine-silt size. Lack of clay mineralization even in weathered rocks having been exposed above the ice sheet prior to 4 Ma ago indicates that hydrolysis or carbonation of rock minerals has been insignificant during the past 4 Ma. The final products of weathering are due mainly to salt action and reflect the parent lithology. Resistant fine-grained granite forms strongly oxidized tors carved with tafoni, or fields of mushroom-like boulders overlying the fractured bedrock. Less resistant rocks, like biotite gneiss and amphibolite, produce stone pavements underlain by saline, silty soils up to 30–40 cm thick, the thickness of which corresponds to the maximum thaw depth.     

Contrasting soils and landscapes of the Piedmont and Coastal Plain, eastern United States

The Piedmont and Coastal Plain physiographic provinces comprise 80 percent of the Atlantic Coastal states from New Jersey to Georgia. The provinces are climatically similar. The soil moisture regime is udic. The soil temperature regime is typically thermic from Virginia through Georgia, although it is mesic at altitudes above 400 m in Georgia and above 320 m in Virginia. The soil temperature regime is mesic for the Piedmont and Coastal Plain from Maryland through New Jersey. The tightly folded, structurally complex crystalline rocks of the Piedmont and the gently dipping “layer-cake” clastic sedimentary rocks and sediments of the Coastal Plain respond differently to weathering, pedogenesis, and erosion. The different responses result in two physiographically contrasting terrains; each has distinctive near-surface hydrology, regolith, drainage morphology, and morphometry.The Piedmont is predominantly an erosional terrain. Interfluves are as narrow as 0.5 to 2 km, and are convex upward. Valleys are as narrow as 0.1 to 0.5 km and generally V-shaped in cross section. Alluvial terraces are rare and discontinuous. Soils in the Piedmont are typically less than 1 m thick, have less sand and more clay than Coastal Plain soils, and generally have not developed sandy epipedons. Infiltration rates for Piedmont soils are low at 6–15 cm/h. The soil/saprolite, soil/rock, and saprolite/rock boundaries are distinct (can be placed within 10 cm) and are characterized by ponding and/or lateral movement of water. Water movement through soil into saprolite, and from saprolite into rock, is along joints, foliation, bedding planes and faults. Soils and isotopic data indicate residence times consistent with a Pleistocene age for most Piedmont soils.The Coastal Plain is both an erosional and a constructional terrain. Interfluves commonly are broader than 2 km and are flat. Valleys are commonly as wide as 1 km to greater than 10 km, and contain numerous alluvial and estuarine terrace sequences that can be correlated along valleys for tens of kilometers. Coastal Plain soils are typically as thick as 2 to 8 m, have high sand content throughout, and have sandy epipedons. These epipedons consist of both A and E horizons and are 1 to 4 m thick. In Coastal Plain soils, the boundaries are transitional between the solum and the underlying parent material and between weathered and unweathered parent material. Infiltration rates for Coastal Plain soils are typically higher at 13–28 cm/h, than are those for Piedmont soils. Indeed, for unconsolidated quartz sand, rates may exceed 50 cm/h. Water moves directly from the soil into the parent material through intergranularpores with only minor channelization along macropores, joints, and fractures. The comparatively high infiltration capacity results in relatively low surface runoff, and correspondingly less erosion than on the Piedmont uplands.Due to differences in Piedmont and Coastal Plain erosion rates, topographic inversion is common along the Fall Zone; surfaces on Cenozoic sedimentary deposits of the Coastal Plain are higher than erosional surfaces on regolith weathered from late Precambrian to early Paleozoic crystalline rocks of the Piedmont. Isotopic, paleontologic, and soil data indicate that Coastal Plain surficial deposits are post-middle Miocene to Holocene in age, but most are from 5 to 2 Ma. Thus, the relatively uneroded surfaces comprise a Pliocene landscape. In the eastern third of the Coastal Plain, deposits that are less than 3.5 Ma include alluvial terraces, marine terraces and barrier/back-barrier complexes as morphostratigraphic units that cover thousands of square kilometers. Isotopic and soil data indicate that eastern Piedmont soils range from late Pliocene to Pleistocene in age, but are predominantly less than 2 Ma old. Thus, the eroded uplands of the Piedmont “peneplain” comprise a Pleistocene landscape.     

Hillslope gullying in the Solway Firth — Morecambe Bay region, Great Britain: Responses to human impact and/or climatic deterioration?

In the Solway Firth — Morecambe Bay region of Great Britain there is evidence for heightened hillslope instability during the late Holocene (after 3000 cal. BP). Little or no hillslope geomorphic activity has been identified occurring during the early Holocene, but there is abundant evidence for late Holocene hillslope erosion (gullying) and associated alluvial fan and valley floor deposition. Interpretation of the regional radiocarbon chronology available from organic matter buried beneath alluvial fan units suggests much of this geomorphic activity can be attributed to four phases of more extensive gullying identified after 2500–2200, 1300–1000, 1000–800 and 500 cal. BP. Both climate and human impact models can be evoked to explain the crossing of geomorphic thresholds: and palaeoecological data on climatic change (bog surface wetness) and human impact (pollen), together with archaeological and documentary evidence of landscape history, provide a context for addressing the causes of late Holocene geomorphic instability. High magnitude storm events are the primary agent responsible for gully incision, but neither such events nor cooler/wetter climatic episodes appear to have produced gully systems in the region before 3000 cal. BP. Increased gullying after 2500–2200 cal. BP coincides with population expansion during Iron Age and Romano-British times. The widespread and extensive gullying after 1300–1000 cal. BP and after 1000–800 cal. BP coincides with periods of population expansion and a growing rural economy identified during Norse times, 9–10th centuries AD, and during the Medieval Period, 12–13th centuries AD. These periods were separated by a downturn associated with the ‘harrying of the north’ AD 1069 to 1070. The gullying episode after 500 cal. BP also coincides with increased anthropogenic pressure on the uplands, with population growth and agricultural expansion after AD 1500 following 150 years of malaise caused by livestock and human (the Black Death) plagues, poor harvests and conflicts on the Scottish/English border. The increased susceptibility to erosion of gullies is a response to increased anthropogenic pressure on upland hillslopes during the late Holocene, and the role of this pressure appears crucial in priming hillslopes before subsequent major storm events. In particular, the cycles of expansion and contraction in both population and agriculture appear to have affected the susceptibility of the upland landscape to erosion, and the hillslope gullying record in the region, therefore, contributes to understanding of the timing and spatial pattern of human exploitation of the upland landscape.     

Analysis of soil surface component patterns affecting runoff generation. An example of methods applied to Mediterranean hillslopes in Alicante (Spain)

Spatial patterns of soil surface components (vegetation, rock fragments, crusts, bedrock outcrops, etc.) are a key factor determining hydrological functioning of hillslopes. A methodological approach to analyse the patterns of soil surface components at a detailed scale is proposed in this paper. The methods proposed are applied to two contrasting semi-arid Mediterranean hillslopes, and the influence of soil surface component patterns on the runoff response of the slopes was analysed. A soil surface components map was derived from a high resolution photo-mosaic obtained in the field by means of a digital camera. Rainfall simulation experimental data were used to characterise the hydrological behaviour of areas with a specific pattern of soil surface components by means of the parameters of the Horton equation. Plot runoff data were extrapolated at the hillslope scale based on the soil surface component maps and their hydrological characterisation. The results show that in both slopes runoff generation is concentrated up- and downslope, with a water accepting area in the centre of both slopes disrupting the hydrological connectivity at the slope scale. This reinfiltration patch at the centre of the slope is related to the type of soil surface component and its spatial pattern. Herbaceous vegetation and ‘on top rock fragments’ increase the infiltration capacity of soils at the centre of the slope. In contrast, embedded rock fragments, rock outcrops, as well as crusted surfaces located in the upper and lower slopes favour runoff generation in these areas. In addition, a general pattern of water contribution areas downslope is apparent on both slopes. The south-facing slope shows a higher hydrological connectivity and more runoff. 55% of the surface of the south-facing slope produces runoff at the end of a 1 hour rainfall event and 17.3% of the surface is covered by a runoff depth between 0.5 and 1 mm. While on the north-facing slope only 38% of the surface produces runoff under the same conditions. Longitudinal connectivity of runoff is higher at the south-facing slope where more runoff-generating surfaces appear and where the vegetation pattern favours the connectivity of bare areas.     

The role of soil processes in determining mechanisms of slope failure and hillslope development in a humid-tropical forest eastern Puerto Rico

Translational failures, with associated downslope earthflow components and shallow slides, appear to be the primary mechanism of hillslope denudation in the humid tropical forests of the mountains of eastern Puerto Rico. In-situ weathering of quartz diorite and marine-deposited volcaniclastics produces residual soil (saprolite; up to 21 m deep) / weathered rock profiles. Discontinuous zones of contrasting density and permeability particularly in quartz-diorite slopes at 0.5 m, and between 3 and 7 m, create both pathways and impedances for water that can result in excess pore pressures and, ultimately, aid in determining the location of failure planes and magnitudes of slope failures. In combination with relict fractures which create planes of weakness within the saprolite, and the potential significance of tensile stresses in the upper zone of saprolite (hypothesized to be caused by subsurface soil creep), shear failure can then occur during or after periods of heavy rainfall.Results of in-situ shear-strength testing show negative y-intercepts on the derived Mohr-Coulomb failure envelopes (approximately 50% of all tests) that are interpreted as apparent tensile stresses. Observation of tension cracks 1–2 m deep support the test data. Subsurface soil creep can cause extension of the soil and the development of tensile stresses along upper-slope segments. Shear-strength data support this hypothesis for both geologic types. Apparent values of maximum and mean tensile stress are greatest along upper slopes (16.5 and 6.29 kPa). Previously documented maximum rates of downslope movement coincided with local minima of shear strength, and the shear-strength minimum for all tests was located near 0.5 m below land surface, the shallow zone of contrasting permeabilities. These results indicate that subsurface soil creep, a slow semi-continuous process, may exert a profound influence on rapid, shallow slope failures in saprolitic soils.Data indicate that cove slopes in quartz diorite tend to be the most unstable when saturation levels reach 75%. Deep failures (7 m deep) appear the most critical but not the most frequent because pore pressure build-up will occur more rapidly in the upper perched zone of translocated clays before reaching the lower zone between 3 and 7 m. Frequent shallow failures could reduce the probability of deeper failures by removing overburden and reducing shear stress at depth. Deep failures are more likely to result from storm events of great duration and intensity.Sixty-six ‘naturally occurring’ and more than 100 ‘road-related’ landslides were mapped. Forest elevations exceed 1000 m, but the majority of these failures were found between 600 and 800 m in elevation. This appears to be the area where there is sufficient concentration of subsurface water to result in excess pore pressures. The high percentage of slope failures in the 600–800-m range, relative to the percentage at higher elevations, suggests that differences in soil-water processes are responsible for the form of these mountain slopes. Steep linear segments are maintained at higher elevations. Slope angles are reduced in the 600–800-m range by frequent shallow slides, creating a largely concave surface. In combination, slope segments above 800 m, and those between 600 and 800 m, produce the characteristic form of the mountains of eastern Puerto Rico.     

Dating of small gully formation and establishing erosion rates in old gullies under forest by means of anatomical changes in exposed tree roots (Southern Poland)

Small gullies occur in forested gully systems on the undulating loess plateau in southern Poland. The old gully hillslopes are mainly covered with 200-year old beech trees in contrast with the surface of the summit plateau, which is cultivated agricultural land. Beech roots are exposed in the gullies through erosion. Wood vessels in the root tree rings divide into early wood and late wood and, after the roots are exposed, start to make fewer vessels. These anatomical changes in root tree rings allow us to date erosion episodes.Small gullies form in a different manner on the valley floor and on hillslopes. In valley bottoms, erosion features are often formed at some distance from one another, and in time small gullies combine to form a single, longer one. Depending on local conditions, such as the hillslope profile, hillslopes may exhibit headward erosion or may be eroded downwards. Hilllope gullies may be transformed into side valleys as a result of gradual widening and deepening.Dating the exposure of roots indicates that small gullies had already formed in the valley system by 1949. Intensive gully erosion was recorded between 1984 and 2002, during intense precipitation in 1984 and, of particular note, during the extraordinary flood of 1997 which affected all of Central Europe. The mean rate of small gully erosion in the old gully system studied is 0.63 m/year. On hillslopes the mean gully erosion rate is 0.21–0.52 m/year, and on the valley bottoms 0.18–1.98 m/year. High bottom erosion rates resulted from the emergence of long gullies during the erosion episodes in 1984 and 1997. Sheet flow in valley floors intensifies at times of heavy rainfall which causes long gullies to form.Taking into account the fact that conditions favoured erosion, the rate at which the old gullies under forest were transformed should be considered slow. New side gullies form slowly within the valley and it appears that if erosion progressed at the rate observed, new side valleys would take a few hundred years to form.     

Influence of rock strength on the valley morphometry of Big Creek, central Idaho, USA

Analysis of valley morphometry and bedrock strength along Big Creek, central Idaho, shows that valley floor width is strongly controlled by bedrock. We performed statistical analysis of Schmidt hammer rock strength as a function of lithology and aspect and of valley morphometry as a function of rock strength. Rock strength is significantly greater on the south side of the valley and in Eocene granodiorites. Rock strength is weakest in Eocene volcanic tuffs. Valley floor width depends negatively on weakest valley-side rock strength, and hillslope gradient on the north side of the valley depends positively on rock strength. Stream gradient does not depend on rock strength. Valley floor width appears to be controlled by bedrock strength on the weaker side of the valley, which was generally the north (south-facing) side. We speculate that a higher degree of weathering via freeze–thaw cycles contributes to lower strength on the north side. The positive dependence of hillslope gradient on rock strength on the north side provides evidence that differential weathering across lithologies determines the gradient that can be maintained as lateral migration of the stream erodes valley walls. These results suggest that in situ rock strength exerts strong influences on some measures of valley morphometry by modulating hillslope mass wasting processes and limiting lateral erosion.     

A Karst Connection model for Grand Canyon, Arizona, USA

Along the Asturian coast of northern Spain an uplifted wave-cut platform extends for  100 km east–west. The steep cliff which bounds the gently seaward-dipping platform to the north increases in height from 30 m in the west to 100 m in the east and reflects the overall eastward increase in platform elevation. The southern edge of the 2–4 km-wide platform runs along the foothills of the Cantabrian Mountains, as constrained by a high-resolution digital elevation model. The marine platform, which was carved into deformed Paleozoic bedrock with abundant quartzite beds, is largely covered by weathered marine and continental sediments. Quartzite samples from flat bedrock outcrops which are currently not covered by sediment or soil yield cosmogenic nuclide concentrations (²¹Ne, ¹⁰Be and ²⁶Al) that demonstrate a long and complex exposure history, including periods of burial with partial or complete shielding from cosmic rays. The combination of multiple cosmogenic nuclides yields a minimum age of 1–2 Ma for the platform. Taking into account (i) the horizontal and vertical extent of the platform, (ii) the high resistance to erosion of the quartzitic bedrock, and (iii) published data on the magnitude of past sea level fluctuations, we suggest that the wave-cut platform formed in the Pliocene. Subvertical faults cutting the platform at high angles to the coastline offset the southern edge of the platform by 20 to 40 m and reactivate the pre-existing anisotropy in the Paleozoic bedrock. Uplift and crustal deformation of the coastal region have occurred after platform formation in the Pliocene and may still be active. The slow deformation of the northern edge of the Iberian plate including the Cantabrian Mountains may result from the ongoing slow convergence at an incipient subduction zone extending along the coast of northern Spain.