2 search hits
-
Soil-Landscape Modelling in an Andean Mountain Forest Region in Southern Ecuador
(2010)
-
Mareike Ließ
- Soil-landscapes are diverse and complex due to the interaction of pedogenetic, geo-morphological and hydrological processes. The resulting soil profile reflects the balance of these processes in its properties. Early conceptual models have by now resulted into quantitative soil-landscape models including soil variation and its unpredictability as a key soil attribute. Soils in the Andean mountain rainforest area of southern Ecuador are influenced by hillslope processes and landslides in particular. The lack of knowledge on the distribution of soils and especially physical soil properties to understand slope failure, resulted in the study of this particular soil-landscape by means of statistical models relating soil to terrain attributes, i.e. predictive soil mapping. A 24 terrain classes comprising sampling design for soil investigation in mountainous areas was developed to obtain a representative dataset for statistical modelling. The soils were investigated by 56 profiles and 315 auger points. The Reference Soil Groups (RSGs) Histosol, Stagnosol, Umbrisol, Cambisol, Leptosol and Regosol were identified according to the World Reference Base for Soil Resources (WRB). While soil profiles and auger points were described in their horizon composition, thickness, soil cohesion, bulk density and texture were analysed in soil profiles only. The prediction of soil parameters was carried out with Classification and regression tree (CART) and Random Forest (RF) method. At this, prediction uncertainty was addressed with hundredfold model runs based on different random Jackknife partitions. Problems with the prediction of the RSGs, likely caused by inconsequence within the WRB (absolute and relative values as decision criteria), resulted in the proposal of “incomplete soil classification”, which relates the thickness of the diagnostic WRB horizons to the upper 100 soil centimetres. Histosol and Stagnosol have been distinguished as dominant RSGs within the inves-tigation area. While Histosol probability depended on hydrological parameters; highest Stagnosol probability was predicted on slopes < 40° and above 2146 m a.s.l. Whether the first mineral soil horizon displays stagnic properties or not, likely depends on physical soil properties in addition to terrain parameters. Incomplete soil classification resulted in histic and stagnic soil parts dominating the first 100 cm of the soil volume for most of the research area. Comparing CART and Random Forest (RF) in their model performance to predict topsoil texture and bulk density as well as mineral soil thickness by hundredfold model runs with random Jackknife partitions, RF predictions resulted more powerful. Altitude a.s.l. was the most important predictor for all three soil parameters. Increasing sand/ clay ratios with increasing altitude, on steep slopes and with overland flow distance to the channel network are caused by shallow subsurface flow removing clay particles downslope. Deeper soil layers are not influenced by the same process and therefore showed different texture properties. Terrain parameters could only explain the spatial distribution of topsoil properties to a limited extent, subsoil properties could not be predicted at all. Other parameters that likely influence soil properties within the investigation area are parent material and landslides. Strong evidence was found that topsoil horizons did not form from the bedrock underlying the soil profile. Parent material changes within short distance and often within one soil profile. Landslides have a strong influence on soil-landscape formation in shifting soil and rock material. Soil mechanical and hydrological properties in addition to terrain steepness were hypothesized to be the major factors in causing soil slides. Thus, the factor of safety (FS) was calculated as the soil shear ratio that is necessary to maintain the critical state equilibrium on a potential sliding surface. The depth of the failure plane was assumed at the lower boundary of the stagnic soil layer or complete soil depth, depending on soils being stagnic or non-stagnic. The FS was determined in dependence of soil wetness referring to 0.001, 0.01, 0.1 and 3 mm/h net rainfall rate. Sites with a FS ≥ 1 at 3 mm/h (complete saturation) were classified as unconditionally stable, sites with a FS < 1 at 0.001 mm/h as unconditionally unstable. The latter coincided quite well with landslide scars from a recent aerial photograph.
-
Water flow paths in soils of an undisturbed and landslide affected mature montane rainforest in South Ecuador
(2010)
-
Folkert Christian Bauer
- The number of previous hydrological studies concerning water flow paths in tropical montane rainforest is small. However, due to the increasing pressure of deforestation and land use change comprehensive knowledge of these natural ecosystems is needed if sustainable land use strategies should keep negative effects of human impacts on water flow paths as low as possible. In this context, present work addresses the identification, characterisation, and modelling of water flow paths in soils of an undisturbed and landslide affected natural Andean forest ecosystem in the south of Ecuador whose deforestation rate is one of the highest in South America. In an investigation area situated in the Andes of South Ecuador, in gentler slopes and altitudes above 2100 m ASL mainly Stagnosols and Histosols with stagnic colour pattern and low to negligible rock fragment content prevail. With increasing altitude the abundance of these soils increase, while the presence of Cambisols and Regosols is most pronounced below 2100 m ASL and clearly correlated with the slope angle. Therefore, these soils were mainly encountered in steeper, particularly landslide affected sites often resembling a melange of fine soil and high contents of rock fragments. Aside the investigation of the influence of the rock fragment content on soil hydrological- and physical parameters such as the relationship between rock fragment content and saturated hydraulic conductivity of the soil, present study aims particularly to investigate flow paths of water in soils of landslide affected and unaffected hillslopes. Therefore, we employed conventional field- and laboratory methods, dye tracer experiments including an appropriate image processing technique, as well as statistical models. Results show that both rock fragment content and bulk density control significantly, but not largely the saturated hydraulic conductivity of the mineral soils. Dye tracer experiments and soil parameters document a deeper percolation in the landslide affected hillslopes than in the landslide unaffected hillslopes, where we found preferential flow in root channels with low soil matrix interaction as dominant flow mechanism. A surface near quasi impervious layer along the interface between topsoil and subsoil limits percolation of the water giving the prerequisites of a lateral shallow subsurface flow along the interface between topsoil and organic layer. This is in line with previous studies performed in the same investigation area which already proved indirectly the existence of this flow. However, in none of these studies the shallow subsurface flow was assigned to certain slope inclinations or altitudes. Due to a recently published digital soil map and the results we obtained from the landslide unaffected sites, we know that particularly in hillslopes of less than 30 degrees above 2100 m ASL prerequisites are given for spatially extended shallow subsurface flow. However, even if these prerequisites are not evident for the landslide affected hillslopes, we cannot exclude the possibility of shallow subsurface flow occurrence here since soil cover of the steep terrain is relatively shallow while rainfall is high throughout the year. Therefore, and given that key parameters such as permeability of subsoil and bedrock, interception and evaporation remain unclear or were investigated exclusively such as the spatial variability of the saturated hydraulic conductivity, we conducted a series of virtual experiments in order to assess the potential occurrence of shallow subsurface flow in Cambisols below 2100 m ASL. In these experiments we also included the organic layer being highly abundant in the investigation area, whose hydraulic parameters were estimated by means of inverse numerical modelling. The virtual experiments were based on a two dimensional finite element model representing a steep forested hillslope transect of approx. 54 m length. Aside soil properties, evapotranspiration and interception, the model included the spatial variability of the saturated hydraulic conductivity, the pressure head and their spatial trends. The results of virtual experiment series show that a sound evidence of the key parameters aforementioned is obligate if process conceptualisation regarding shallow subsurface flow generation, but also landslide initiation, solute and matter transport is in the spotlight.
