Southern Appalachian hillslope erosion rates measured by soil and detrital radiocarbon in hollows

Understanding the dynamics of sediment generation and transport on hillslopes provides important constraints on the rate of sediment output from orogenic systems. Hillslope sediment fluxes are recorded by organic material found in the deposits infilling unchanneled convergent topographic features called hollows. This study describes the first hollow infilling rates measured in the southern Appalachian Mountains. Infilling rates (and bedrock erosion rates) were calculated from the vertical distribution of radiocarbon ages at two sites in the Coweeta drainage basin, western North Carolina. At each site we dated paired charcoal and silt soil organic matter samples from five different horizons. Paired radiocarbon samples were used to bracket the age of the soil material in order to capture the range of complex soil forming processes and deposition within the hollows. These dates constrain hillslope erosion rates of between 0.051 and 0.111 mm yr^− 1. These rates are up to 4 times higher than spatially-averaged rates for the Southern Appalachian Mountains making creep processes one of the most efficient erosional mechanisms in this mountain range. Our hillslope erosion rates are consistent with those of forested mountain ranges in the western United States, suggesting that the mechanisms (dominantly tree throw) driving creep erosion in both the western United States and the Southern Appalachian Mountains are equally effective.     

Landscape evolution and glaciation of the Rwenzori Mountains, Uganda: Insights from numerical modeling

The Rwenzori Mountains are a high alpine mountain chain, about 40 × 80 km in size, just north of the equator in the western branch of the East African Rift System in Africa. The central part of the mountain chain is located in Uganda, and the highest peak, the Margherita Peak with 5119 m, lies on the border to the Democratic Republic of Congo. Topography is very pronounced, with steeply incised valleys and clear glacial landforms in the upper part of the mountain chain. The Rwenzori Mountains are an unusually high mountain chain located in the extensional setting of the East African Rift System, and the large elevation poses a challenging problem for geodynamists to explain.We have used the landscape evolution model ULTIMA THULE, which combines hillslope diffusion, fluvial erosion, and glacial abrasion and is driven by a climate driver, simulating the variations in temperature, precipitation, and relief over several glacial cycles. With a simulation time of 800 ka, we test the hypothesis of climate-tectonic interactions on the uplift of the Rwenzori Mountains.Our results show that a moderate cooling of around 6° causes widespread glaciation of the high mountain regions as observed during the peak glacial phases, and that morphological processes degrading the landscape allow for a tectonic uplift rate of around 0.5 mm a^− 1.     

Relating accretion and erosion at an exposed tidal wetland to the bottom shear stress of combined current-wave action

Sediment dynamics have an important influence on the morphological evolution of tidal wetlands, which consist of mudflats and salt marshes. To understand the nature of sediment behavior under combined current-wave action at an exposed tidal wetland, we measured the waves, currents, water depths, bed-level changes, and sediment properties at a mudflat-salt marsh transition on the Yangtze Delta, China, during five consecutive tides under onshore winds of ~ 8 m/s, and calculated the bed shear stresses due to currents (τ_c), waves (τ_w), combined current-wave action (τ_cw), and the critical shear stress for erosion of the bottom sediment (τ_ce). The bed shear stresses under combined current-wave action (τ_cw) were approximately five times higher on the mudflat (up to 1.11 N/m²; average 0.27 N/m²) than on the salt marsh (up to 0.14 N/m²; average, 0.06 N/m²). On the mudflat, τ_cw was larger than the critical erosion shear stress (τ_ce = 0.103 N/m²) for 70% of the period of submergence, whereas τ_cw was always lower than τ_ce at the salt marsh site (τ_ce = 0.116 N/m²). This result indicates that the sediment dynamics on the mudflat were dominated by erosion, whereas at the salt marsh they were governed by deposition, which is in agreement with the observed bed-level change during the study period (− 3.3 mm/tide on the mudflat and 3.0 mm/tide on the salt marsh). A comparison of τ_cw values calculated using the [van Rijn, 1993] and [Soulsby, 1995] models for bed shear stresses under combined current-wave action indicates that both models are applicable to the present case and effectively predict the bottom shear stress under combined current-wave action. Overall, we conclude that τ_cw in combination with τ_ce is useful in assessing the hydrodynamic mechanisms that underlie the morphological evolution of exposed tidal wetlands.     

The influence of upper-crust lithology on topographic development in the central Coast Ranges of California

A fundamental geological tenet is that as landscapes evolve over graded to geologic time, geologic structures control patterns of topographic distribution in mountainous areas such that terrain underlain by competent rock will be higher than terrain underlain by incompetent rock. This paper shows that in active orogens where markedly weak and markedly strong rocks are juxtaposed along contacts that parallel regional structures, relatively high topography can form where strain is localized in the weak rock. Such a relationship is illustrated by the topography of the central Coast Ranges between the Pacific coastline and the San Andreas fault zone (SAFZ), and along the length of the Gabilan Mesa (the “Gabilan Mesa segment” of the central Coast Ranges). Within the Gabilan Mesa segment, the granitic upper crust of the Salinian terrane is in contact with the accretionary-prism mélange upper crust of the Nacimiento terrane along the inactive Nacimiento fault zone. A prominent topographic lineament is present along most of this lithologic boundary, approximately 50 to 65 km southwest of the SAFZ, with the higher topography formed in the mélange on the southwest side of the Nacimiento fault.This paper investigates factors influencing the pattern of topographic development in the Gabilan Mesa segment of the central Coast Ranges by correlating shortening magnitude with the upper-crust compositions of the Salinian and Nacimiento terranes. The fluvial geomorphology of two valleys in the Gabilan Mesa, which is within the Salinian terrane, and alluvial geochronology based on optically-stimulated luminescence (OSL) age estimates, reveal that the magnitude of shortening accommodated by down-to-the-southwest tilting of the mesa since 400 ka is less than 1 to 2 m. Our results, combined with those of previous studies, indicate that at least 63% to 78% of late-Cenozoic, northeast-southwest directed, upper-crustal shortening across the Gabilan Mesa segment has been accommodated within the Nacimiento terrane. This is significant because perpendicular to orogenic strike the Nacimiento terrane constitutes less than ¼ of the distance between the coast and the SAFZ, and the other ¾ (or greater) of the distance between the coast and the SAFZ is underlain by the granitic upper crust of the Salinian terrane. We propose that strain and mountain building are localized within the Nacimiento terrane because it consists predominantly of the relatively weak Franciscan Complex mélange, and because the upper crust of the Salinian terrane is composed of relatively strong granitic rocks. Our hypothesis is supported by the distribution of post-seismic surface uplift associated with the 2003, 6.5 M_W San Simeon earthquake, which mimics the topography of the southwestern part of the Gabilan Mesa segment of the central Coast Ranges.     

Neotectonic evolution of the northern Laborec drainage basin (northeastern part of Slovakia)

This study investigates the influence of neotectonic activity on river and basin patterns in a mountainous area located in the northeastern part of the Carpathian Belt (the Laborecká vrchovina and Bukovské vrchy Mts. in eastern Slovakia). This area evolved within the accretionary wedge of the Carpathians during the Neogene, and it was alsowas affected by Middle to Late Miocene thrusting of the External Carpathians. Morphometric analysis, longitudinal and transverse river valley profiles, analysis of basin and valley symmetries, and investigation of alluvial terraces were carried out on the northern Laborec River and its tributaries. This was done to detect a possible relationship between their river courses and any ongoing neotectonic activity, which is otherwise difficult to detect by methods of structural geology because of the poorly exposed area.The general topography of the basin is characterized by a stepwise inclination to the SW as a result of differential uplift and subsidence. The reorganization of the river network in the Laborec drainage basin was influenced by tectonic activity along the NE-SW up to N-S fault structures during the neotectonic phase (Pliocene-Quaternary). The movement along these fault structures is predominantly normal to transtensive. The obtained data assumes that the region is under approximately NE-SW oriented S_H compression and NW-SE trending S_h tension. The Laborec drainage basin is characterized by a very high degree of asymmetry that sharply increases from the upper to the lower courses of the river. The right-bank tributaries of the Laborec River are < 12 km in length; however, the left-bank tributaries such as Vydraňka, Ol'šava, Výrava, Udava, and Cirocha Streams are up to 50 km long with a high potential of headward erosion and capturing. The valley asymmetry is also very variable in the upper and lower portions of the basin. Based on these presented results, the ancient river thalweg was located along the axis of the Hostovice-Habura depression, and it was captured by the Ol'šava, Výrava, and Udava Streams. The asymmetric pattern of the drainage basin is the result of active tectonics, the continual subsidence of the Transcarpathian Basin, and by the uplift of the Laborecká vrchovina and Bukovské vrchy Mts. These events caused rejuvenation of the headward erosion of streams in the southern part. Favorable lithology was also essential in the process of river capture.     

Uplift, exhumation and precipitation: tectonic and climatic control of Late Cenozoic landscape evolution in the northern Sierras Pampeanas, Argentina

Deciphering the evolution of mountain belts requires information on the temporal history of both topographic growth and erosion. The exhumation rate of a mountain range undergoing shortening is related to the erodability of the uplifting range as well as the efficiency of erosion, which partly depends on the available precipitation. Young, rapidly deposited sediments have low thermal conductivity and are readily eroded, in contrast to underlying resistant basement rocks that have a higher thermal conductivity. Apatite fission‐track thermochronology can quantify cooling; thermal models constrain the relationship between this cooling and exhumation. By utilizing geological relations for a datum, we can examine the evolution of rock uplift, surface uplift and exhumation. In the northern Sierras Pampeanas of Argentina, a young sedimentary basin that overlay resistant crystalline basement prior to rapid exhumation provides an ideal setting to examine the effect of contrasting thermal and erosional regimes. There, tectonically active reverse‐fault‐bounded blocks partly preserve a basement peneplain at elevations in excess of 4500 m. Prior to exhumation, the two study areas were covered by 1000 and 1600 m of recently deposited sediments; this sequence begins with shallow marine deposits immediately overlying the regional erosion surface. Apatite fission‐track data were obtained from vertical transects in the Calchaquíes and Aconquija ranges. At Cumbres Calchaquíes, erosion leading to the development of the peneplain commenced in the Cretaceous, probably as a result of rift‐shoulder uplift. In contrast, Sierra Aconquija cooled rapidly between 5.5 and 4.5 Myr. At the onset of this rapid exhumation, the sediment was quickly removed, causing fast cooling, but relatively slow rates of surface uplift. Syntectonic conglomerates were produced when faulting exposed resistant bedrock; this change in rock erodability led to enhanced surface uplift rates, but decreased exhumation rates. The creation of an orographic barrier after the range had attained sufficient elevation further decreased exhumation rates and increased surface uplift rates. Differences in the magnitude of exhumation at the two transects are related to both differences in the thickness of the sedimentary basin prior to exhumation and differences in the effective precipitation due to an orographic barrier in the foreland and hence differences in the magnitude of headward erosion.     

A frost “buzzsaw” mechanism for erosion of the eastern Southern Alps, New Zealand

In the Southern Alps, New Zealand, large gradients in precipitation (< 1 to 12 m year^− 1) and rock uplift (< 1 to 10 mm year^− 1) produce distinct post-glacial geomorphic domains in which landslide-driven sediment production dominates in the wet, rapid-uplift western region, and rockfall controls erosion in the drier, low-uplift eastern region. Because the western region accounts for < 25% of the active orogen, the dynamics of erosion in the extensive eastern region are of equal importance in estimating the relative balance of uplift and erosion across the Southern Alps. Here, we assess the efficacy of frost cracking as the primary rockfall mechanism in the eastern Southern Alps using air photo and topographic analysis of scree slopes, cosmogenic radionuclide dating of headwalls, paleo-climate data, and a numerical model of headwall temperature. Currently, active scree slopes occur at a relatively uniform mean elevation ( 1450 m) and their distribution is independent of hillslope aspect and rock type, consistent with the notion that frost cracking (which is maximized between − 3 and − 8 °C) may control rockfall erosion. Headwall erosion rates of 0.3 to 0.9 mm year^− 1, measured using in-situ ¹⁰Be and ²⁶Al in the Cragieburn Range, confirm that rockfall erosion is active in the late Holocene at rates that roughly balance rock uplift. Models of the predicted depth of frost activity are consistent with the scale of fractures and scree blocks in our field sites. Also, vegetated, paleo-scree slopes are ubiquitous at elevations lower than active scree slopes, consistent with the notion that lower temperatures during the last glacial advance induced pervasive rockfall erosion due to frost cracking. Our modeling suggests temporally-averaged peak frost cracking intensity occurs at 2300 m a.s.l., the approximate elevation of the highest peaks in the central Southern Alps, suggesting that the height of these peaks may be limited by a “frost buzzsaw.”     

Small valley glaciers and the effectiveness of the glacial buzzsaw in the northern Basin and Range, USA

The glacial buzzsaw hypothesis suggests that efficient erosion limits topographic elevations in extensively glaciated orogens. Studies to date have largely focussed on regions where large glaciers (tens of kilometres long) have been active. In light of recent studies emphasising the importance of lateral glacial erosion in lowering peaks and ridgelines, we examine the effectiveness of small glaciers in limiting topography under both relatively slow and rapid rock uplift conditions. Four ranges in the northern Basin and Range, Idaho, Montana, and Wyoming, USA, were chosen for this analysis. Estimates of maximum Pleistocene slip rates along normal faults bounding the Beaverhead–Bitterroot Mountains (~ 0.14 mm y^− 1), Lemhi Range (~ 0.3 mm y^− 1) and Lost River Range (~ 0.3 mm y^− 1) are an order of magnitude lower than those on the Teton Fault (~ 2 mm y^− 1). We compare the distribution of glacial erosion (estimated from cirque floor elevations and last glacial maximum (LGM) equilibrium line altitude (ELA) reconstructions) and fault slip rate with three metrics of topography in each range: the along-strike maximum elevation swath profile, hypsometry, and slope-elevation profiles. In the slowly uplifting Beaverhead–Bitterroot Mountains, and Lemhi and Lost River Ranges, trends in maximum elevation parallel ELAs, independent of variations in fault slip rate. Maximum elevations are offset ~ 500 m from LGM ELAs in the Lost River Range, Lemhi Range, and northern Beaverhead–Bitterroot Mountains, and by ~ 350 m in the southern Beaverhead–Bitterroot Mountains, where glacial extents were less. The offset between maximum topography and mean Quaternary ELAs, inferred from cirque floor elevations, is ~ 350 m in the Lost River and Lemhi Ranges, and 200–250 m in the Beaverhead–Bitterroot Mountains. Additionally, slope-elevation profiles are flattened and hypsometry profiles show a peak in surface areas close to the ELA in the Lemhi Range and Beaverhead–Bitterroot Mountains, suggesting that small glaciers efficiently limit topography. The situation in the Lost River Range is less clear as a glacial signature is not apparent in either slope-elevation profiles or the hypsometry. In the rapidly uplifting Teton Range, the distribution of ELAs appears superficially to correspond to maximum topography, hypsometry, and slope-elevations profiles, with regression lines on maximum elevations offset by ~ 700 and ~ 350 m from the LGM and mean Quaternary ELA respectively. However, Grand Teton and Mt. Moran represent high-elevation “Teflon Peaks” that appear impervious to glacial erosion, formed in the hard crystalline bedrock at the core of the range. Glacier size and drainage density, rock uplift rate, and bedrock lithology are all important considerations when assessing the ability of glaciers to limit mountain range topography. In the northern Basin and Range, it is only under exceptional circumstances in the Teton Range that small glaciers appear to be incapable of imposing a fully efficient glacial buzzsaw, emphasising that high peaks represent an important caveat to the glacial buzzsaw hypothesis.     

A comparison of permafrost prediction models along a section of Trail Ridge Road, Rocky Mountain National Park, Colorado, USA

The distribution of mountain permafrost along Trail Ridge Road (TRR) in Rocky Mountain National Park, Colorado, was modeled using ‘frost numbers’ and a ‘temperature of permafrost model’ (TTOP) in order to assess the accuracy of prediction models. The TTOP model is based on regional observations of air temperature and heat transfer functions involving vegetation, soil, and snow; whereas the frost number model is based on site-specific ratios of ground temperature measurements of frozen and thawed degree-days. Thirty HOBO© temperature data loggers were installed near the surface as well as at depth (30 to 85 cm). From mid-July 2008 to 2010, the mean annual soil temperature (MAST) for all surface sites was − 1.5 °C. Frost numbers averaged 0.56; TTOP averaged − 1.8 °C. The MAST was colder on western-facing slopes at high elevations. Surface and deeper probes had similar MASTs; however, deeper probes had less daily and seasonal variation. Another model developed at the regional scale based on proxy indicators of permafrost (rock glaciers and land cover) classified 5.1 km² of permafrost within the study area, whereas co-kriging interpolations of frost numbers and TTOP data indicated 2.0 km² and 4.6 km² of permafrost, respectively. Only 0.8 km² were common among all three models. Three boreholes drilled within 2 m of TRR indicate that permafrost does not exist at these locations despite each borehole being classified as containing permafrost by at least one model. Addressing model uncertainty is important because nutrients stored within frozen or frost-affected soils can be released and impact alpine water bodies. The uncertainty also exposes two fundamental problems: empirical models designed for high latitudes are not necessarily applicable to mountain permafrost, and the presence of mountain permafrost in the alpine tundra of the Colorado Front Range has not been validated.     

How does alluvial sedimentation at range fronts modify the erosional dynamics of mountain catchments?

At the geological time scale, the way in which the erosion of drainage catchments responds to tectonic uplift and climate changes depends on boundary conditions. In particular, sediment accumulation and erosion occurring at the edge of mountain ranges should influence the base level of mountain catchments, as well as sediment and water discharges. In this paper, we use a landform evolution model (LEM) to investigate how the presence of alluvial sedimentation at range fronts affects catchment responses to climatic or tectonic changes. This approach is applied to a 25 km × 50 km domain, in which the central part is uplifted progressively to simulate the growth of a small mountain range. The LEM includes different slope and river processes that can compete with each other. This competition leads to ‘transport‐limited’, ‘detachment‐limited’ or ‘mixed’ transport conditions in mountains at dynamic equilibrium. In addition, two end‐member algorithms (the channellized‐flow and the sheet‐flow regimes) have been included for the alluvial fan‐flow regime. The three transport conditions and the two flow algorithms represent six different models for which the responses to increase of rock uplift rate and/or cyclic variation of the precipitation rate are investigated. Our results indicate that addition of an alluvial apron increases the long‐term mountain denudation. In response to uplift, mountain rivers adapt their profile in two successive stages; first by propagation of an erosion wave and then by slowly increasing their channel gradients. During the second stage, the erosion rate is almost uniform across the catchment area at any one time, which suggests that dynamic equilibrium has been reached, although the balance between erosion and rock uplift rates has not yet been achieved. This second stage is initiated by the uplift of the mountain river outlets because of sedimentation aggradation at the mountain front. The response time depends on the type of water flow imposed on the alluvial fans domains (× by 1.5 for channelized flow regime and by 10 for the sheet flow one). Cyclic variations of precipitation rate generate cyclic incisions in the alluvial apron. These incision pulses create knick‐points in the river profile in the case of ‘detachment‐limited’ and ‘mixed’ river conditions, which could be mistaken for tectonically induced knick‐points. ‘Transport‐limited’ conditions do not create such knick‐points, but nevertheless trigger erosion in catchments. The feedbacks linked to sedimentation and erosion at range front can therefore control catchment incision or aggradation. In addition, random river captures in the range front trigger auto‐cyclic erosion pulses in the catchment, capable of generating incision–aggradation cycles.     

Quaternary landscape evolution and erosion rates for an intramontane Neogene basin (Guadix–Baza basin, SE Spain)

The landscape evolution in Neogene intramontane basins is a result of the interaction of climatic, lithologic, and tectonic factors. When sedimentation ceases and a basin enters an erosional stage, estimating erosion rates across the entire basin can offer a good view of landscape evolution. In this work, the erosion rates in the Guadix–Baza basin have been calculated based on a volumetric estimate of sediment loss by river erosion since the Late Pleistocene. To do so, the distribution of a glacis surface at ca. 43 kyr, characterised by a calcrete layer that caps the basin infilling, has been reconstructed. To support this age, new radiometric data of the glacis are presented. The volume of sediment loss by water erosion has been calculated for the entire basin by comparing the reconstructed geomorphic surface and the present-day topography. The resulting erosion rates vary between 4.28 and 6.57 m³ ha^− 1 yr^− 1, and are the consequence of the interaction of climatic, lithologic, topographic, and tectonic factors. Individual erosion rates for the Guadix and Baza sub-basins (11.80 m³ ha^− 1 yr^− 1 and 1.77 m³ ha^− 1 yr^− 1 respectively) suggest different stages of drainage pattern evolution in the two sub-basins. We attribute the lower values obtained in the Baza sub-basin to the down-throw of this sub-basin caused by very recent activity along the Baza fault.     

Comparison of conceptual landscape metrics to define hillslope-scale sediment delivery ratio

The aim of this study was to evaluate four metrics to define the spatially variable (regionalised) hillslope sediment delivery ratio (HSDR). A catchment model that accounted for gully and streambank erosion and floodplain deposition was used to isolate the effects of hillslope gross erosion and hillslope delivery from other landscape processes. The analysis was carried out at the subcatchment (~ 40 km²) and the cell scale (400 m²) in the Avon-Richardson catchment (3300 km²), south-east Australia. The four landscape metrics selected for the study were based on sediment travel time, sediment transport capacity, flux connectivity, and residence time. Model configurations with spatially-constant or regionalised HSDR were calibrated against sediment yield measured at five gauging stations. The impact of using regionalised HSDR was evaluated in terms of improved model performance against measured sediment yields in a nested monitoring network, the complexity and data requirements of the metric, and the resulting spatial relationship between hillslope erosion and landscape factors in the catchment and along hillslope transects. The introduction of a regionalised HSDR generally improved model predictions of specific sediment yields at the subcatchment scale, increasing model efficiency from 0.48 to > 0.6 in the best cases. However, the introduction of regionalised HSDR metrics at the cell scale did not improve model performance. The flux connectivity was the most promising metric because it showed the largest improvement in predicting specific sediment yields, was easy to implement, was scale-independent and its formulation was consistent with sedimentological connectivity concepts. These properties make the flux connectivity metric preferable for applications to catchments where climatic conditions can be considered homogeneous, i.e. in small-medium sized basins (up to approximately 3000 km² for Australian conditions, with the Avon-Richardson catchment being at the upper boundary). The residence time metric improved model assessment of sediment yields and enabled accounting for climatic variability on sediment delivery, but at the cost of greater complexity and data requirements; this metric might be more suitable for application in catchments with important climatic gradients, i.e. large basins and at the regional scale. The application of a regionalised HSDR metric did not increase data or computational requirements substantially, and is recommended to improve assessment of hillslope erosion in empirical, semi-lumped erosion modelling applications. However, more research is needed to assess the quality of spatial patterns of erosion depicted by the different landscape metrics.     

Woodland expansion’s influence on belowground carbon and nitrogen in the Great Basin U.S.

Vegetation changes associated with climate shifts and anthropogenic disturbance can have major impacts on biogeochemical cycling and soils. Much of the Great Basin, U.S. is currently dominated by sagebrush (Artemisia tridentate (Rydb.) Boivin) ecosystems. Sagebrush ecosystems are increasingly influenced by pinyon (Pinus monophylla Torr. & Frém and Pinus edulis Engelm.) and juniper (Juniperus osteosperma Torr. and Juniperus occidentalis Hook.) expansion. Some scientists and policy makers believe that increasing woodland cover in the intermountain western U.S. offers the possibility of increased organic carbon (OC) storage on the landscape; however, little is currently known about the distribution of OC on these landscapes, or the role that nitrogen (N) plays in OC retention. We quantified the relationship between tree cover, belowground OC, and total below ground N in expansion woodlands at 13 sites in Utah, Oregon, Idaho, California, and Nevada, USA. One hundred and twenty nine soil cores were taken using a mechanically driven diamond tipped core drill to a depth of 90 cm. Soil, coarse fragments, and coarse roots were analyzed for OC and total N. Woodland expansion influenced the vertical distribution of root OC by increasing 15-30 cm root OC by 2.6 Mg ha⁻¹ and root N by 0.04 Mg ha⁻¹. Root OC and N increased through the entire profile by 3.8 and 0.06 Mg ha⁻¹ respectively. Woodland expansion influenced the vertical distribution of soil OC by increasing surface soil (0-15 cm) OC by 2.2 Mg ha⁻¹. Woodland expansion also caused a 1.3 Mg ha⁻¹ decrease in coarse fragment associated OC from 75-90 cm. Our data suggests that woodland expansion into sagebrush ecosystems has limited potential to store additional belowground OC, and must be weighed against the risk of increased wildfire and exotic grass invasion.     

Decoding temporal and spatial patterns of fault uplift using transient river long profiles

We present detailed observations of rivers crossing active normal faults in the Central Apennines, Italy, where excellent constraints exist on the temporal and spatial history of fault movement. We demonstrate that rivers with drainage areas > 10 km² and crossing faults that have undergone an increase in throw rate within the last 1 My, have significant long-profile convexities. In contrast, channels that cross faults that have had a constant-slip rate for 3 My have concave-up profiles and have similar concavities and steepness indices to rivers that do not cross any active fault structures. This trend is consistent across the Central Apennines and cannot be explained by appeal to lithology or regional base level change. The data challenge the belief that active faulting must always be reflected in river profiles; instead, the long-profile convexities are best explained as a transient response of the river system to a change in tectonic uplift rate. Moreover, for these rivers we demonstrate that the height of the profile convexity, as measured from the fault, scales with the magnitude of the uplift rate increase on the fault; and we establish that this relationship holds for throw rate variation along strike for the same fault segment, as well as between faults. These findings are shown to be consistent with predictions of channel response to changing uplift rate rates using a detachment-limited fluvial erosion model, and they illustrate that analysis of the magnitude of profile convexities has considerable predictive potential for extracting tectonic information. We also demonstrate that the migration rate of the profile convexities varies from 1.5–10 mm/y, and is a function of the slip rate increase as well as the drainage area. This is consistent with n > 1 for the slope exponent in a classical detachment-limited stream-power erosion law, but could potentially be explained by incorporating an erosion threshold or an explicit role for sediment in enhancing erosion rates. Finally, we show that for rivers in extensional settings, where the response times to tectonic perturbation are long (in this case > 1 My), attempts to extract tectonic uplift rates from normalised steepness indices are likely to be flawed because topographic steady state has not yet been achieved.     

A model for post‐orogenic development of a mountain range and its foreland

Decaying mountain ranges often show a surprisingly dynamic pattern of landscape evolution. Although one might expect a simple, monotonic decline in relief over time, evidence from several inactive mountain ranges shows alternating sequences of deposition and erosion in the associated basins, suggesting variations in relief and exhumation rate in the ranges themselves. Examples include the Southern Rocky Mountains, the Pyrenees, the European Alps and the Atlas Mountains. In this paper, we explore the possible origins of post‐orogenic landscape dynamics using a simple mathematical model of a mountain range and an adjacent foreland basin. The analysis highlights the importance of mass balance. In particular, a switch from basin exhumation to renewed sedimentation requires either an increase in sediment influx from the range or a decrease in sediment outflux beyond the basin margin. Although it is widely understood that post‐orogenic changes in erosion and sediment flux can have multiple causes (including climate change, regional tectonic uplift or tilting, or exhumation of variable lithologies), an important implication of our analysis is that the impact of such changes must differ in sign or magnitude between the range and the basin to be recorded. This requirement places an important constraint on viable explanations for alternating sequences of deposition and erosion in a decaying mountain‐basin pair.     

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.     

Patterns of rock fragment cover generated by tillage erosion

Intensively cultivated areas in the upper part of the Guadalentin catchment (southeast Spain) show a systematic spatial pattern of surface rock fragment cover (R_c). The objective of this paper is to quantify and to explain this spatial rock fragment cover pattern. Therefore, a map of an intensively cultivated area of 5 km² was digitised, and for each pixel total topographic curvature was calculated. Next, rock fragment cover was determined photographically at 35 sites with a range of total slope curvatures. A linear relation between total curvature and rock fragment cover was found, except for narrow concavities. It was hypothesised that this pattern can be explained by a significant net downslope movement of rock fragments and fine earth by tillage. The displacement distances of rock fragments by tillage with a duckfoot chisel were measured by monitoring the displacement of tracers (painted rock fragments and aluminium cubes) on 5 sites having different slopes. The rare of tillage erosion for one tillage pass with a duckfoot chisel, expressed by the diffusion constant (k), equals 282 kg/m for up and downslope tillage and only 139 kg/m for contour tillage. Nomograms indicate that mean denudation rates in almond groves due to tillage erosion (3 to 5 tillage passes per year) can easily amount to 1.5–2.6 mm/year for contour tillage and up to 3.6–5.9 mm/year for up- and downslope tillage for a field, 50 m long and having a slope of 20%. These figures are at least one order of magnitude larger than reported denudation rates caused by water erosion in similar environments. Hence tillage erosion contributes significantly to land degradation. The downslope soil flux induced by tillage not only causes considerable denudation on topographic convexities (hill tops and spurs) and upper field boundaries but also an important sediment accumulation in topographic concavities (hollows and valley bottoms) and at lower field boundaries. Kinetic sieving (i.e. the upward migration of rock fragments) by the tines of the duckfoot chisel also concentrates the largest rock fragments in the topsoil in such a way that a rock fragment mulch develops in narrow valleys and at the foot of the slopes. These results clearly indicate that tillage erosion is the main process responsible for the observed rock fragment cover pattern in the study area. Since the study area is representative for many parts of southern Spain where almond groves have expanded since 1970, the results have a wider application. They show to what extent intensive tillage of steep slopes has contributed to the increase in soil degradation, to changes in hillslope morphology (i.e. strong denudation of convexities, development of lynchets and rapid infilling of narrow valley bottoms) and to the development of rock fragment cover patterns which control the spatial variability of the hydrological and water erosion response within such landscapes.     

Regional relief characteristics and denudation pattern of the western Southern Alps, New Zealand

The Southern Alps of New Zealand are the topographic expression of active oblique continental convergence of the Australian and Pacific plates. Despite inferred high rates of tectonic and climatic forcing, the pattern of differential uplift and erosion remains uncertain. We use a 25-m DEM to conduct a regional-scale relief analysis of a 250-km long strip of the western Southern Alps (WSA). We present a preliminary map of regional erosion and denudation by overlaying mean basin relief, a modelled stream-power erosion index, river incision rates, historic landslide denudation rates, and landslide density. The interplay between strong tectonic and climatic forcing has led to relief production that locally attains 2 km in major catchments, with mean values of 0.65–0.68 km. Interpolation between elevations of major catchment divides indicates potential removal of l0¹–10³ km³, or a mean basin relief of 0.51–0.85 km in the larger catchments. Local relief and inferred river incision rates into bedrock are highest about 50–67% of the distance between the Alpine fault and the main divide. The mean regional relief variability is ± 0.5 km.Local relief, valley cross-sectional area, and catchment width correlate moderately with catchment area, and also reach maximum values between the range front and the divide. Hypsometric integrals show scale dependence, and together with hypsometric curves, are insufficient to clearly differentiate between glacial and fluvial dominated basins. Mean slope angle in the WSA (ψ = 30°) is lower where major longitudinal valleys and extensive ice cover occur, and may be an insensitive measure of regional relief. Modal slope angle is strikingly uniform throughout the WSA (φ = 38–40°), and may record adjustment to runoff and landsliding. Both ψ and φ show non-linear relationships with elevation, which we attribute to dominant geomorphic process domains, such as fluvial processes in low-altitude valley trains, surface runoff and frequent landsliding on montane hillslopes, “relief dampening” by glaciers, and rock fall/avalanching on steep main-divide slopes.     

Quaternary bedrock erosion and landscape evolution in the Sør Rondane Mountains, East Antarctica: Reevaluating rates and processes

Rates and processes of rock weathering, soil formation, and mountain erosion during the Quaternary were evaluated in an inland Antarctic cold desert. The fieldwork involved investigations of weathering features and soil profiles for different stages after deglaciation. Laboratory analyses addressed chemistry of rock coatings and soils, as well as ¹⁰Be and ²⁶Al exposure ages of the bedrock. Less resistant gneiss bedrock exposed over 1 Ma shows stone pavements underlain by in situ produced silty soils thinner than 40 cm and rich in sulfates, which reflect the active layer thickness, the absence of cryoturbation, and the predominance of salt weathering. During the same exposure period, more resistant granite bedrock has undergone long-lasting cavernous weathering that produces rootless mushroom-like boulders with a strongly Fe-oxidized coating. The red coating protects the upper surface from weathering while very slow microcracking progresses by the growth of sulfates. Geomorphological evidence and cosmogenic exposure ages combine to provide contrasting average erosion rates. No erosion during the Quaternary is suggested by a striated roche moutonnée exposed more than 2 Ma ago. Differential erosion between granite and gneiss suggests a significant lowering rate of desert pavements in excess of 10 m Ma^− 1. The landscape has been (on the whole) stable, but the erosion rate varies spatially according to microclimate, geology, and surface composition.     

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.