Slope Scale Research Articles

VARIABILIDADE ESPACIAL DA GERAÇÃO E CONECTIVIDADE DE SEDIMENTOS EM BACIAS DE UHES

The knowledge of the spatial variability of the generation and transport of sediments allows the identification with considerable precision critical areas regarding the generation and flow of detrital materials in atersheds. This paper aims to relate the spatial variability of soil loss and sediment connectivity in a contribution basin of the Batalha Hydroelectric Power Plant reservoir in the municipality of Cristalina (GO). The methodology comprised the bivariate spatial correlation between production estimates (MUSLE) and the Connectivity Index (CI) applied at a slope scale in the basin's contribution area. The results indicate the configuration of four spatial patterns: low production and low connectivity index in most of the area, comprising the higher and flatter portions; low production and high connectivity index in the vicinity of springs and drainage channels; high production and low connectivity index in the higher, steeper and more distant portions; and high production and high connectivity, predominating in the steepest portions and closest to the channels or even in those portions that are more distant but are connected to the channel, in the latter situation via more intense surface runoff flow lines. In this sense, the importance of this methodology is highlighted in the identification of critical areas, which are priority areas for the implementation of mitigating measures of the impacts resulting from hydric erosive processes. Keywords: Water Erosion, Sediment Contribution, Critical Areas, Water Reservoirs.

Snow Level From Post‐Processing of Atmospheric Model Improves Snowfall Estimate and Snowpack Prediction in Mountains

AbstractIn mountains, the precipitation phase greatly varies in space and time and affects the evolution of the snow cover. Snowpack models usually rely on precipitation‐phase partitioning methods (PPMs) that use near‐surface variables. These PPMs ignore conditions above the surface thus limiting their ability to predict the precipitation phase at the surface. In this study, the impact on snowpack simulations of atmospheric‐based PPMs, incorporating upper atmospheric information, is tested using the snowpack scheme Crocus. Crocus is run at 2.5‐km grid spacing over the mountains of southwestern Canada and northwestern United States and is driven by meteorological fields from an atmospheric model at the same resolution. Two atmospheric‐based PPMs were considered from the atmospheric model: the output from a detailed microphysics scheme and a post‐processing algorithm determining the snow level and the associated precipitation phase. Two ground‐based PPMs were also included as lower and upper benchmarks: a single air temperature threshold at 0°C and a PPM using wet‐bulb temperature. Compared to the upper benchmark, the snow‐level based PPM improved the estimation of snowfall occurrence by 5% and the simulation of snow water equivalent (SWE) by 9% during the snow melting season. In contrast, due to missing processes, the microphysics scheme decreased performances in phase estimate and SWE simulations compared to the upper benchmark. These results highlight the need for detailed evaluation of the precipitation phase from atmospheric models and the benefit for mountain snow hydrology of the post‐processed snow level. The limitations to drive snowpack models at slope scale are also discussed.

Assessment of probability of failure on rainfall-induced shallow landslides at slope scale using a physical-based model and fuzzy point estimate method

Shallow hillslope failure caused by rainfall is characterized by complex soil hydrology and mechanical behavior. It is important to understand the hydraulic behavior of hillslopes and quantify the effect of the uncertainty of mechanical parameters on hillslope stability for forewarning and hillslope management. Intra-hole deformation and displacement were record for the hillslope of the Babaoliao collapse site in the Chiayi County, as a case study. The fuzzy point estimation method and physical-based model were combined with the local factor of safety (LFS) theory to calculate the internal local factor of safety of the hillslope. A reliability analysis was then performed to determine the failure probability at different depths. Historical rainfall events were used to validate the model and predict the development of the failure probability for different rainfall patterns with the same warning rainfall. The results revealed that the failure probability model could effectively predict the area of hillslope instability and its changes over time and space. Different rainfall patterns affected the infiltration flux, leading to the difference in hillslope failure time. The delayed rainfall pattern had a significant impact on the time of slope instability, and shallow collapse was most likely to occur earlier. This study can be used as a reference for developing future hillslope warnings.

Effects of material migration on the spatial distribution of topsoil moisture at the slope scale

Soil moisture is strongly spatially variable and scale dependent, but there is a lack of research focusing on factors that link soil actual evaporation to the spatial distribution of topsoil moisture at the slope scale. External factors (e.g., meteorological) have similar effects across entire slopes; hence, this paper instead explores the internal factors controlling differential soil actual evaporation. We consider the effects of historical rainfall on soil erosion, and adopt a combined approach encompassing slope investigations, a slope rainfall simulation experiment, slope erosion analyzed by 3D laser terrain scanning, and evaporation experiments. Along the runoff direction, we show that the relative topsoil clay content increases, but slope surface deformation and failure lead to a decrease in the relative content of clay in deposits located on the lower of slope. The original thickness of slope surface decreases in the upper part of the slope, due to particle migration, while the slope surface relatively far away from the upper part of the slope tends to be thicker, due to the accumulation of slope deposits. The actual evaporation of the topsoil relatively far from the upper part of the slope is significantly weaker (63%–72% probability), which may be related to the high clay content, dense soil structure, and high content of clay minerals. The topsoil evaporation also tends to weaken as overburden thickness increases. Here, we suggest that marked differences in the original thickness of the slope surface and in the properties of the topsoil caused by the migration and deposition of soil particles under successive years of rainfall erosion, are important factors controlling differences in topsoil actual evaporation in different parts of the slope. This understanding lays the foundation for the analysis of hydrological processes in different parts of a slope, which is crucial for the study of shallow slope mechanical behavior.

The Spatiotemporal Characteristics of Flow-Sediment Relationships in a Hilly Watershed of the Chinese Loess Plateau.

The flow–sediment relationship is important to understand soil erosion and sediment transport in severely eroded areas, such as Loess Plateau. Previous research focused on the variation and driving forces of runoff and sediment at the different scales in a watershed. However, the variations of the flow–sediment relationship on multispatial scales (slope, subgully, gully, and watershed scales) and multitemporal scales (annual, flood events, and flood process) were less focused. Taking the Peijiamao watershed, which includes whole slope runoff plot (0.25 ha, slope scale), branch gully (6.9 ha, subgully scale), gully (45 ha, gully scale), and watershed (3930 ha, watershed scale), four different geomorphic units located at the Chinese Loess Plateau, as the research site, a total of 31 flood events from 1986 to 2008 were investigated, and two flood process data were recorded across all the four geomorphic units. The results showed that on the annual timescale, the average sediment transport modulus and runoff depth at four scales exhibited a linear relationship, with determination coefficients of 0.81, 0.72, 0.74, and 0.77, respectively. At the flood event timescale, the relationships between sediment transport modulus and runoff depth at the gully and watershed scales could also be fitted with a linear relationship with high determination coefficients (from 0.77 to 0.99), but the determination coefficient at the slope scale was only 0.37 at the event scale. On the single rainfall event timescale, the flow–sediment relationship at the slope scale showed a figure-eight hysteretic pattern while those relationships at larger scales showed an anticlockwise loop hysteretic pattern. Under the same flow condition, the suspended sediment concentrations during the falling stage were significantly higher than those during the rising stage. Moreover, the difference was bigger as the spatial scale increased due to the wash loads in the downstream gullies, which favored the occurrence of hyper-concentration flow. The results of the study could provide useful insights into the temporal–spatial scale effects of sediment transport and their internal driving mechanisms at the watershed scale.

Relationship between community composition and water infiltration of biological soil crusts in the hilly Loess Plateau, China

Biological soil crusts (biocrusts) are the mixed community composed of different ratios of cyanobacteria, mosses, and lichens at the slope scale in the Hilly Loess Plateau region. Biocrusts significantly affect water infiltration in this area. The relationship between infiltration rate and community structure of mixed biocrusts is unknown, which would hinder the assessment of soil permeability of biocrusts at the slope scale. We measured the stable infiltration rate of cyanobacteria, moss, and mixed biocrusts with different proportions of cyanobacteria and moss including moss coverage of <15%, 15%-30%, 30%-45%, 45%-60% and >60%, respectively. The principal component analysis and path analysis were used to understand the influencing factors of stable infiltration rate of mixed biocrusts, and to clarify the relationship between the stable infiltration rate and the community structure of mixed biocrusts. The results showed that the saturated hydraulic conductivity of cyanobacteria and moss crusts was 0.66 mm·min-1 and 2.40 mm·min-1, respectively. The stable infiltration rates of mixed biocrusts with moss coverage <15% to >60% were 0.80-2.30 mm·min-1. The stable infiltration rate of mixed biocrusts at the slope scale depended on moss coverage and its improvement on soil pore structure, with the correlation coefficients being 0.636 (P=0.011) and 0.835 (P=0.000) respectively. Herein, the saturated hydraulic conductivity and coverage of cyanobacteria and moss confirmed the weighted prediction of water infiltration volume (y) i.e., a significant correlation (r=0.945) with the measured water infiltration volume (x) of mixed biocrusts. The linear fitting of measured and predicted water infiltration volume of mixed biocrusts was y=0.85x (R2=0.98, P<0.05). This study clarified the relationship between water infiltration of mixed biocrust community composition and individual biocrust composition, which provided a scientific basis for accurately evaluating the hydrological process of biocrusts in this area.

Route Stability Index (RSI): an Index for the Assessment of Rockfall-Related Hazards in Rock Slopes Equipped for Sport Climbing

A multidisciplinary geomechanical method to evaluate and quantify rockfall hazards in rock slopes equipped for sport climbing is presented. This method exploits close-range survey data obtained from terrestrial laser scanner (TLS) and unmanned aerial vehicle-based digital photogrammetry (UAV-DP) and integrates them with traditional scanline survey data collected in the field. Based on the acquired information, a kinematic analysis is performed to predict the location of structurally controlled rock failure mechanisms. An index denoted as route stability index (RSI) is then proposed to quantify the associated hazard along sport climbing routes in the study area. In addition, hazard is also assessed at the overall slope scale. The method was successfully applied at an abandoned climbing site historically affected by rockfalls, 5 km north of Florence (Italy). RSI can be adopted as part of the planning and monitoring strategies for the management of risk related to outdoor climbing activities, which are becoming an important part of Italian and European geo-tourism.

The diversity of spiral galaxies explained

ABSTRACT An extensive catalogue of spatially resolved galaxy rotation curves (RCs) and multiband optical light profiles for 1752 observed spiral galaxies is assembled to explore the drivers of diversity in galaxy structural parameters, RC shapes, and stellar mass profiles. Similar data were extracted from the ‘Numerical Investigation of a Hundred Astrophysical Objects’ galaxy simulations to identify any differences between observations and simulations. Several parameters, including the inner slope $\mathcal {S}$ of a RC, were tested for diversity. Two distinct populations are found in observed and simulated galaxies: (i) blue, low-mass spirals with stellar mass M⋆ ≲ 109.3 M⊙ and roughly constant $\mathcal {S}$; and (ii) redder, more massive and more diverse spirals with rapidly increasing $\mathcal {S}$. In all cases, the value of $\mathcal {S}$ seems equally contributed by the baryonic and non-baryonic (dark) matter. Diversity is shown to increase mildly with mass. Numerical simulations reproduce well most baryon-dominated galaxy parameter distributions, such as the inner stellar mass profile slope and baryonic scaling relations, but they struggle to match the full diversity of observed galaxy RCs (through $\mathcal {S}$) and most dark matter-dominated parameters. To reproduce observations, the error broadening of the simulation’s intrinsic spread of RC metrics would have to be tripled. The differences in various projections of observed and simulated scaling relations may reflect limitations of current subgrid physics models to fully capture the complex nature of galaxies. For instance, active galactic nuclei (AGNs) are shown to have a significant effect on the shapes of simulated RCs. The inclusion of AGN feedback brings simulated and observed inner RC shapes into closer agreement.

The Combined Use of 137Cs Measurements and Zr-Methods for Estimating Soil Erosion and Weathering in Karst Areas of Southwestern China

Soil physical erosion and chemical weathering quantification at the slope scale are important to reveal the material cycle of the ecosystem in the karst region, because of the high heterogeneity due to the complex habitats. The Zr-based mass balance method has been widely used to quantify long-term physical erosion and chemical weathering at the slope scale, but the method is still in the exploratory research stage for quantifying short-term physical erosion and chemical weathering. In this paper, sloping fields (nine sampling sites and two sloping fields) in Zunyi, within the Guizhou karst region, were studied. We attempted to quantify the short-term physical erosion and chemical weathering rates by 137Cs combined with the Zr-based mass balance method, and an ideal distribution curve of the Zr concentration in the soil surface layer of a karst region is proposed. The results showed the following: (1) The average soil erosion rate on the slope of the study area is 580 t/(km2·a), which is equivalent to 14% of the average value of the Wujiang River basin in which it is located. This shows that the spatial distribution of soil erosion in this area varies significantly. (2) The Zr concentration in the erosion profile (EUC (1)) corresponds to 48% of that in the deposition profile (DUC (3)). This indicates that physical erosion in the study area survives chemical weathering, which is also consistent with the relationship model hypothesis. In addition, the vertical distribution characteristics of Zr concentration in all profiles are basically consistent with the ideal hypothesis curve. (3) The chemical weathering rate of the topsoil has been preliminarily estimated to be around 30 t/(km2·a), and the ratio between the physical soil erosion and chemical weathering is 20:1. The results show that the physical erosion is dominant over the chemical weathering of topsoil, and the chemical weathering rate was proportional to the physical erosion. This study provides a new method for quantifying short-term soil erosion and weathering erosion at the slope scale in karst regions, which is important for regional ecological restoration and sustainable development.

A Downscaling Intercomparison Study: The Representation of Slope- and Ridge-Scale Processes in Models of Different Complexity

Spatially distributed meteorological information at the slope scale is relevant for many processes in complex terrain, yet information at this sub-km spatial resolution is difficult to obtain. While downscaling to kilometer resolutions is well described in literature, moving beyond the kilometer scale is not. In this work, we present a methodical comparison of three downscaling methods of varying complexity, that are used to downscale data from the Numerical Weather Prediction model COSMO-1 at 1.1 km horizontal resolution to 250 and 50 m over a domain of highly complex terrain in the Swiss Alps. We compare WRF, a dynamical atmospheric model; ICAR, a model of intermediate complexity; and TopoSCALE, an efficient topography-based downscaling scheme. Point-scale comparisons show similar results amongst all three models w.r.t. mean-error statistics, but underlying dynamics are different. Ridge-flow interactions show reasonable agreement between WRF and ICAR at 250 m model resolution. However, at 50 m resolution WRF is able to simulate complex flow patterns that ICAR cannot. Validation against Lidar data suggests that only WRF is able to capture preferential deposition of snow. Based on these findings and the significant reduction in computational costs, ICAR is a cost efficient alternative to WRF at the 250 m resolution. TopoScale performs very well in point-scale comparisons, but it is unclear if this can be attributed to the model itself or to the forcing data and the observations assimilated therein. Further study is required to quantify this effect.

A Lagrangian heuristic for the multicommodity capacitated location problem with balancing requirements

We consider the multicommodity capacitated location problem with balancing requirements. This problem has applications in transportation and was introduced in the context of planning the movements of containers between ports and terminals, identified as depots in the problem formulation. We propose a Lagrangian relaxation based on a reformulation of the problem that allows the subproblem to be decomposed by depot, giving rise to a series of continuous knapsack subproblems. We optimize the Lagrangian dual by a bundle method and we develop a Lagrangian heuristic based on slope scaling. Numerical results are presented on a set of large-scale benchmark instances. We obtain good quality solutions within short computing times using the Lagrangian heuristic approach.

Estimation of Wave-Breaking Index by Learning Nonlinear Relation Using Multilayer Neural Network

Estimating wave-breaking indexes such as wave height and water depth is essential to understanding the location and scale of the breaking wave. Therefore, numerous wave-flume laboratory experiments have been conducted to develop empirical wave-breaking formulas. However, the nonlinearity between the parameters has not been fully incorporated into the empirical equations. Thus, this study proposes a multilayer neural network utilizing the nonlinear activation function and backpropagation to extract nonlinear relationships. Existing laboratory experiment data for the monochromatic regular wave are used to train the proposed network. Specifically, the bottom slope, deep-water wave height and wave period are plugged in as the input values that simultaneously estimate the breaking-wave height and wave-breaking location. Typical empirical equations employ deep-water wave height and length as input variables to predict the breaking-wave height and water depth. A newly proposed model directly utilizes breaking-wave height and water depth without nondimensionalization. Thus, the applicability can be significantly improved. The estimated wave-breaking index is statistically verified using the bias, root-mean-square errors, and Pearson correlation coefficient. The performance of the proposed model is better than existing breaking-wave-index formulas as well as having robust applicability to laboratory experiment conditions, such as wave condition, bottom slope, and experimental scale.

EQUIVALENT GEOLOGICAL STRENGTH INDEX (GSI) APPROACH WITH APPLICATION TO ROCK MASS SLOPE STABILITY

A considerable amount of slope stability analysis has been observed in jointed rock masses in which the GSI (Geological Strength Index) estimated at the outcropping level is considered input data to define the rock mass strength. However, this procedure is unsuitable when the rock outcrop scale and the slope scale are significantly different (e.g. open-pit slopes), resulting in an overestimated rock mass strength. For this reason, and in the absence of criteria to modify the GSI based on the scale effects, in this research, a new GSI version is proposed, called GSIe or “equivalent GSI”. To define an expression for obtaining the GSIe in terms of the rock mass properties, comparative stability analyses were conducted in a series of hypothetical slopes using two approaches: the first considers the rock mass as a discontinuous medium of rock blocks separated by discontinuities; the second considers the rock mass as an equivalent continuous medium characterized by an equivalent GSI. For the adequate equivalent GSI value, evaluated in each analyzed slope, the safety factor and the failure surface are similar in both approaches. In conformity with the results, a GSIe formulation in terms of the slope height, the spacing, the intact rock strength, the persistence, and the joint conditions has been proposed. Finally, the formulation was validated by applying it in five cases of mining slopes where the failure occurred.

Mitigation of soil organic carbon mineralization by soil redistribution - An erosion-deposition plot study under natural rainfall over five years

As the essential driving force of soil redistribution, runoff and sediment are not often considered when quantifying and integrating the effects of soil erosion and deposition on soil CO2 emissions. Therefore, in this study, variations in soil CO2 emissions from erosion-deposition plots were regularly monitored in China’s Loess Plateau (2015–2019). The cumulative soil CO2 emissions from the depositional zones with slope gradients of 5°, 10°, and 20° were increased by 0.4–16.7%, 20.1–32.6%, and 31.9–51.5%, respectively, than those of their respective eroding slopes. Relative to the 5° eroding slopes, the cumulative soil CO2 emissions decreased by 2.8–13.5% and 11.3–15.6% on the steeper 10° and 20° eroding slopes, respectively. Conversely, cumulative soil CO2 emissions increased by 1.0–7.9% and 6.9–13.3% in the depositional zones of 10° and 20°, respectively, as compared with that in the 5° depositional zones. Considering both the eroding slopes and depositional zones, the total amount of CO2 emissions from the 10° and 20° erosion-deposition plots were 114 g CO2-C y−1 and 177 g CO2-C y−1 lower than those from the 5° erosion-deposition plots, respectively. This can be attributed to the combined effects of the amounts of runoff and sediment displacement on soil CO2 emissions from eroding slopes and depositional zones. Our results indicate that, soil erosion and deposition together have great potential to mitigate soil organic carbon mineralization and preserve soil carbon pools in slope lands. Furthermore, soil organic carbon mineralization at the slope scale is sensitive to displaced runoff and sediment. Therefore, possible variations in soil-atmospheric carbon exchanges induced by erosion-displaced runoff and sediment should be properly accounted for when attributing carbon fluxes in regions dominated by sloping landscapes.

Cultural Heritage and Rockfalls: Analysis of Multi-Scale Processes Nearby the Lucus Angitiae Archaeological Site (Central Italy)

Archaeological areas in the mountain region of central Italy can be seriously threatened by geological hazards, and efforts are required to preserve cultural heritage. The Lucus Angitiae is a pre-Roman site located along the western edge of the Fucino Basin, the largest continental depression of central Apennines. The carbonate slope overhanging the area is affected by active rockfall processes from two main rock escarpments. In this paper, rockfall assessment was pursued through a 3D kinematic modelling, performed by adopting a probabilistic approach. Specific attention was dedicated to the choice and calibration of the input data, based on field evidence and a literature review. Two different sizes of wedge-shaped rock blocks were identified on rock escarpments, and specific stability analyses were performed. Sensitivity analyses accounting for possible triggering factors, such as water pressure increase and seismic action, were also carried out, together with an investigation of the seismological characteristics of the area. The results of the numerical simulations were used to design effective countermeasures in the framework of a mitigation plan for protection of the archaeological site. Finally, clues of gravity-driven slope deformations at the slope scale were documented, framing the rockfall process in a wider geological scenario.

Threshold Definition for Monitoring Gapa Landslide under Large Variations in Reservoir Level Using GNSS

The triggering threshold is one of the most important parameters for landslide early warning systems (EWSs) at the slope scale. In the present work, a velocity threshold is recommended for an early warning system of the Gapa landslide in Southwest China, which was reactivated by the impoundment of a large reservoir behind Jinping’s first dam. Based on GNSS monitoring data over the last five years, the velocity threshold is defined by a novel method, which is implemented by the forward and reverse double moving average of time series. As the landslide deformation is strongly related to the fluctuations in reservoir water levels, a crucial water level is also defined to reduce false warnings from the velocity threshold alone. In recognition of the importance of geological evolution, the evolution process of the Gapa landslide from topping to sliding is described in this study to help to understand its behavior and predict its potential trends. Moreover, based on the improved Saito’s three-stage deformation model, the warning level is set as “attention level”, because the current deformation stage of the landslide is considered to be between the initial and constant stages. At present, the early warning system mainly consists of six surface displacement monitoring sites and one water level observation site. If the daily recorded velocity in each monitoring site exceeds 4 mm/d and, meanwhile, the water level is below 1820 m above sea level (asl), a warning of likely landslide deformation accelerations will be released by relevant monitoring sites. The thresholds are always discretely exceeded on about 3% of annual monitoring days, and they are most frequently exceeded in June (especially in mid-June). The thresholds provide an efficient and effective way for judging accelerations of this landslide and are verified by the current application. The work presented provides critical insights into the development of early warning systems for reservoir-induced large-scale landslides.

Determining the area corrections affecting the map areas in GIS applications

Abstract Nowadays, there are many area-based Geographic Information Systems (GIS) applications such as real estate valuation, land tax, farming support and cost–benefit analysis. Areas used in such applications are calculated by means of two-dimensional plane geometry. However, the computed area value is not the exact area value in the terrain. In order to calculate the exact area value of a parcel, area corrections due to various factors must be taken into account. These factors are selection of projection, slope of the terrain, elevation of the terrain and scale of the map. Selection of projection and slope of terrain are available; elevation of the terrain and scale of map are not available in all GIS software. In this study, the effect of area corrections on the area value calculated from the map is examined with sample applications and the results are presented to the GIS users. According to the results, GIS users should select the equal area projection. In addition, scale of map, elevation and slope of terrain should be taken into account in the area calculation where land measurements are not possible.

Predictive Capabilities of 2D and 3D Block Propagation Models Integrating Block Shape Assessed from Field Experiments

Block propagation models have been used for years for rockfall hazard assessment. However, the calibration of model parameters that allow the simulations to accurately predict rockfall trajectories for a given study site remains a key issue. This research aims at investigating the predictive capabilities of block propagation models after a preliminary calibration phase. It is focused on models integrating the shape of blocks since, despite their sound physical bases, they remain less used than lumped-mass approaches due to their more recent popularisation. Benefiting from both a recently built model integrating block shape, usable in 2D and 3D, and from recent experimental results at the slope scale, we first performed a calibration based on the use of the 2D model, and then we evaluated the predictive capabilities of the calibrated model in 2D and in 3D using the remaining part of the experimental results. The calibrated model simulations predict the main characteristics of the propagation, that is the preferential deposit zones and the ranges of velocities at specific locations. Good matches between simulations and experimental results in both the calibration and validation phases emphasizes the wide applicability of the model: after a calibration phase on a sufficient number of different soil types, the model may be used in a predictive manner. The good match between 2D and 3D simulations also highlights the ease-of-use of the model for field applications, as the 2D model produces sufficiently accurate results while also being easier and faster to calibrate. As classically observed for block propagation models, the model is not sufficient to predict the details of the velocity and stopping points but provides accurate predictions of the global ranges of these quantities, in particular of the extreme values. To lift these limitations in terms of predictive capabilities, more advanced calibration procedures based on optimization techniques constitute a promising path forward.

The Geo-Hydro-Mechanical Properties of a Turbiditic Formation as Internal Factors of Slope Failure Processes

Similar to many inner areas of Southern Europe, the Daunia Apennines are affected by widespread landsliding, often consisting of slow, deep-seated movements. Recurrent acceleration of these landslides causes damage to buildings and infrastructures, severely biasing the socio-economic development of the region. Most landslides in the area of study occur within clayey units of turbiditic flysch formations, often severely disturbed by tectonic thrust and previous landsliding. The Faeto Flysch (FAE) is one of the most widespread turbiditic formations in the Daunia Apennines and is representative of the tectonised geological formations involved in slope failure. This work, by examining the landslide processes occurring at four pilot sites, aims at connecting the observed mechanisms to the geo-hydro-mechanical setup of FAE in the slopes. It is found that the soil portion of FAE consists of highly plastic clays, resulting in low intrinsic shear strength, and hence controls the initiation and progression of failure in the slopes, as such representing an internal predisposing factor to landsliding. In addition, the presence of fractured rock strata confers a high permeability at the slope scale, with respect to that of the soil matrix. This results in severe piezometric levels in the slope, which represent another internal predisposing factor to failure, and in the ability to induce significant seasonal pore water pressure oscillations down to great depths, connected to rainfall infiltration, thus triggering the recurrent acceleration of the landslides.

Determination of the representative elementary area (REA) of biocrusts: A case study from the Hilly Loess Plateau region, China

Biological soil crusts (biocrusts) are mixed communities of cyanobacteria, lichens and mosses in different ratios, contributing to important ecological functions in arid and semiarid regions worldwide. Biocrusts are spatially variable, and the variability in biocrust composition and coverage is scale-dependent. The following question can then be asked: What is the appropriate spatial scale for observing ecological functions? Without clarifying this issue, we cannot fully understand the ecological functions of biocrusts. The key to answering this question is to determine a threshold area, or representative elementary area (REA). Accordingly, we analyzed red–green–blue (RGB) images of 90 biocrust plots (2.0 m × 2.0 m) from nine revegetated grasslands in the Hilly Loess Plateau region of China. The variability in biocrust composition and coverage across the plot sizes was studied by gradually expanding the plot size from 0.01 m2 to 4.00 m2. The results showed that as the plot sizes increased from 0.01 m2 to 4.00 m2, the number of biocrust types logarithmically increased. Biocrust patches of a particular type (such as moss, cyanobacteria or lichen) were often characteristic of a plot size of 0.01 m2, whereas plot sizes larger than 0.25 m2 supported mixed biocrusts of multiple patch types. The variability in coverage of mixed biocrusts logarithmically decreased with increasing plot sizes. The coverage of mixed biocrusts maintained an approximately constant value after a certain critical plot size was reached (1.00 m2 in this study). Our data indicated that REAs of mixed biocrusts exist at the slope scale. The REAs of mixed biocrusts were 0.5–1.0 m2, 1.5–2.5 m2 and 3.0–3.8 m2, with alpha (α) values of 0.1, 0.05 and 0.01, respectively. The size of the REAs on the north-facing slope was larger than that on the south-facing slope, and the patch density of biocrusts had an important influence on the REAs. The results of this study could provide a method for determining the REAs of mixed biocrusts and guide surveys and experimental layouts.