Time stability and variability of Electronically Scanned Thinned Array Radiometer soil moisture during Southern Great Plains hydrology experiments

Variability and time‐stability analysis for field‐scale (800 m) Electronically Scanned Thinned Array Radiometer soil moisture within a satellite scale footprint (∼ 50 km) were quantified using observations from the Southern Great Plains Hydrology Experiment 1997 and 1999 (SGP97 and SGP99). The pixels' time‐stability properties were examined with respect to soil, vegetation and topographic parameters in order to determine which physical parameters can be used to identify good candidate observation locations for validating soil moisture from satellite observations and global‐scale model output. The results show that the time‐stability concept remains valid at the satellite scale. The root mean square error values were 1·47, 1·51, 1·93 and 2·32% for the 1st, 2nd, 50th and 100th most stable fields, respectively. The most stable locations had sand and clay percentages consistent with sandy loam soils and moderate to high normalized difference vegetation index values. Neither land cover nor topography properties could be used to identify potentially stable fields in the study region. Copyright © 2010 John Wiley & Sons, Ltd.     

Discriminating residual and colluvial soils using topographic analysis in a small,forested catchment in Japan

Soil moisture plays an important role in hydrology. Understanding factors (such as topography, vegetation, and meteorological conditions) that influence spatio‐temporal variability in soil moisture, and how this influence is manifested, is important for understanding hydrological processes. A number of distributed (quasi‐)physical hydrological models have been developed to investigate this subject. Previous studies have shown that the spatial differences in the distribution of soil types (residual and colluvial soils) dominantly reflect spatio‐temporal fluctuations in soil moisture and runoff. We present a methodology for assessing the spatial distribution of residual and colluvial soils, which differ with respect to their physical characteristics, in a 0·88 km² forested catchment with complex topography and a complex land‐use history. Our method is based on penetration resistance profile data; in this data set, each data point represents soil physical characteristics within an area of about 25 m². If the spatial distribution of soils under similar meteorological, geological, historical land use, and other conditions could be characterized on the basis of similarity in topographic features, then the spatial distribution of soil could be predicted based on relationships between various topographic indices (e.g. topographic index and local slope). We tested whether our model correctly assessed the reference data. The model's results were 90·5% correct for residual soils and 87·3% correct for colluvial soils. Further studies will quantify the relationships between topographic features of land covered by residual and colluvial soils and changes in spatio‐temporal variations in the catchment (e.g. vegetation and land use) as a function of geology or meteorology. Copyright © 2007 John Wiley & Sons, Ltd.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

Two‐dimensional hydraulic flood modelling using a finite‐element mesh decomposed according to vegetation and topographic features derived from airborne scanning laser altimetry

Airborne scanning laser altimetry (LiDAR) is an important new data source that can provide two‐dimensional river flood models with spatially distributed floodplain topography for model bathymetry, together with vegetation heights for parameterization of model friction. Methods are described for improving such models by decomposing the model's finite‐element mesh to reflect floodplain vegetation features such as hedges and trees having different frictional properties to their surroundings, and significant floodplain topographic features having high height curvatures. The decomposition is achieved using an image segmentation system that converts the LiDAR height image into separate images of surface topography and vegetation height at each point. The vegetation height map is used to estimate a friction factor at each mesh node. The spatially distributed friction model has the advantage that it is physically based, and removes the need for a model calibration exercise in which free parameters specifying friction in the channel and floodplain are adjusted to achieve best fit between modelled and observed flood extents. The scheme was tested in a modelling study of a flood that occurred on the River Severn, UK, in 1998. A satellite synthetic aperture radar image of flood extent was used to validate the model predictions. The simulated hydraulics using the decomposed mesh gave a better representation of the observed flood extent than the more simplistic but computationally efficient approach of sampling topography and vegetation friction factors on to larger floodplain elements in an undecomposed mesh, as well as the traditional approach using no LiDAR‐derived data but simply using a constant floodplain friction factor. Use of the decomposed mesh also allowed velocity variations to be predicted in the neighbourhood of vegetation features such as hedges. These variations could be of use in predicting localized erosion and deposition patterns that might result in the event of a flood. Copyright © 2003 John Wiley & Sons, Ltd.     

A new model to analyse the impact of woody riparian vegetation on the geotechnical stability of riverbanks

We present a geotechnical stability analysis for the planar failure of riverbanks, which incorporates the effects of root reinforcement and surcharge for mature stands of woody riparian vegetation. The analysis relies on a new method of representing the root distribution in the soil, which evaluates the effects of the vegetation's position on the bank. The model is used in a series of sensitivity analyses performed for a wide range of bank morphological (bank slope and height) and sedimentological (bank cohesion and friction angle) conditions, enabling discrimination of the types of bank environment for which vegetation has an effect on bank stability. The results indicate that woody vegetation elements have a maximal impact on bank stability when they are located at the ends of the incipient failure plane (i.e. at the bank toe or at the intersection of the failure plane with the floodplain) and that vegetation has a greater effect on net bank stability when it is growing on low, shallow, banks comprised of weakly cohesive sediments. However, the magnitude of these effects is limited, with vegetation typically inducing changes (relative to non‐vegetated banks) in simulated factors of safety of less than 5%. If correct, this suggests that the well documented effects of vegetation on channel morphology must be related to alternative process mechanisms (such as the interaction of vegetation with river flows) rather than the mechanical effects of vegetation on bank failure, except in special cases where the equivalent non‐vegetated bank has a highly marginal stability status. Copyright © 2007 John Wiley & Sons, Ltd.     

Testing the influence of topography and material properties on catchment‐scale soil moisture patterns using remotely sensed vegetation patterns in a humid temperate catchment,northern Britain

In order to evaluate the relationship between the apparent complexity of hillslope soil moisture and the emergent patterns of catchment hydrological behaviour and water quality, we need fine‐resolution catchment‐wide data on soil moisture characteristics. This study proposes a methodology whereby vegetation patterns obtained from high‐resolution orthorectified aerial photographs are used as an indicator of soil moisture characteristics. This enables us to examine a set of hypotheses regarding what drives the spatial patterns of soil moisture at the catchment scale (material properties or topography). We find that the pattern of Juncus effusus vegetation is controlled largely by topography and mediated by the catchment's material properties. Characterizing topography using the topographic index adds value to the soil moisture predictions relative to slope or upslope contributing area (UCA). However, these predictions depart from the observed soil moisture patterns at very steep slopes or low UCAs. Copyright © 2012 John Wiley & Sons, Ltd.     

Biotic controls on shallow translational landslides

In undisturbed tropical montane rainforests massive organic layers accommodate the majority of roots and only a small fraction of roots penetrate the mineral soil. We investigated the contribution of vegetation to slope stability in such environments by modifying a standard model for slope stability to include an organic layer with distinct mechanical properties. The importance of individual model parameters was evaluated using detailed measurements of soil and vegetation properties to reproduce the observed depth of 11 shallow landslides in the Andes of southern Ecuador. By distinguishing mineral soil, organic layer and above‐ground biomass, it is shown that in this environment vegetation provides a destabilizing effect mainly due to its contribution to the mass of the organic layer (up to 973 t ha^{? 1} under wet conditions). Sensitivity analysis shows that the destabilizing effect of the mass of soil and vegetation can only be effective on slopes steeper than 37.9°. This situation applies to 36% of the study area. Thus, on the steep slopes of this megadiverse ecosystem, the mass of the growing forest promotes landsliding, which in turn promotes a new cycle of succession. This feedback mechanism is worth consideration in further investigations of the impact of landslides on plant diversity in similar environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the role of climate,vegetation and pedogenesis in shallow translational hillslope failure

Recent research in geomorphology has considered the significance of progressive pedogenesis and climatic change to slope failure initiation for the Holocene, using physically based models. To date, the significance of vegetation change to slope stability has been largely unexplored through modelling, since available physically based models cannot consider vegetation effects directly. To address the existing deficiency this paper adapts, parameterizes and applies a physically based model of slope hydrology and stability to the combined effect of vegetation change and progessive pedogenesis on slope failure initiation. There is considerable debate in the literature concerning the relative significance of climatic change and vegetation modification to slope failure initiation in the Holocene. This paper uses the model to provide additional evidence for situations in which either climatic or vegetation change is significant to slope failure, depending on the prevailing degree of soil development. The results indicate that young podsols appear to be stable under all the climatic and vegetation conditions considered, but mature podsols may be susceptible to failure. Both climate and vegetation influence slope stability, but their relative significance depends on the stage of soil development. In particular, the stability of young soils is influenced considerably by vegetation, while climate assumes greater significance in mature soils. It is recognized that this conclusion is limited to freely draining podsol profiles, and that more research is needed to consider other soil type and vegetation combinations.     

The use of mycorrhiza for eco‐engineering measures in steep alpine environments: effects on soil aggregate formation and fine‐root development

Steep erosion‐prone and vegetation‐free slopes are widespread in alpine areas and are often discussed since they have a high socio‐economic damage potential. We present an eco‐engineering approach to test whether a mycorrhizal inoculum improves the establishment of hedge brush layers and in turn soil structural stability on a steep, coarse‐grained vegetation‐free slope in the eastern Swiss Alps. We established (i) mycorrhizal and (ii) non‐mycorrhizal treated eco‐engineered research plots on a field experimental scale, covering a total area of approximately 1000 m² on an east‐northeast (ENE) exposed slope, where many environmental parameters can be regarded as homogeneous. After a full vegetation period, we quantified soil aggregate stability, the formation of water stable aggregates and the fine‐root development. Our results illustrate that the establishment of brush layers without mycorrhizal inoculum increased aggregate stability significantly. Against our expectation and glasshouse experiments, the addition of mycorrhizal inoculum did not have a statistically significant effect after one vegetation period although it tended to increase aggregate stability. Analogously, root length density (RLD) tended to be higher at the non‐mycorrhizal treated site. Aggregate stability was significantly correlated with RLD. Studies on a bigger field experimental scale are inevitable, complement glasshouse studies and lead to a better understanding for a successful application of sustainable eco‐engineering measures in alpine environments. Based on our results and considering the fact that the response time in natural ecosystems may be slower than in laboratory approaches, we conclude that long‐term field studies are necessary to validate results gained through laboratory experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

Prediction of vegetation anomalies over an inland river basin in north‐western China

Vegetation cover plays an important role in linking the atmosphere, water, and land and is deemed as a key indicator in the terrestrial ecological system. Therefore, it is of great importance to monitor vegetation dynamics and understand the mechanisms of vegetation change, including that driven by climate change. This study examines (a) the evolution of vegetation dynamics over the Heihe River Basin in the typical arid zone in north‐western China using nonparametric Mann–Kendall test and Thiel Sen's slope; (b) the relationships between remotely sensed vegetation indices (normalized difference vegetation index [NDVI] and enhanced vegetation index [EVI]) and hydroclimatic variables based on correlation analysis; and (c) the prediction of vegetation anomalies using a multiple linear regression model. For the analysis, the Moderate Resolution Imaging Spectroradiometer NDVI/EVI product and the gridded daily meteorological data at a spatial resolution of 0.125° over the period 2001–2010 are considered. The results indicate that vegetation cover improved over a large proportion during 2001–2010, with a significant trend towards warm and wet, characterized by an increase in average annual temperature and precipitation by 0.042 °C/year and 5.8 mm/year, respectively. We test the feasibility of NDVI and EVI in quantifying the responses of vegetation anomaly to climate change and develop a statistical model to predict vegetation dynamics in the basin. The NDVI‐based model is found to be more reliable than the EVI‐based model, partly due to the vegetation characteristics and geomorphologic properties of the study region. The proposed model performs well when there is no lag time between meteorological factors and vegetation indices for grassland and cropland, whereas 1‐month lead time prediction is found to be best for forest. The soil water content is introduced as an extra explanatory variable, which effectively improves the prediction accuracy for different land use types. In general, the predictive ability of the proposed model is stable and satisfactory, and the model can provide useful early warning information for regional water resources management under changing climate.     

A channel head stability criterion for first‐order catchments

The present analysis derives a stability criterion for long‐term equilibrium channel heads. The concept of finite perturbation analysis is presented, during which the surface is subjected to perturbations of a finite amplitude and resulting changes in flow path structure and slope are computed. Based on these quantities the analysis predicts whether the perturbed location is going to erode, be filled in or remain steady. The channel head is defined geometrically as the focus point of converging flow lines at the bottom of hollows. It is demonstrated that stability at the channel head grows out of the competition between the rate of flow path convergence and the degree of profile concavity. Analytical functions are derived to compute channel head‐contributing area and ‐slope, flow path convergence and profile concavity as a function of perturbation depth, distance from the crest and the initial slope. In a numerical model these quantities point to the long‐term equilibrium channel head position, which is shown to depend also on the width to length ratio of hollows. It is also demonstrated that the equilibrium channel head position is sensitive to the base‐level lowering/non‐dimensional slope length ratio and to the slope of the initial topography. Morphometrical measurements both in the field and on simulated topographies were used to test the theoretical predictions.     

A hillslope hydrology approach for catchment‐scale slope stability analysis

Regional analysis of slope stability is often constrained by availability of data. Model requirements for input data cannot be met at the desired spatial resolution because data are either site‐speci?c or non‐existent. Faced with these dif?culties it has often been the practice to assume that certain parameters are uniform throughout the area of interest. An alternative approach proposed here allows a more detailed discrimination of slope stability conditions. Based on the principles of hillslope hydrology, hydrologic information can be generated at suf?cient resolution to allow higher resolution slope stability analysis. Measurements from an instrumented network in a small area have been used to establish index‐based models for topographic and climate‐related controls of piezometric response. The ability to relate groundwater levels to rainfall and topographic parameters provides a means of up‐scaling to larger catchments and ultimately the opportunity to generate a catchment‐wide prediction of the distribution, magnitude and frequency of rainstorm‐generated groundwater levels. The example provided in this study uses the topography index of TOPMODEL in GIS to predict the spatial patterns of groundwater elevation for seasonal soil moisture conditions and given rainfall inputs. This allows modelling of catchment‐wide response of soil water to rainstorms with different return periods (representing different magnitudes), and is an essential prerequisite for a probabilistic regional slope stability analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

Which way do you lean? Using slope aspect variations to understand Critical Zone processes and feedbacks

Soil‐mantled pole‐facing hillslopes on Earth tend to be steeper, wetter, and have more vegetation cover compared with adjacent equator‐facing hillslopes. These and other slope aspect controls are often the consequence of feedbacks among hydrologic, ecologic, pedogenic, and geomorphic processes triggered by spatial variations in mean annual insolation. In this paper we review the state of knowledge on slope aspect controls of Critical Zone (CZ) processes using the latitudinal and elevational dependence of topographic asymmetry as a motivating observation. At relatively low latitudes and elevations, pole‐facing hillslopes tend to be steeper. At higher latitudes and elevations this pattern reverses. We reproduce this pattern using an empirical model based on parsimonious functions of latitude, an aridity index, mean‐annual temperature, and slope gradient. Using this empirical model and the literature as guides, we present a conceptual model for the slope‐aspect‐driven CZ feedbacks that generate asymmetry in water‐limited and temperature‐limited end‐member cases. In this conceptual model the dominant factor driving slope aspect differences at relatively low latitudes and elevations is the difference in mean‐annual soil moisture. The dominant factor at higher latitudes and elevations is temperature limitation on vegetation growth. In water‐limited cases, we propose that higher mean‐annual soil moisture on pole‐facing hillslopes drives higher soil production rates, higher water storage potential, more vegetation cover, faster dust deposition, and lower erosional efficiency in a positive feedback. At higher latitudes and elevations, pole‐facing hillslopes tend to have less vegetation cover, greater erosional efficiency, and gentler slopes, thus reversing the pattern of asymmetry found at lower latitudes and elevations. Our conceptual model emphasizes the linkages among short‐ and long‐timescale processes and across CZ sub‐disciplines; it also points to opportunities to further understand how CZ processes interact. We also demonstrate the importance of paleoclimatic conditions and non‐climatic factors in influencing slope aspect variations. Copyright © 2017 John Wiley & Sons, Ltd.     

Do not only connect: a model of infiltration‐excess overland flow based on simulation

The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.     

Floodplain friction parameterization in two‐dimensional river flood models using vegetation heights derived from airborne scanning laser altimetry

Two‐dimensional (2‐D) hydraulic models are currently at the forefront of research into river flood inundation prediction. Airborne scanning laser altimetry is an important new data source that can provide such models with spatially distributed floodplain topography together with vegetation heights for parameterization of model friction. The paper investigates how vegetation height data can be used to realize the currently unexploited potential of 2‐D flood models to specify a friction factor at each node of the finite element model mesh. The only vegetation attribute required in the estimation of floodplain node friction factors is vegetation height. Different sets of flow resistance equations are used to model channel sediment, short vegetation, and tall and intermediate vegetation. The scheme was tested in a modelling study of a flood event that occurred on the River Severn, UK, in October 1998. A synthetic aperture radar image acquired during the flood provided an observed flood extent against which to validate the predicted extent. The modelled flood extent using variable friction was found to agree with the observed extent almost everywhere within the model domain. The variable‐friction model has the considerable advantage that it makes unnecessary the unphysical fitting of floodplain and channel friction factors required in the traditional approach to model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Long‐term modelling of landslides for different forest management practices

Long‐term effects of different forest management practices on landslide initiation and volume were analysed using a physically based slope stability model. The watershed‐based model calculates the effects of multiple harvesting entries on slope stability by accounting for the cumulative impacts of a prior vegetation removal on a more recent removal related to vegetation root strength and tree surcharge. Four sequential clearcuts and partial cuts with variable rotation lengths were simulated with or without leave areas and with or without understorey vegetation in a subwatershed of Carnation Creek, Vancouver Island, British Columbia. The combined in?nite slope and distributed hydrologic models used to calculate safety factor revealed that most of the simulated landslides were clustered within a 5 to 17 year period after initial harvesting in cases where suf?cient time (c. 50 years) lapsed prior to the next harvesting cycle. Partial cutting produced fewer landslides and reduced landslide volume by 1·4‐ to 1·6‐fold compared to clearcutting. Approximately the same total landslide volume was produced when 100 per cent of the site was initially clearcut compared to harvesting 20 per cent of the area in successive 10 year intervals; a similar ?nding was obtained for partial cutting. Vegetation leave areas were effective in reducing landsliding by 2‐ to 3‐fold. Retaining vigorous understorey vegetation also reduced landslide volume by 3·8‐ to 4·8‐fold. The combined management strategies of partial cutting, increasing rotation length, provision of leave areas, and retention of viable understorey vegetation offer the best alternative for minimizing landslide occurrence in managed forests. Copyright © 2003 John Wiley & Sons, Ltd.     

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

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

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Effects of slope angle and aspect on plant cover and species richness in a humid Mediterranean badland

Soil erosion is one of the most severe land degradation processes in the Mediterranean region. Although badlands occupy a relatively small fraction of the Mediterranean area, their erosion rates are very high. Many studies have investigated to what extent vegetation controls soil erosion rates. This study, however, deals with the impact of erosion on vegetation establishment. In semi‐arid badlands of the Mediterranean, soil water availability constitutes the main limiting factor for vegetation development. As a consequence, south‐facing slopes are typically less vegetated due to a very large water stress. However, these findings do not necessarily apply to humid badlands. The main objective of this paper is to determine the topographic thresholds for plant colonization in relation to slope aspect and to assess the spatial patterns of vegetation cover and species richness. We surveyed 179 plots on highly eroded badland slopes in the Central Pyrenees. We defined four aspect classes subdivided into slope angle classes. Colonization success was expressed in terms of vegetation cover and species richness. Slope angle thresholds for plant colonization were identified for each slope aspect class by means of binary logistic regressions. The results show that a critical slope angle exists below which plants colonize the badland slopes. Below this critical slope angle, plant cover and species richness increase with a decreasing slope angle. The largest critical slope angles in humid badlands are observed on south‐facing slopes, which contrasts with the results obtained in semi‐arid badlands. North‐facing slopes however are characterized by a reduced overall vegetation cover and species richness, and lower topographic threshold values. The possible underlying processes responsible for this slope‐aspect discrepancy in vegetation characteristics are discussed in terms of environmental variables that control regolith development, weathering and erosion processes. Moreover, possible restoration strategies through the use of vegetation in highly degraded environments are highlighted. Copyright © 2014 John Wiley & Sons, Ltd.