Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Use of stable water isotopes to assess sources and influences of slope groundwater on slope failure

This study employs stable oxygen and hydrogen isotopes as natural tracers to assess the headwater of a landslide next to a drainage divide and the importance of the slope's headwater in the study area. The study is undertaken near Wu‐She Township in the mountains of central Taiwan. Because a reservoir is located on the other side of the divide, this study evaluates the relationship between the reservoir water and headwater of the landslide as well. Over a 1‐year period, water samples from September 2008 to September 2009, including local precipitation (LP), Wu‐She Reservoir's water (WSRW), slope groundwater (SGW), upper‐reach stream water (USTW), and down‐reach stream water (DSTW), were analysed for deuterium (δD) and oxygen (δ¹⁸O) stable isotopes. Results indicate that WSRW is the predominant component in SGW: approximately 70% of SGW originates from WSRW and 30% from LP based on a two end‐member mass‐balance mixing model for δ¹⁸O. The similar two end‐member mixing model is also employed to assess the contributions of USTW and SGW to DSTW. Model results indicate that SGW is the major source of DSTW with a contribution of about 67%. Accordingly, about 47% of DSTW sources from the WSRW. In short, owing to reservoir leakage, WSRW contributes the greater part of both SGW and DSTW. Plentiful WSRW in SGW threatens the stability of the slope in the divide area. To avoid subsequent continuous slope failure, necessary mitigation steps are required. Copyright © 2011 John Wiley & Sons, Ltd.     

Internal soil moisture and piezometric responses to rainfall-induced shallow slope failures

Better knowledge regarding internal soil moisture and piezometric responses in the process of rainfall-induced shallow slope failures is the key to an effective prediction of the landslide and/or debris flow initiation. To this end, internal soil moisture and piezometric response of 0.7-m-deep, 1.5-m-wide, 1.7-m-high, and 3.94-m-long semi-infinite sandy slopes rested on a bi-linear impermeable bedrock were explored using a chute test facility with artificial rainfall applications. The internal response time defined by the inflection point of the soil moisture and piezometric response curves obtained along the soil–bedrock interface were closely related to some critical failure states, such as the slope toe failure and extensive slope failures. It was also found that the response times obtained at the point of abrupt bedrock slope decrease can be used as indicators for the initiation of rainfall-induced shallow slope failures. An investigation of spatial distributions of soil water content, ω (or degrees of saturation, S_r), in the slope at critical failure states shows that the 0.2 m – below – surface zone remains unsaturated with S_r 40–60%, regardless of their distances from the toe and the rainfall intensity. Non-uniform distributions of ω (or S_r) along the soil–bedrock interface at critical failure states were always associated with near-saturation states (S_r 80–100%) around the point of bedrock slope change or around the transient ‘toe’ upstream of the slumped mass induced by the retrogressive failure of the slope. These observations suggest the important role of the interflow along the soil–bedrock interface and the high soil water content (or high porewater pressure) around the point of bedrock slope deflection in the rainfall-induced failure of sandy slopes consisting of shallow impermeable bedrocks. The present study proposes an ‘internal response time’ criterion to substantiate the prediction of rainfall-induced shallow slope failures. It is believed that the ‘internal response time’ reflects the overall characteristics of a slope under rainfall infiltration and can be as useful as the conventional meteorology-based threshold times. The ‘internal response time’ theory can be generalized via numerical modeling of slope hydrology, slope geology and slope stability in the future.     

Finite element analysis of steep excavation slope failure by CFS theory

The distribution of Coulomb failure stress (CFS) change in the steep excavation slope is calculated by finite element method in this paper, and the failure mechanics under different conditions have been investigated. Comparing the CFSs before and after the slope excavation (stress loading and unloading processes), the dangerous internal zone and the most likely failure external area are attained. Given the shear cracks on the top surface while tensile stress or cracks along the toe of the slope, we analyze the high cutting-angle steep slope in Kaixian county of the Three Gorges Reservoir region. We bring forward that the peak value of CFS after excavation can reach to the order of 0.1 MPa, which is greatly higher than that of before. Our preliminary results are useful for optimizing the reinforcement structure during the steep slope stabilization engineering.     

Rock slope failure in a recently deglaciated permafrost rock wall at Piz Kesch (Eastern Swiss Alps), February 2014

In February 2014, a rock pillar with a volume of around 150 000  m³ collapsed at Piz Kesch in the Eastern Swiss Alps. A reconstruction of the conditions prior to the event and of the event itself is presented on the basis of different sources of data. The methods applied include photogrammetry, terrestrial laser scanning, structural geological analysis, examination of meteorological data, carbon‐14 (¹⁴C) dating of organic material in permafrost ice from a tension crack and numerical modelling of likely modes of failure. Despite a complete lack of in situ measurements in the rock wall prior to the event and of direct observations during the event, the available data allow the determination of the approximate timing of the event as well as the structural predisposition, the probable mode of failure and the timescale of several millennia involved in the triggering of the failure of the rock pillar. The interdisciplinary analysis of this event contributes towards understanding the complex interaction of processes involved in large rock slope failures currently occurring in warming mountain permafrost regions. Copyright © 2016 John Wiley & Sons, Ltd.     

Field studies on the influence of rainfall intensity,vegetation cover and slope length on soil moisture infiltration on typical watersheds of the Loess Plateau,China

Soil moisture is a key process in the hydrological cycle. During ecological restoration of the Loess Plateau, soil moisture status has undergone important changes, and infiltration of soil moisture during precipitation events is a key link affecting water distribution. Our study aims to quantify the effects of vegetation cover, rainfall intensity and slope length on total infiltration and the spatial variation of water flow. Infiltration data from the upper, middle and lower slopes of a bare slope, a natural grassland and an artificial shrub grassland were obtained using a simulated rainfall experiment. The angle of the study slope was 15° and rainfall intensity was set at 60, 90, 120, 150, and 180 mm/hr. The effect these factors have on soil moisture infiltration was quantified using main effect analysis. Our results indicate that the average infiltration depth (ID) of a bare slope, a grassland slope and an artificial shrub grassland slope was 46.7–73.3, 60–80, and 60–93.3 cm, respectively, and average soil moisture storage increment was 3.5–5.7, 5.0–9.4, and 5.7–10.2 mm under different rainfall intensities, respectively. Heavy rainfall intensity and vegetation cover reduced the difference of soil infiltration in the 0–40 cm soil layer, and rainfall intensity increased surface infiltration differences on the bare slope, the grassland slope and the artificial shrub grassland slope. Infiltration was dominated by rainfall intensity, accounting for 63.03–88.92%. As rainfall continued, the contribution of rainfall intensity to infiltration gradually decreased, and the contribution of vegetation cover and slope length to infiltration increased. The interactive contribution was: rainfall intensity * vegetation cover > vegetation cover * slope length > rainfall * slope length. In the grass and shrub grass slopes, lateral flow was found at a depth of 23–37 cm when the slope length was 5–10 m, this being related to the difference in soil infiltration capacity between different soil layers formed by the spatial cross-connection of roots.     

Multivariate analysis of soil moisture history for a hillslope

In this study, the spatial distribution of measured soil moisture was analyzed on the platform of multivariate modeling. Soil moisture time series for two seasons were selected and used for analysis to reveal similarities and differences in soil moisture responses for a few rainfall events. The development of a soil moisture transport process that considers the representative element volume and uncertainty of soil media provides the hydrological basis for time series modeling. The systematic procedure of Box–Jenkins with noise modeling was used to delineate the final models for all monitoring points. The physical basis of mass balance and the continuity in inflow contribution, as well as statistical criteria, were used in the model selection procedure. Heuristic approaches provide the spatial distribution of selected models along the transect of a hillside. Comparative analysis for two different depths and seasons provide an understanding of the variation in soil moisture transfer processes at the hillslope scale. Differences in soil moisture models for both depths and seasons are associated with eco-hydrological processes. The relationships between distributed topographic features and modeling results were explored to configure dominant hydrological processes for each season.     

Stability analysis of probabilistic soil moisture dynamics

We study the dynamics of soil moisture in the presence of precipitation, evapotranspiration and leakage. This study utilizes a probabilistic framework, and considers the evolution of the probability density function of soil moisture. Our main objective is to assess when the stationary solution is a physically realizable state of the system. Mathematically, this is equivalent to analyzing when the stationary solutions are stable. We show the following main results: For a general rainfall and evapotranspiration model, in a multiple plant framework, the stationary solution is a stable attractor in L₁. However, the stationary solution is at best conditionally stable in L_∞. This has important implications in terms of (1) the validity of moment analysis of the transient problem, (2) the physical importance of transient solutions, and (3) the numerical solution of the stationary problem as a late time solution of the transient problem.     

Responses of soil moisture at different topographic positions to rainfall events along a steep subtropical forested hillslope

Understanding the dynamic response of soil moisture to rainfall is crucial for describing hydrological processes at the hillslope scale. However, because of sparse monitoring coupled with the complexity of water movement and steep topography, the findings of rainfall-related soil moisture dynamics have not always been consistent, indicating a need for systematic investigations of soil moisture dynamics and infiltration patterns following rainfall inputs at multiple topographic positions along a hillslope. This study aimed to examine the nature of these responses by characterizing and quantifying the response amplitude, rate and time for 37 large rainfall events at 25 combinations of topographic positions and soil depths along a steep forested hillslope. Our results showed that soil moisture responses under different rainfall patterns could be attributed to one or the other rainfall characteristics, such as rainfall intensity and amount. However, soil moisture dynamics at different hillslope positions after rainfall varied widely due to the controls of soil properties, topography, and non-equilibrium flow. Preferential flow was more evident under dry initial soil conditions than under wet initial soil conditions. Findings of this study reveal that the dynamic response patterns of soil moisture to rainfall do not always follow topographic controls, which can improve our understanding of water cycling related to the infiltration process at the hillslope scale, and support water resources management in subtropical mountain ecosystems.     

Initiation of earthquake-induced slope failure: influence of topographical and other site specific amplification effects

Increased structural damage caused byearthquakes on hilltops and along ridgeshas often been related to amplification ofground motion due to the presence oftopography. However, comparison betweenobservations and numerical modeling hasshown that amplification is only partlydependent on the prominent surfacemorphology. Strong effects are also inducedby soft layers, such as weathered rockmaterial or colluvium, covering thetopographies.Numerous seismically triggered landslidesare reported to occur in the same materialsthat are likely to amplify ground motions.Therefore, it can be suspected that groundmotion dynamics significantly contribute tothe observed slopes failures. Thispotential relationship is the subject ofthe present case study, the Ananevorockslide in the northeastern Tien Shanmountains. The survey included geophysicalprospecting, earthquake recordings andstructural analyses of the rock fabric. Onthe basis of the field data, observedamplification effects could be related tothe local geological conditions andparticularly to the surface morphology andto the presence of low-velocity layers –deeply weathered rocks – on the top of thebedrock. Surface layer- andtopography-dependent amplification has alsobeen studied numerically by 2D and 3Dfinite element modeling of ground motiondynamics. The present paper focuses onadditional effects that may be induced bythe presence of a fault zone and of thelandslide scarp. Further, observed andcomputed ground motion dynamics areconnected with slope failuresusceptibility: 2D numerical simulationsreveal that strain localization is closelyrelated to wave amplification in surficiallow-velocity layers, particularly belowconvex surface morphologies.     

强震作用下层状岩质斜坡破坏的大型振动台试验研究

以“5.12”汶川地震背景为基础,对不同岩性组合水平层状岩质边坡进行了大型振动台模型试验研究.文中介绍了振动台试验过程,通过对试验现象的观察处理分析,考察了边坡的破坏现象特征、变形破坏模式.实验结果表明,地震力作用下边坡的变形破坏程度、特征以及稳定性不仅与地震波的类型、加载方向、频率、振幅有关,同时还受边坡的岩体性质、...     

Spatial and temporal dynamics of soil moisture after rainfall events along a slope in Regosols of southwest China

High frequency seasonal drought in purple soils (Regosols in FAO taxonomy) of the hilly upland areas of Sichuan basin, China, is one of the key restrictive factors for crop production. In order to manage irrigation and fertilizer application in these soils effectively, the soil water content in a sloped plot with 60 cm soil depth was measured by neutron probe devices to investigate the soil moisture regime during the 1998 rainy season after various amounts of rainfall events. The results showed that variation of soil moisture along the slope positions was highest in the top soil layer during the period of sporadic rainfall that did not induce any runoff. The coefficients of variation of soil moisture at various slope positions (upper, middle, and lower) are 17·36%, 8·95%, 10·25%, 8·58%, 8·05% and 9·21% at the 10 cm, 20 cm, 30 cm, 40 cm, 50 cm and 60 cm soil depths respectively. When surface runoff occurred, the soil moisture dynamics at various positions on the plot were then very different. Soil water content decreased more rapidly on the upper slope than on the middle and lower slope positions. When both surface runoff and throughflow occurred, the soil moisture dynamics in the various layers showed a stable period (soil water content is near constant as time elapses) that lasted about 1 week. Also, the pattern of moisture dynamics is ‘decreasing–stabilization–decreasing’. Thus, irrigation and fertilization management according to the spatial and temporal features of soil moisture dynamics on sloped land can increase the water and fertilizer utilization efficacy by reducing their losses during the stable period. Copyright © 2006 John Wiley & Sons, Ltd.     

Analysis of soil moisture dynamics beneath olive trees

Olive cultivation is a widespread land use in Mediterranean climates. The proper implementation of soil and water conservation practices in groves requires detailed knowledge of the governing hydrological processes. In this work topsoil moisture dynamics under wet and dry conditions and across a small catchment was investigated in the inter row (IR) and directly under the olive tree canopies (UC). We do this using a sensor network (11 stations) and a simple bucket model which was calibrated (June, 2011–2012) and validated (June, 2012–2013). During most of the year the normalized soil moisture contents (s) were greater in the IR than under UC, with an average normalized soil moisture difference of 0.12. The difference between UC and IR normalized soil moisture followed a seasonal pattern, reaching a maximum near 0.30 during spring. An analysis of the normalized soil moisture probability density functions (pdfs) was bimodal, showing characteristic dominant wet and dry soil moisture states, with the highest probability densities for the dry state. Overall the spatial variability of soil moisture was lower UC than in the IR. This was a result of the soil moisture buffering capacity of the canopy with respect to rainfall and evaporation, in addition to observed differences in soil properties. Hourly soil moisture data were successfully modelled (R² > 0.85), both UC and in the IR, yet with the inclusion of a simple formulation for canopy interception for the former. The results provide insight into how olive trees change hydrological processes in their neighbourhood. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of rainfall and slope on runoff,soil erosion and rill development: an experimental study using two loess soils

Runoff generation and soil loss from slopes have been studied for decades, but the relationships among runoff, soil loss and rill development are still not well understood. In this paper, rainfall simulation experiments were conducted in two neighbouring plots (scale: 1 m by 5 m) with four varying slopes (17.6%, 26.8%, 36.4% and 46.6%) and two rainfall intensities (90 and 120 mm h^?1) using two loess soils. Data on rill development were extracted from the digital elevation models by means of photogrammetry. The effects of rainfall intensity and slope gradient on runoff, soil loss and rill development were different for the two soils. The runoff and soil loss from the Anthrosol surface were generally higher than those from the Calcaric Cambisol surface. Higher rainfall intensity produced less runoff and more sediment for almost each treatment. With increasing slope gradient, the values of cumulative runoff and soil loss peaked, except for the treatments with 90 mm h^?1 rainfall on the slopes with Anthrosol. With rainfall duration, runoff discharge decreased for Anthrosol and increased for Calcaric Cambisol for almost all the treatments. For both soils, sediment concentration was very high at the onset of rainfall and decreased quickly. Almost all the sediment concentrations increased on the 17.6% and 26.8% slopes and peaked on the 36.4% and 46.6% slopes. Sediment concentrations were higher on the Anthrosol slopes than on the Calcaric Cambisol slopes. At 90 mm h^?1 rainfall intensity, increasingly denser rills appeared on the Anthrosol slope as the slope gradient increased, while only steep slopes (36.4% and 46.6%) developed rills for the Calcaric Cambisol soil. The contributions of rill erosion ranged from 36% to 62% of the cumulative soil losses for Anthrosol, while the maximum contribution of rill erosion to the cumulative soil loss was only 37.9% for Calcaric Cambisol. Copyright © 2014 John Wiley & Sons, Ltd.     

Enhanced rock‐slope failure following ice‐sheet deglaciation: timing and causes

The temporal pattern of rock‐slope failures (RSFs) following Late Pleistocene deglaciation on tectonically stable terrains is controversial: previous studies variously suggest (1) a rapid response due to removal of supporting ice (‘debuttressing’), (2) a progressive decline in RSF frequency, and (3) a millennial‐scale delay before peak RSF activity. We test these competing models through beryllium‐10 (¹⁰Be) exposure dating of five closely‐spaced quartzite RSFs on the Isle of Jura, Scotland, to establish the relationship between timing of failure and those of deglaciation, episodes of rapid warming and periods of rapid glacio‐isostatic uplift. All five dated RSFs occurred at least 720–2240 years after deglaciation, with the probability of failure peaking ~2 ka after deglaciation, consistent with millennial‐scale delay model (3). This excludes debuttressing as an immediate cause of failure, though it is likely that time‐dependent stress release due to deglacial unloading resulted in progressive development of failure planes within the rock. Thaw of permafrost ice in joints is unlikely to have been a prime trigger of failure as some RSFs occurred several centuries after the onset of interstadial warming. Conversely, the timespan of the RSFs coincides with the period of maximum glacio‐isostatic crustal uplift, suggesting that failure was triggered by uplift‐driven seismic events acting on fractured rock masses. Implications of this and related research are: (1) that retreat of the last Pleistocene ice sheets across tectonically‐stable mountainous terrains was succeeded by a period of enhanced rock‐slope failure due to deglacial unloading and probably uplift‐driven seismicity; (2) that the great majority of RSFs in the British Isles outside the limits of Loch Lomond Stadial (= Younger Dryas) glaciation are of Lateglacial (pre‐Holocene) age; and (3) numerous RSFs must also have occurred inside Loch Lomond Stadial (LLS) glacial limits, but that runout debris was removed by LLS glaciers. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of soil moisture hysteresis on the functioning of capillary barriers

Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 John Wiley & Sons, Ltd.     

Stochastic modelling of soil moisture dynamics in a grassland of Qilian Mountain at point scale

Stochastic modeling of soil moisture dynamics is crucial to the quantitative understanding of plant responses to water stresses,hydrological control of nutrient cycling processes,water competition among plants,and some other ecological dynamics,and thus has become a hotspot in ecohydrology at present.In this paper,we based on the continuously monitored data of soil moisture during 2002―2005 and daily precipitation date of 1992―2006,and tried to make a probabilistic analysis of soil moisture dynamics at point scale in a grassland of Qilian Mountain by integrating the stochastic model improved by Laio and the Monte Carlo method.The results show that the inter-annual variations for the soil moisture patterns at different depths are very significant,and that the coefficient of variance(CV) of surface soil moisture(20 cm) is almost continually kept at about 0.23 whether in the rich or poor rainy years.Interestingly,it has been found that the maximal CV of soil moisture has not always appeared at the surface layer.Comparison of the analytically derived soil moisture probability density function(PDF) with the statistical distribution of the observed soil moisture data suggests that the stochastic model can reasonably describe and predict the soil moisture dynamics of the grassland in Qilian Mountain at point scale.By extracting the statistical information of the historical precipitation data in 1994―2006,and inputting them into the stochastic model,we analytically derived the long-term soil moisture PDF without considering the inter-annual climate fluctuations,and then numerically derived the one when considering the inter-annual fluctuation effects in combination with a Monte-Carlo procedure.It was found that,though the peak position of the probability density distribution significantly moved towards drought when considering the disturbance forces,and its width was narrowed,accordingly its peak value was increased,no significant bimodality was observed in the soil moisture dynamics under the given intensity of random fluctuation disturbance.     

Spatio-temporal evolution and time-stable characteristics of soil moisture within remote sensing footprints with varying soil, slope, and vegetation

Air-borne passive microwave remote sensors measure soil moisture at the footprint scale, a scale of several hundred square meters or kilometers that encompasses different characteristic combinations of soil, topography, vegetation, and climate. Studies of within-footprint variability of soil moisture are needed to determine the factors governing hydrologic processes and their relative importance, as well as to test the efficacy of remote sensors. Gridded ground-based impedance probe water content data and aircraft-mounted Electronically Scanned Thinned Array Radiometer (ESTAR) pixel-average soil moisture data were used to investigate the spatio-temporal evolution and time-stable characteristics of soil moisture in three selected (LW03, LW13, LW21) footprints from the Southern Great Plains 1997 (SGP97) Hydrology Experiment. Better time-stable features were observed within a footprint containing sandy loam soil than within two pixels containing silty loam soil. Additionally, flat topography with split wheat/grass land cover produced the largest spatio-temporal variability and the least time stability in soil moisture patterns. A comparison of ground-based and remote sensing data showed that ESTAR footprint-average soil moisture was well calibrated for the LW03 pixel with sandy loam soil, rolling topography, and pasture land cover, but improved calibration is warranted for the LW13 (silty loam soil, rolling topography, pasture land) and LW21 (silty loam soil, flat topography, split vegetation of wheat and grass land with tillage practice) pixels. Footprint-scale variability and associated nonlinear soil moisture dynamics may prove to be critical in the regional-scale hydroclimatic models.