Coupled numerical–analytical approach to landscape evolution modeling

The Earth's topography is shaped by surface processes that operate on various scales. In particular, river processes control landscape dynamics over large length scales, whereas hillslope processes control the dynamics over smaller length scales. This scale separation challenges numerical treatments of landscape evolution that use space discretization. Large grid spacing cannot account for the dynamics of water divides that control drainage area competition, and erosion rate and slope distribution. Small grid spacing that properly accounts for divide dynamics is computationally inefficient when studying large domains. Here we propose a new approach for landscape evolution modeling that couples irregular grid‐based numerical solutions for the large‐scale fluvial dynamics and continuum‐based analytical solutions for the small‐scale fluvial and hillslope dynamics. The new approach is implemented in the landscape evolution model DAC (divide and capture). The geometrical and topological characteristics of DAC's landscapes show compatibility with those of natural landscapes. A comparative study shows that, even with large grid spacing, DAC predictions fit well an analytical solution for divide migration in the presence of horizontal advection of topography. In addition, DAC is used to study some outstanding problems in landscape evolution. (i) The time to steady‐state is investigated and simulations show that steady‐state requires much more time to achieve than predicted by fixed area calculations, due to divides migration and persistent reorganization of low‐order streams. (ii) Large‐scale stream captures in a strike‐slip environment are studied and show a distinct pattern of erosion rates that can be used to identify recent capture events. (iii) Three tectono‐climatic mechanisms that can lead to asymmetric mountains are studied. Each of the mechanisms produces a distinct morphology and erosion rate distribution. Application to the Southern Alps of New Zealand suggests that tectonic advection, precipitation gradients and non‐uniform tectonic uplift act together to shape the first‐order topography of this mountain range. Copyright © 2014 John Wiley & Sons, Ltd.     

Complexities of landscape evolution during incision through layered stratigraphy with contrasts in rock strength

Variation in the erodibility of rock units has long been recognized as an important determinant of landscape evolution but has been little studied in landscape evolution models. We use a modified version of the Channel‐Hillslope Integrated Landscape Development (CHILD) model, which explicitly allows for variations in rock strength, to reveal and explore the remarkably rich, complex behavior induced by rock erodibility variations in even very simple geologic settings with invariant climate and tectonics. We study the importance of relative contrasts in erodibility between just two units, the order of these units (whether hard rocks overlie soft or soft rocks overlie hard) and the orientation of the contact between the two units. We emphasize the spatial and temporal evolution of erosion rates, which have important implications for basin analysis, detrital mineral records, and the interpretation of cosmogenic isotope concentrations in detrital samples. Results of the landscape evolution modeling indicate that the stratigraphic order of units in terms of erodibility, the gross orientation of the contact (i.e. dipping away or toward the outlet of the landscape) and the contact dip angle all have measurable effects on landscape evolution, including significant spatial and temporal variations in erosion rates. Steady‐state denudation conditions are unlikely to develop in landscapes with significant contrasts in rock strength in horizontal to moderately tilted rock layers, at least at the scale of the entire landscape. Additionally, our results demonstrate that there is no general relation between rock erodibility and erosion rates in natural settings. Although rock erodibility directly controls the erosion rate constant in our models, it is not uncommon for higher erosion rates to occur in the harder, less erodible rock. Indeed erosion rates may be either greater or less than the rock uplift rate (invariant in time and space in our models) in both hard and soft rocks, depending on the local geology, topography, and the pattern of landscape evolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Erosion of an active fault scarp leads to drainage capture in the Amazon region,Brazil

Far from the continental margin, drainage basins in Central Amazonia should be in topographic steady state; but they are not. Abandoned remnant fluvial valleys up to hundreds of square kilometers in size are observed throughout Amazonia, and are evidence of significant landscape reorganization. While major Late Miocene drainage shifts occurred due to initiation of the transcontinental Amazon River, local landscape change has remained active until today. Driven either by dynamic topography, tectonism, and/or climatic fluctuations, drainage captures in Amazonia provide a natural experiment for assessing the geomorphic response of low‐slope basins to sudden, capture related base‐level falls. This paper evaluates the timing of geomorphic change by examining a drainage capture event across the Baependi fault scarp involving the Cuieiras and Tarumã‐Mirim River basins northwest of the city of Manaus in Brazil. A system of capture‐related knickpoints was generated by base‐level fall following drainage capture; through numerical modeling of their initiation and propagation, the capture event is inferred to have occurred between the middle and late Pleistocene, consistent with other studies of landscape change in surrounding areas. In low‐slope settings like the Amazon River basin, base‐level fall can increase erosion rates by more than an order of magnitude, and moderate to large river basins can respond to episodes of base‐level fall over timescales of tens to hundreds of thousands of years. Copyright © 2013 John Wiley & Sons, Ltd.     

Landscape Evolution Modelling of naturally dammed rivers

Natural damming of upland river systems, such as landslide or lava damming, occurs worldwide. Many dams fail shortly after their creation, while other dams are long‐lived and therefore have a long‐term impact on fluvial and landscape evolution. This long‐term impact is still poorly understood and landscape evolution modelling (LEM) can increase our understanding of different aspects of this response. Our objective was to simulate fluvial response to damming, by monitoring sediment redistribution and river profile evolution for a range of geomorphic settings. We used LEM LAPSUS, which calculates runoff erosion and deposition and can deal with non‐spurious sinks, such as dam‐impounded areas. Because fluvial dynamics under detachment‐limited and transport‐limited conditions are different, we mimicked these conditions using low and high erodibility settings, respectively. To compare the relative impact of different dam types, we evaluated five scenarios for each landscape condition: one scenario without a dam and four scenarios with dams of increasing erodibility. Results showed that dam‐related sediment storage persisted at least until 15 000 years for all dam scenarios. Incision and knickpoint retreat occurred faster in the detachment‐limited landscape than in the transport‐limited landscape. Furthermore, in the transport‐limited landscape, knickpoint persistence decreased with increasing dam erodibility. Stream capture occurred only in the transport‐limited landscape due to a persisting floodplain behind the dam and headward erosion of adjacent channels. Changes in sediment yield variation due to stream captures did occur but cannot be distinguished from other changes in variation of sediment yield. Comparison of the model results with field examples indicates that the model reproduces several key phenomena of damming response in both transport‐limited and detachment‐limited landscapes. We conclude that a damming event which occurred 15 000 years ago can influence present‐day sediment yield, profile evolution and stream patterns. Copyright © 2014 John Wiley & Sons, Ltd.     

Butte detachment: how pre‐rift geological structure and drainage integration drive escarpment evolution at rifted continental margins

The erosional pattern of passive margins often follows the fabric of ancient, compressional geological structures exposed by the topographic energy of rifting. As erosion cuts into these belted outcrop systems they impose initial and boundary conditions that steer drainage recession into the plateau edge and control escarpment‐forming conditions. Pattern therefore controls process. Although generic surface process models predict scarp patterns and retreat in settings devoid of geological heterogeneity, they tend to do so only at isolated locations and for periods shorter than the lifespan of the escarpments. Thus, to focus on relatively narrow strike‐perpendicular swaths of passive margin topography misses important aspects of drainage integration, which involves mobile drainage basin boundaries shifting across but also along the strike of inherited geological structures and through continental‐scale bioclimatic zones. Space‐for‐time substitution along three passive margin escarpments (Blue Ridge, Western Ghats, Eastern Ghats) reveals the significance of escarpment jumps and the detachment of topographic outliers, here generically termed ‘buttes’, as key processes of escarpment evolution. The examples show that these continental escarpments are strongly patterned after pre‐rift structural and lithological heterogeneities. As seaward sloping drainages cut into the rift margin, they extend their drainage heads in a non‐uniform and unsteady fashion. As a result escarpments can form, be destroyed, reform, and leave topographic vestiges (buttes) of the retreating escarpment. Given the pre‐rift geological heterogeneities, there are no a priori reasons why escarpment landscape change should be uniform, steady or self‐similar. Copyright © 2010 John Wiley & Sons, Ltd.     

Dry aeolian dust accumulation in rocky deserts: ­a medium‐term field experiment based on short‐term wind tunnel simulations

The spatial pattern of medium‐term (a few months) dry aeolian dust accumulation in rocky deserts is predicted using short‐term deposition and erosion experiments in a wind tunnel. The predictions are tested in a field experiment set up in the northern Negev Desert of Israel. The results show that superimposing wind tunnel deposition and erosion maps usually leads to correct predictions of medium‐term dust accumulation. The predictions are somewhat less confident near the inflection lines of windward hillslopes, where small‐scale irregularities in the local topography make it difficult to locate the exact position of the areas of little accumulation. Elsewhere in the topography predictions are good, and the method works satisfactorily. Highest accumulation occurs on concave windward slopes and, to a lesser extent, on slopes parallel to the wind. Little accumulation occurs on the convex windward slopes and in dust separation bubbles. The smallest accumulation rates are observed immediately upwind of the top of pronounced hills and on leeslopes. The rate of dry dust accumulation measured during the field experiment varied from 17 to 93 g m⁻² a⁻¹, depending on the topographic position of the accumulation plots. For most plots, it was of the order of 30–60 g m⁻² a⁻¹. Copyright © 2000 John Wiley & Sons, Ltd.     

Landscape response to tectonic and climatic forcing in the foredeep of the southern Apennines,Italy: insights from Quaternary stratigraphy,quantitative geomorphic analysis,and denudation rate proxies

We present new data about the morphological and stratigraphic evolution and the rates of fluvial denudation of the Tavoliere di Puglia plain, a low‐relief landscape representing the northernmost sector of the Pliocene‐Pleistocene foredeep of the southern Apennines. The study area is located between the easternmost part of the southern Apennine chain and the Gargano promontory and it is characterized by several orders of terraced fluvial deposits, disconformably overlying lower Pleistocene marine clay and organized in a staircase geometry, which recorded the emersion and the long‐term incision history of this sector since mid‐Pleistocene times. We used the spatial and altimetric distribution of several orders of middle to late Pleistocene fluvial terraces in order to perform paleotopographic reconstruction and GIS‐aided eroded volumes estimates. Then, we estimated denudation rates on the basis of the terraces chronostratigraphy, supported by published OSL and AAR dating. Middle to upper Pleistocene denudation rates estimated by means of such an approach are slightly lower than 0.1 mm yr^‐1, in good agreement with short‐term data from direct and indirect evaluation of suspended sediment yield. The analysis of longitudinal river profiles using the stream power erosion model provided additional information on the incision rates of the studied area. Middle to late Quaternary uplift rates (about 0.15 mm yr^‐1), calculated on the basis of the elevation above sea level of marine deposits outcropping in the easternmost sector of the study area, are quite similar to the erosion rates average value, thus suggesting a steady‐state fluvial incision. The approach adopted in this work has demonstrated that erosion rates traditionally obtained by quantitative geomorphic analysis and k_sn estimations can be successfully integrated to quantify rates of tectonic or geomorphological processes of a landscape approaching steady‐state equilibrium. Copyright © 2014 John Wiley & Sons, Ltd.     

Linking landslides,hillslope erosion,and landscape evolution

This Virtual Issue highlights 10 recent innovative, unconventional, or otherwise significant contributions to Earth Surface Processes and Landforms that help advance the state‐of‐the‐art in research on linkages between landslides, hillslope erosion, and landscape evolution. The selected studies address this feedback within a temporal spectrum that ranges from the event to the millennial scale, thus underscoring the importance of detailed field observations, high‐resolution digital topographic data, and geochronological methods for increasing our capability of quantifying landslide processes and hillslope erosion. Copyright © 2009 John Wiley & Sons, Ltd.     

Soil formation and weathering in a permafrost environment of the Swiss Alps: a multi‐parameter and non‐steady‐state approach

Spatially discontinuous permafrost conditions frequently occur in the European Alps. How soils under such conditions have evolved and how they may react to climate warming is largely unknown. This study focuses on the comparison of nearby soils that are characterised by the presence or absence of permafrost (active‐layer thickness: 2–3 m) in the alpine (tundra) and subalpine (forest) range of the Eastern Swiss Alps using a multi‐method (geochemical and mineralogical) approach. Moreover, a new non‐steady‐state concept was applied to determine rates of chemical weathering, soil erosion, soil formation, soil denudation, and soil production. Long‐term chemical weathering rates, soil formation and erosion rates were assessed by using immobile elements, fine‐earth stocks and meteoric ¹⁰Be. In addition, the weathering index (K + Ca)/Ti, the amount of Fe‐ and Al‐oxyhydroxides and clay minerals characteristics were considered. All methods indicated that the differences between permafrost‐affected and non‐permafrost‐affected soils were small. Furthermore, the soils did not uniformly differ in their weathering behaviour. A tendency towards less intense weathering in soils that were affected by permafrost was noted: at most sites, weathering rates, the proportion of oxyhydroxides and the weathering stage of clay minerals were lower in permafrost soils. In part, erosion rates were higher at the permafrost sites and accounted for 79–97% of the denudation rates. In general, soil formation rates (8.8–86.7 t/km²/yr) were in the expected range for Alpine soils. Independent of permafrost conditions, it seems that the local microenvironment (particularly vegetation and subsequently soil organic matter) has strongly influenced denudation rates. As the climate has varied since the beginning of soil evolution, the conditions for soil formation and weathering were not stable over time. Soil evolution in high Alpine settings is complex owing to, among others, spatio‐temporal variations of permafrost conditions and thus climate. This makes predictions of future behaviour very difficult. Copyright © 2016 John Wiley & Sons, Ltd.     

The shape of patterns to come: from initial formation to long‐term evolution

Many studies focus on the emergence and development of rhythmic landscape patterns. In this contribution we explore the different behaviors found as patterns evolve; the trajectories that patterns exhibit as they transit from infinitesimal‐amplitude perturbation to a statistically steady state (or in some cases to continued statistical evolution). The variety of behaviors observed, either through field and laboratory experiments or numerical modeling, can be reduced to four classes: (a) simple stabilization where predictions based on the initial growth of small perturbations corresponds with the characteristics of patterns observed in nature; (b) significant pattern coarsening en route to saturated wavelength, where non‐linear interactions between finite‐amplitude pattern elements change the geometric properties of a pattern as it approaches steady‐state; (c) perpetual coarsening where the wavelength associated with the emerging pattern continues to grow over time and is only limited by physical boundaries or external constrains; (d) slow evolution toward a different attractor, a novel behavior observed in numerical modeling that involves profound temporal changes in pattern characteristics. Within these classes we also observe generalizable non‐linear behaviors: dependence on initial conditions, the emergence of pattern‐scale variables such as pattern defects, and the presence of multiple stable states. Predicting the shape of patterns to come remains a challenge – one that we suggest requires a range of modeling approaches to address both initial instabilities and the emergent properties of evolving patterns, which involve disparate forms of non‐linear interactions. Consideration of generic system behaviors at the pattern scale could enhance future pattern formation studies, facilitating appropriate pairings of analysis approaches and pattern‐evolution modes. Copyright © 2013 John Wiley & Sons, Ltd.     

Transient hillslope erosion in slow evolution landscapes

Transient evolution and adjustment to changing tectonic and climatic boundary conditions is an essential attribute of landscapes, and characterizing transient behavior is a key to understanding their dynamics and history. Developing new approaches to detect such transience has been explored by various methods, in particular to identify landscape response to Late Cenozoic and Quaternary climatic changes. Such studies have often focused on regions of high relief and/or active tectonic activity where interferences between tectonic and climatic signals might complicate the interpretation of the observations. We investigated the case of the hillslopes of the Serra do Cipó quartzitic range in SE Brazil in order to detect and quantify transience in a tectonically quiescent landscape over 100-ka timescales. We determined hilltop curvature from a high-resolution digital surface model derived from Pléiades imagery and measured cosmogenic nuclide (¹⁰Be and ²⁶Al) concentrations at these hilltop sites. We compare both observations with predictions of hillslope diffusion theory, observing a distinctive signature of an acceleration of denudation. We performed a joint inversion of topographic and isotopic data to retrieve an evolution of the hillslope sediment transport coefficient through time. The timing of the increase in denudation cannot be unequivocally associated with a single climatic event but is consistent with important, climatically modulated fluctuations in precipitation and erosion in this area during the Middle and Late Pleistocene.     

Plateau reduction by drainage divide migration in the Eastern Cordillera of Colombia defined by morphometry and 10Be terrestrial cosmogenic nuclides

Catchment‐wide erosion rates were defined using ¹⁰Be terrestrial cosmogenic nuclides for the Eastern Cordillera of the Colombian Andes to help determine the nature of drainage development and landscape evolution. The Eastern Cordillera, characterized by a smooth axial plateau bordered by steep flanks, has a mean erosion rate of 11 ± 1 mm/ka across the plateau and 70 ± 10 mm/ka on its flanks, with local high rates >400 mm/ka. The erosional contrast between the plateau and its flanks was produced by the increase in the orogen regional slope, derived from the progressive shortening and thickening of the Eastern Cordillera. The erosion rates together with digital topographic analysis show that the drainage network is dynamic and confirms the view that drainage divides in the Eastern Cordillera are migrating towards the interior of the mountain belt resulting in progressive drainage reorganization from longitudinal to transverse‐dominated rivers and areal reduction of the Sabana de Bogotá plateau. Copyright © 2016 John Wiley & Sons, Ltd.     

Test the topographic steady state in an active mountain belt

Landscape evolution is the window to the link between deep earth and surface processes. One of the key issues in landscape evolution is to characterize the topographic steady state of mountain belts. The Taiwan mountain belt is an extraordinary case due to its extremely high uplift and denudation rates. The uplift of Taiwan Island is caused by the oblique collision between the Luzon Arc and the East Asian continent. In this case, the mountain building process in the north always occurs earlier than that in the south,which causes the spatial distribution of steady-state regions.The East Central Range receives much research attention with the presence of river basins that mainly distribute along the trajectory of the collision propagation. Normally, based on analyses of geomorphic parameters, the whole Central Range, or at least part of it, should be at a topographic steady state. However, the balance between uplift rates and denudation rates that exist in these regions is seldomly tested. In this contribution, we make a comprehensive literature review on the uplift and denudation rates derived from various approaches, including sediment yields, in-situ cosmogenic nuclide ~(10) Be, incision of river channels,thermochronology, and GPS observations. This literature review reveals that the topographic steady state may prevail in the northern and middle parts of the East Central Range.However, an obvious inconsistency in denudation rates calculated by different methods prevents us from better constraining the topographic steady state in some regions of this mountain range.     

Modelling long‐term (300 ka) upland catchment response to multiple lava damming events

Landscapes respond in complex ways to external drivers such as base level change due to damming events. In this study, landscape evolution modelling was used to understand and analyse long‐term catchment response to lava damming events. PalaeoDEM reconstruction of a small Turkish catchment (45 km²) that endured multiple lava damming events in the past 300 ka, was used to derive long‐term net erosion rates. These erosion rates were used for parameter calibration and led to a best fit parameter set. This optimal parameter set was used to compare net erosion landscape time series of four scenarios: (i) no uplift and no damming events; (ii) no uplift and three damming events; (iii) uplift and no damming events; and (iv) uplift and three damming events. Spatial evolution of net erosion and sediment storage of scenario (iii) and (iv) were compared. Simulation results demonstrate net erosion differences after 250 000 years between scenarios with and without dams. Initially, trunk gullies show less net erosion in the scenario with damming events compared with the scenario without damming events. This effect of dampened erosion migrates upstream to smaller gullies and local slopes. Finally, an intrinsic incision pulse in the dam scenario results in a higher net erosion of trunk gullies while decoupled local slopes are still responding to the pre‐incision landscape conditions. Sediment storage differences also occur on a 100 ka scale. These differences behaved in a complex manner owing to different timings of the migration of erosion and sediment waves along the gullies for each scenario. Although the specific spatial and temporal sequence of erosion and deposition events is sensitive to local parameters, this model study shows the manner in which past short‐lived events like lava dams have long‐lasting effects on catchment evolution. Copyright © 2014 John Wiley & Sons, Ltd.     

Drainage basin sensitivity to tectonic and climatic forcing: implications of a stochastic model for the role of entrainment and erosion thresholds

Long‐term average rates of channel erosion and sediment transport depend on the frequency–magnitude characteristics of ?ood ?ows that exceed an erosion threshold. Using a Poisson model for rainfall and runoff, analytical solutions are developed for average rates of stream incision and sediment transport in the presence of such a threshold. Solutions are derived and numerically tested for three erosion/transport formulas: the Howard–Kerby shear‐stress incision model, the Bridge–Dominic sediment transport model, and a generic shear‐stress sediment transport model. Results imply that non‐linearity resulting from threshold effects can have a ?rst‐order impact on topography and patterns of dynamic response to tectonic and climate forcing. This non‐linearity becomes signi?cant when fewer than about half of ?ood events are capable of detaching rock or sediment. Predicted morphology and uplift‐gradient scaling is more closely consistent with observations and laboratory experiments than conventional slope‐linear or shear‐linear erosion laws. These results imply that particle detachment thresholds are not details that can be conveniently ignored in long‐term landscape evolution models. Copyright © 2004 John Wiley & Sons, Ltd.     

Slope–channel coupling between pipes,gullies and tributary channels in the Mocatán catchment badlands,Southeast Spain

The dispersive nature of the highly sodic silts of the Triassic‐rich unit of the Góchar formation plays a fundamental role in the erosion of the Mocatán catchment badlands in Almería, where a rejuvenating pipe and incised channel network is rapidly evacuating slope materials. Referring to concepts of medium‐ and long‐term landscape evolution, and incorporating contemporary thoughts on the role of connectivity in system evolution from the geomorphological literature, this paper attempts to develop a conceptual model of the way geologic, topographic, material property and ecological factors combine to explain the complex geomorphological evolution of the site. An electronic distance measurement (EDM) survey was undertaken using a Leica TC3100, to produce a detailed topographic map. This database was supplemented by geomorphological, geological and ecological data derived from ground survey and remote sensing, and further morphometric analysis undertaken. Preferred orientations of channel segments, and the topographic distribution of pipe‐roof‐collapse features and outfalls in relation to known stratigraphic controls, suggests that, once coupled to the slope‐base channel, pipe networks develop in a systematic, sequential way. A wave of incision moving up the main channel reconnects channels with slopes, and the resulting increased hydraulic gradients on sideslopes encourage extensive deep pipe development for the first time in these materials. Once major pipe development is possible, three‐dimensional pipe networks enlarge and then collapse to form an extensive, partially coupled steep‐sided gully network. From this perspective, the coupling of the pipe to a rejuvenating channel is a significant intrinsic threshold event and the main reason that badlands have developed locally in these dispersive materials. It is concluded that erosion in this landscape will only be suppressed after a considerable period of slope‐base stability, which could allow a gradual loss of sodium from the surface by leaching or organic exchanges. Both theoretical and management implications are explored. Copyright © 2007 John Wiley & Sons, Ltd.     

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

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

The topographic signature of a major typhoon

In August 2009, the typhoon Morakot, characterized by a cumulative rainfall up to 2884 mm in about three days, triggered thousands of landslides in Taiwan. The availability of LiDAR surveys before (2005) and after (2010) this event offers a unique opportunity to investigate the topographic signatures of a major typhoon. The analysis considers the comparison of slope–area relationships derived by LiDAR digital terrain models (DTMs). This approach has been successfully used to distinguish hillslope from channelized processes, as a basis to develop landscape evolution models and theories, and understand the linkages between landscape morphology and tectonics, climate, and geology. We considered six catchments affected by a different degree of erosion: three affected by shallow and deep‐seated landslides, and three not affected by erosion. For each of these catchments, 2 m DTMs were derived from LiDAR data. The scaling regimes of local slope versus drainage area suggested that for the catchments affected by landslides: (i) the hillslope‐to‐valley transitions morphology, for a given value of drainage area, is shifted towards higher value of slopes, thus indicating a likely migration of the channelized processes and erosion toward the catchment boundary (the catchment head becomes steeper because of erosion); (ii) the topographic gradient along valley profiles tends to decrease progressively (the valley profile becomes gentler because of sediment deposition after the typhoon). The catchments without any landslides present a statistically indistinguishable slope–area scaling regime. These results are interesting since for the first time, using multi‐temporal high‐resolution topography derived by LiDAR, we demonstrated that a single climate event is able to cause significant major geomorphic changes on the landscape, detectable using slope–area scaling analysis. This provides new insights about landscape evolution under major climate forcing. Copyright © 2015 John Wiley & Sons, Ltd.     

EVALUATION OF TOPOGRAPHIC INDICES FOR EPHEMERAL-GULLY EROSION ASSESSMENT

1 INTRODUCTION Erosion caused by ephemeral flows is a frequent phenomenon in nature and contributes to the shape of the landscape. This type of erosion may cause great soil losses in agricultural areas, which are quickly transferred to the watershed outlets through the rill and gully network (Bennett et al., 2000; Poesen et al., 2003). Concentrated flow erosion is controlled by the erodibility of surface materials, climate, soil use and management, and watershed topography. Several metho…