Coupling between hydrogeology and deformation of mountainous rock slopes: Insights from La Clapière area (southern Alps,France)

Meteoric infiltration influence on large mountainous rock slopes stability is investigated by comparing hydrogeologic and gravitational structures from detailed mapping of the ‘La Clapière’ slope. The slope infiltrated waters are trapped in a perched aquifer that is contained in deposits inside tensile cracks of the upper part of the slope. Flow rates of 0.4 to 0.8 l?s^?1 from the perched aquifer to the landslide cause landslide accelerations. Numerical modeling shows that a 0.75 l?s^?1 infiltration yield increases conditions for toppling with failure through tilting of large rock volumes from the perched aquifer bottom down to the foot of the slope. To cite this article: Y. Guglielmi et al., C. R. Geoscience 337 (2005).     

Temporal evolution of weathered cataclastic material in gravitational faults of the La Clapiere deep-seated landslide by mechanical approach

After a few years of research, the observation and the analysis of the deep-seated landslides suggest that these are mainly controlled by tectonic structures, which play a dominant role in the deformation of massif slopes. The La Clapière deep-seated landslide (Argentera Mercantour massif) is embedded in a deep-seated gravitational slope deformation affecting the entire slope, and characterized by specific landforms (trenches, scarps??). Onsite, the tangential displacement direction of the trenches and the scarps are controlled by the tectonic structures. The reactivation of the inherited fault in gravitational faults create a gouge material exposed to an additional mechanical and chemical weathering as well as an increased of leaching. The displacement of these reactivated faults gets increasingly important around the area of the La Clapière landslide and this since 3.6?ka BP. In this study, mechanical analysis and grain size distributions were performed and these data were analysed according to their proximity the La Clapiere landslide and times of initiation of the landslide by ¹⁰Be dating. Triaxial test results show that the effective cohesion decreases and the effective angle of internal friction increases from the unweathered area to the weathered area. The whole distribution of the grain size indicates that the further the shear zone is open or developed, the further the residual material loses its finest particles. This paper suggests that the mechanical evolution along the reactivated fault is influenced by the leaching processes. For the first time, we can extract from these data temporal behaviour of the two main mechanical parameters (cohesion and angle of internal friction) from the beginning of the La Clapiere landslide initiation (3.6 ka BP) to now.     

Evaluating the susceptibility of landslide landforms in Japan using slope stability analysis: a case study of the 2016 Kumamoto earthquake

This study analyzed 267 landslide landforms (LLs) in the Kumamoto area of Japan from the database of about 0.4 million LLs for the whole of Japan identified from aerial photos by the National Research Institute for Earth Science and Disaster Resilience of Japan (NIED). Each LL in the inventory is composed of a scarp and a moving mass. Since landslides are prone to reactivation, it is important to evaluate the sliding-recurrence susceptibility of LLs. One possible approach to evaluate the susceptibility of LLs is slope stability analysis. A previous study found a good correlation (R ² = 0.99) between the safety factor (F _s ) and slope angle (α) of F _s  = 17.3α ^?0.843. We applied the equation to the analysis of F _s for 267 LLs in the area affected by the 2016 Kumamoto earthquake (M _j  = 7.3). The F _s was calculated for the following three cases of failure: scarps only, moving mass only, and scarps and moving mass together. Verification with the 2016 Kumamoto earthquake event shows that the most appropriate method for the evaluation of LLs is to consider the failure of scarps and moving mass together. In addition, by analyzing the relationship between the factors of slope of entire landslide and slope of scarp for LLs and comparing the results with the Aso-ohashi landslide, the largest landslide caused by the 2016 Kumamoto earthquake, we also found that morphometric analysis of LLs is useful for forecasting the travel distance of future landslides.     

Geometrical characteristics of earthquake-induced landslides and correlations with control factors: a case study of the 2013 Minxian,Gansu, China, <Emphasis Type="Italic">Mw</Emphasis> 5.9 event

Geometric parameters are useful for characterizing earthquake-triggered landslides. This paper presents a detailed statistical analysis on this issue using the landslide inventory of the 2013, Minxian, China Mw 5.9 earthquake. Based on GIS software and a 5-m resolution DEM, geometric parameters of 635 coseismic landslides (with areas larger than 500 m²) were obtained, including height, length, width, reach angle (arc tangent of the height-length ratio), and aspect ratio (length-width ratio). The fitting relationship of height and length from these data is H = 0.6164L + 0.4589, with an average reach angle of 31.65°. The landslide aspect ratios concentrate in the range of 1.4～2.6, with an average of 2.11. According to the plane geometric shapes and aspect ratios, the landslides are classified into four categories: transverse landslide (LA1, L/W ≤ 0.8), isometric landslide (LA2, 0.8 < L/W ≤ 1.2), longitudinal landslide (LA3, 1.2 < L/W ≤ 3), and elongated landslide (LA4, L/W > 3). Statistics of these four types of landslides versus ten classified control factors (elevation, slope angle, slope aspect, curvature, slope position, distance to drainages, lithology, seismic intensity, peak ground acceleration, and distance to seismogenic fault) are used to examine their possible correlations and the landslide-prone areas, which would be helpful to the landslide disaster mitigation in the affected area.     

Geophysical survey to estimate the 3D sliding surface and the 4D evolution of the water pressure on part of a deep seated landslide

Geophysical surveys were conducted on the very unstable front part of the La Clapière landslide in the French Alps (Alpes Maritimes). The electrical resistivity survey was carried out to obtain, for the first time on this deep-seated landslide, 3D information on the slipping surface and the vertical drained faults. Moreover, we planned to follow within time (6 months) the evolution of the saturated zones (presence of gravitational water) and their percolation into the shearing zones. Our 4D results showed the importance of the complex water channelization within the slope and relation to geological discontinuities.     

Back analysis of a large landslide in a flysch rock mass

Flysch is a sedimentary rock consisting of a rhythmic alternation of hard (limestone, sandstone, siltstone) and weak (marl, mudstone, claystone) layers. Because of the presence of layers with different physical properties, the mechanical characterization of heterogeneous rock masses such as flysch is a real challenge. Different methods have been proposed in the literature to characterize flysch, combining empirical classification indexes with laboratory tests. Most of these methods, however, were specifically designed for tunneling and underground excavations, and their applicability to slope stability problems is not yet fully investigated. In this study, we analyze a large landslide in a cretaceous flysch rock in order to compare the mobilized strength at failure with those predicted by the modified GSI method (Marinos and Hoek, 2001). The landslide occurred in the Savena River basin (Northern Apennines of Italy) on April 6, 2013, with a volume of about 3 million m³. Soon after the failure, geological, geotechnical, and geophysical investigations were carried out to detect the failure mechanism and define the landslide geometry. Back analyses of the failed slope were performed using both limit equilibrium and finite difference methods to estimate the in situ strength of the flysch. The results show that the mobilized rock mass cohesion is very low (c ^'?≈?20?÷?40 kPa) and that the modified GSI method can predict the in situ strength only assuming a disturbance factor D = 1. Moreover, the analysis shows that the linearization criteria proposed in literature to compute the equivalent Mohr-Coulomb parameters remarkably overestimate the rock mass strength.     

The formation and development of debris flows in large watersheds after the 2008 Wenchuan Earthquake

The Wenchuan earthquake has caused abundance of loose materials supplies for debris flows. Many debris flows have occurred in watersheds in area beyond 20 km², presenting characteristics differing from those in small watersheds. The debris flows yearly frequency decreases exponentially, and the average debris flow magnitude increases linearly with watershed size. The rainfall thresholds for debris flows in large watersheds were expressed as I?=?14.7 D ^?0.79 (2 h?<?D?<?56 h), which is considerably higher than those in small watersheds as I?=?4.4 D ^?0.70 (2 h?<?D?<?37 h). A case study is conducted in Ergou, 39.4 km² in area, to illustrate the formation and development processes of debris flows in large watersheds. A debris flow develops in a large watershed only when the rainfall was high enough to trigger the wide-spread failures and erosions on slope and realize the confluence in the watershed. The debris flow was supplied by the widely distributed failures dominated by rill erosions (14 in 22 sources in this case). The intermittent supplying increased the size and duration of debris flow. While the landslide dam failures provided most amounts for debris flows (57 % of the total amount), and amplified the discharge suddenly. During these processes, the debris flow velocity and density increased as well. The similar processes were observed in other large watersheds, indicating this case is representative.     

Volume estimation and stage division of the Mahu landslide in Sichuan Province,China

The Mahu lake, the third deepest lake in China, is located on the west bank of the Jinsha River in Leibo county, Sichuan Province. It is a dammed lake created by an old landslide on the ancient Huanglang river, a tributary on the west bank of the Jinsha River. Previous studies (Wang and Lu in J Mt Res S1:44–47, 2000) suggested that this landslide was caused by an earthquake approximately 372 ka (Middle Pleistocene), during which a few hundreds of million cubic meters of debris were deposited between 1177 and 900 m a.s.l. (above sea level), covering an area of around 15 km². Our further investigations, including geodetic survey, borehole drilling, and field reconnaissance, combining with five chronological data, have made some new discoveries at this site. First, the toe of the landslide extends from 900 m a.s.l. down to 320 m a.s.l., i.e., the local bed elevation of the contemporary Jinsha River. Second, the area of the landslide deposits is 17.3 km² with a volume of 2.38 km³, much larger than the previous estimation. Thus, it should be one of the largest known landslides in China. And the lower elevation of the landslide’s toe also rules out the possibility that it is a hanging valley on the ancient Huanglang river. Our work suggests that this landslide was created by five events according to the overlapping characteristics of the deposits and five chronological data, which are old than 52,600 years, old than 16,000 years, old than 15,500 years, 5800 years, and old than 4200 years, respectively.     

Geological Characteristics and Stability Evaluation of Wanjia Middle School Slope in Wenchuan Earthquake Area

To investigate the formation mechanism and the stability of Wanjia middle school slope in Wenchuan Earthquake Area, the macroscopic geological characteristics and the failure process of the landslide are researched by engineering geology analysis method, limit equilibrium method, and finit element method. The results show that after the Wenchuan Earthquake, retaining walls, houses and other infrastructure on the foot of Wanjia middle school slope were severely destroyed, 10 cm wide tension fracture appeared at the trailing edge of the slope. Wanjia middle school slope is a type of medium-sized soil landslide. The area of the deformation body is about 19,314 m², the total volume of the deformation body is about 23 × 10⁴ m³. There may be two potential sliding surfaces in the unstable slope: shallow and deep landslide. The analysis results of the limit equilibrium method and the finite element method show that: under dead weight, dead weight + rainstorm, dead weight + earthquake conditions, the plastic zone occurs mainly at the middle part or the trailing edge of the slope, and it doesn’t fully cut through the deep landslide body, so the deep landslide is stable. However, under rainstorm or earthquake conditions, the plastic zone almost completely cut thorough the shallow landslide body, it shows that the shallow landslide is in the understable–basic stable state. It is found that the results of finite element method is concordant with the results of the limit equilibrium method (F _s = 1.06–1.29, the shallow landslide is in the basic stable–stable state). The calculation results show that shallow landslides are likely to occur in Wanjia middle school slope during a rainstorm or an earthquake, so monitoring and control of the slope should be strengthened. The shallow landslide should be managed by some measures, such as anti slide pile retaining structures and drainage works, and the dangerous rock bodies on the slope surface should be cleaned up.     

Setting up evaluate indicators for slope control engineering based on spatial clustering analysis

Rainfall-induced landslides (RILs) have been a source of social and economic disruption in the mountainous Baguio area in northern Philippines. Prolonged heavy rainfall usually happens during tropical cyclone and southwest monsoon activity. A pragmatic approach to RIL mitigation is to develop rainfall-based early warning. We implemented a modified regression method to derive the empirical minimum intensity (I)–duration (D) threshold I = 6.46 D ^?0.28 and a normalized ID threshold NI = 0.002 D ^?0.28 for rainfall duration ranging between 24 and 264 h. Using a separate data set to evaluate the applicability of the threshold, 93% of the landslide-triggering rainfall events fell above the derived threshold. RILs also occurred when 24-h rainfall was 0.02–28% of the mean annual precipitation or after accumulating at least 500 mm of rainfall from the onset of the rainy season. The thresholds may be further refined as more landslide data become available in the future.     

Determination of natural radioactivity levels in soil and travertine of the region of Tokat and Sivas,Turkey

In this study, the environmental radioactivity measurements for Tokat and Sivas provinces in the northeast of Turkey were performed. Using gamma ray spectrometry, the activity concentrations of natural radionuclides in soil and travertine samples (²³²Th, ²²⁶Ra, and ⁴⁰K) were determined. The annual effective dose equivalent, the absorbed doses rate in air, the radium equivalent, and the external hazard index were obtained from these activities. The activity concentrations vary from 9.09 to 17.04 Bq kg^?1 for ²³²Th, from 36.53 to 76.95 Bq kg^?1 for ²²⁶Ra, and from 216.56 to 576.59 Bq kg^?1 for ⁴⁰K in soil samples. The activity concentrations in travertines vary from 15.99 to 21.01 Bq kg^?1 for ²³²Th, from 19.89 to 67.71 Bq kg^?1 for ²²⁶Ra, and from 179.89 to 314.43 Bq kg^?1 for ⁴⁰K. The average dose rate in air for soil and travertine samples was 43.41 and 41.05 nGy h^?1 respectively. The obtained results are presented and compared with other studies, and the results of this study are lower than the international recommended value (55 nGy h^?1) given by UNSCEAR ( 2000). The results show that the region has a background radiation level within the natural limits.     

Modeling and mapping of water erosion in northeastern Algeria using a seasonal multicriteria approach

Aurès region remains one of the most exposed areas to water erosion phenomenon in Algeria, because of the strong climatic aggressiveness, the rugged relief, the predominance of sensitive land, and a vegetative cover that does not play its protective role. This article is a part of studies performed to protect agricultural and water infrastructure in this region. The main objective of this study is the cartographic modeling of an erosion hazard at the Oued Chemoura watershed, representative of the Aurès. The modeling approach uses a geographic information system and incorporates the following six factors controlling erosion: slope, friability of substrate, erodibility of soils, land cover, rainfall erosivity, and support practices. Result shows a synthetic map of the soil erosion hazard which locates the most threatened areas and priorities for possible planning interventions. A statistical study on the relationship solid–liquid flow was developed. Measurements conducted at the station of Chemoura, over the period 1969–1994, were exploited for this purpose. The results show a high specific degradation varying between 50 and 360 tons km^?2 season^?1.     

Precursors of instability in a natural slope due to rainfall: a full-scale experiment

A full-scale landslide-triggering experiment was conducted on a natural sandy slope subjected to an artificial rainfall event, which resulted in mobilisation of 130 m³ of soil mass. Novel slope deformation sensors (SDSs) were applied to monitor the subsurface pre-failure movements and the precursors of the artificially triggered landslide. These fully automated sensors are more flexible than the conventional inclinometers by several orders of magnitude and therefore are able to detect fine movements (<?1 mm) of the soil mass reliably. Data from high-frequency measurements of the external bending work, indicating the transmitted energy from the surrounding soil to these sensors, pore water pressure at various depths, horizontal soil pressure and advanced surface monitoring techniques, contributed to an integrated analysis of the processes that led to triggering of the landslide. Precursors of movements were detected before the failure using the horizontal earth pressure measurements, as well as surface and subsurface movement records. The measurements showed accelerating increases of the horizontal earth pressure in the compression zone of the unstable area and external bending work applied to the slope deformation sensors. These data are compared to the pore water pressure and volumetric water content changes leading to failure.     

Thermoelastic properties of chromium oxide Cr<Subscript>2</Subscript>O<Subscript>3</Subscript> (eskolaite) at high pressures and temperatures

A new synchrotron X-ray diffraction study of chromium oxide Cr₂O₃ (eskolaite) with the corundum-type structure has been carried out in a Kawai-type multi-anvil apparatus to pressure of 15 GPa and temperatures of 1873 K. Fitting the Birch–Murnaghan equation of state (EoS) with the present data up to 15 GPa yielded: bulk modulus (K _0,T0), 206 ± 4 GPa; its pressure derivative K′_0,T , 4.4 ± 0.8; (?K _0,T /?T) = ?0.037 ± 0.006 GPa K^?1; a = 2.98 ± 0.14 × 10^?5 K^?1 and b = 0.47 ± 0.28 × 10^?8 K^?2, where α _0,T  = a + bT is the volumetric thermal expansion coefficient. The thermal expansion of Cr₂O₃ was additionally measured at the high-temperature powder diffraction experiment at ambient pressure and α _0,T0 was determined to be 2.95 × 10^?5 K^?1. The results indicate that coefficient of the thermal expansion calculated from the EoS appeared to be high-precision because it is consistent with the data obtained at 1 atm. However, our results contradict α ₀ value suggested by Rigby et al. (Brit Ceram Trans J 45:137–148, 1946) widely used in many physical and geological databases. Fitting the Mie–Grüneisen–Debye EoS with the present ambient and high-pressure data yielded the following parameters: K _0,T0 = 205 ± 3 GPa, K′_0,T  = 4.0, Grüneisen parameter (γ ₀) = 1.42 ± 0.80, q = 1.82 ± 0.56. The thermoelastic parameters indicate that Cr₂O₃ undergoes near isotropic compression at room and high temperatures up to 15 GPa. Cr₂O₃ is shown to be stable in this pressure range and adopts the corundum-type structure. Using obtained thermoelastic parameters, we calculated the reaction boundary of knorringite formation from enstatite and eskolaite. The Clapeyron slope (with \({\text{d}}P/{\text{d}}T = - 0.014\) GPa/K) was found to be consistent with experimental data.     

Structural control on the kinematics of the deep-seated La Clapière landslide revealed by L-band InSAR observations

The objective of this work is to document the deformation pattern of the deep-seated La Clapière landslide for the period 2007–2010 from the combination of L-band synthetic aperture radar (SAR) interferograms, ground-based total station measurements and identification of the slope geomorphological structures. The interferograms are calculated for pairs of ALOS/PALSAR images at a time interval of 46 days. The displacement field derived from the interferograms reveals a non-uniform displacement gradient from the top (subsidence) to the bottom (accumulation). Vertical velocities are calculated from the unwrapped phase values and are in good agreement with ground-based measurements. The results demonstrate the potential of L-band ALOS/PALSAR imagery for the monitoring of active landslides characterized by complex kinematic patterns and by important changes in the soil surface backscattering in time.     

Landslide susceptibility analysis based on ArcGIS and Artificial Neural Network for a large catchment in Three Gorges region,China

A landslide susceptibility map is very important and necessary to efficiently prevent and mitigate the losses brought by natural hazard for a large area. For the purpose of landslide susceptibility analysis for the whole Xiangxi catchment (3,209 km²), Artificial Neural Network (ANN) analysis was applied as the main method. The whole catchment was divided into two parts: the training area and the implementation area. The backwater area (559 km²) of Xiangxi catchment was used as the training area for the ANN method. In the training area the correlations between the landslide distribution and its causative factors, which includes lithology, slope angle, slope curvature and river network, have been analyzed based on the geological map and digital elevation model (DEM). The back-propagation training algorithm in ANN was selected to train the sample data from the training area, which were composed of input data (causative factors) and target output data (landslide occurrence), in order to find the correlations between them. Based on these correlations and input data in the implementation area (causative factors), the network output data were obtained for the implementation area. In the end, a map of landslide susceptibility, which was established by network output data, was presented for Xiangxi catchment. ArcGIS was applied to extract and quantify input information from a DEM for susceptibility analysis and also to present the result visually. As a result, a landslide susceptibility map, in which 70 % of all landslides are rightly classified in the training area (backwater area), was created for Xiangxi catchment.     

Groundwater flow dynamics of weathered hard-rock aquifers under climate-change conditions: an illustrative example of numerical modeling through the equivalent porous media approach in the north-western Pyrenees (France)

A numerical groundwater model of the weathered crystalline aquifer of Ursuya (a major water source for the north-western Pyrenees region, south-western France) has been computed based on monitoring of hydrological, hydrodynamic and meteorological parameters over 3 years. The equivalent porous media model was used to simulate groundwater flow in the different layers of the weathered profile: from surface to depth, the weathered layer (5?·?10^?8?≤?K?≤?5?·? 10^?7 m s^?1), the transition layer (7?·?10^?8?≤?K?≤?1?·? 10^?5 m s^?1, the highest values being along major discontinuities), two fissured layers (3.5?·?10^?8?≤?K?≤?5?·?? 10^?4 m s^?1, depending on weathering profile conditions and on the existence of active fractures), and the hard-rock basement simulated with a negligible hydraulic conductivity (K = 1 10 ^?9 ). Hydrodynamic properties of these five calculation layers demonstrate both the impact of the weathering degree and of the discontinuities on the groundwater flow. The great agreement between simulated and observed hydraulic conditions allowed for validation of the methodology and its proposed use for application on analogous aquifers. With the aim of long-term management of this strategic aquifer, the model was then used to evaluate the impact of climate change on the groundwater resource. The simulations performed according to the most pessimistic climatic scenario until 2050 show a low sensitivity of the aquifer. The decreasing trend of the natural discharge is estimated at about ?360 m³ y^?1 for recharge decreasing at about ?5.6 mm y^?1 (0.8 % of annual recharge).     

Assessment of natural radioactivity,associated radiological health hazards indices and soil-to-crop transfer factors in cultivated area around a fertilizer factory in Onne,Nigeria

Concentration of natural radionuclides in three major staple food crops cultivated around a fertilizer plant in Onne, Rivers State Nigeria and the cultivated soil samples were determined using gamma spectroscopy operated on a Canberra vertical high purity 3″ × 3″ NaI(TI) detector. The average activity concentration of ²³⁸U, ²³²Th and ⁴⁰K was determined, for cassava flour (U 19.3 ± 5.0, Th 11.4 ± 3.3, K 426.9 ± 33.8) Bq kg^?1, for yam flour (U 6.3 ± 1.8, Th 8.4 ± 2.6, K 227.0.9 ± 27.3) Bq kg^?1 while for cocoyam flour (U 7.5 ± 2.7, Th 7.1 ± 2.3, K 195.8 ± 25.83) Bq kg^?1. The mean activity concentration for soil samples is 18.7 ± 3.7 Bq kg^?1, 18.0 ± 3.8 Bq kg^?1 and 308.4 ± 22.4 Bq kg^?1 for ²³⁸U, ²³²Th and ⁴⁰K, respectively. These values obtained show enhanced ⁴⁰K concentration which is attributed to the effluent discharge from a fertilizer plant and its applications to farmlands, but ²³⁸U, ²³²Th values are well within the global average and values reported in some regions and countries of the world. Radiation hazard indices obtained to estimate potential radiological health risk in both foodstuffs and soil samples are well below their permissible limit as set by UNSCEAR [Sources and effects of ionizing radiation (Report to the General Assembly), 2000]. The rate of radionuclides transfer from soil to crops was moderate with mean transfer factors of ²³²Th < ²³⁸U < ⁴⁰K.     

Dynamics of carbon fluxes with responses to vegetation,meteorological and terrain factors in the south-eastern Tibetan Plateau

Tibetan Plateau (TP) is the highest and most extensive plateau in the world and has been known as the roof of the world, and it is sensitive to climate change. The researches of CO₂ fluxes (F _C) in the TP region play a significant role in understanding regional and global carbon balance and climate change. Eddy covariance flux measurements were conducted at three sites of south-eastern TP comprising Dali (DL, cropland ecosystem), LinZhi (LZ, alpine meadow ecosystem) and Wenjiang (WJ, cropland ecosystem); amongst those DL and LZ are located in plateau region, while WJ is in plain region. Dynamics of F _C and influences of vegetation, meteorological (air temperature, photosynthetically active radiation, soil temperature and soil water content) and terrain factors (altitude) were analysed on the basis of data taken during 2008. The results showed that, in the cool sub-season (March, April, October and December), carbon sink appeared even in December with fluxes of (?0.021 to ?0.05) mg CO₂ m^?2 s^?1 and carbon source only in October (0.03 ± 0.0048) mg CO₂ m^?2 s^?1 in DL and WJ site. In LZ site, carbon sink was observed in April: (?0.036 ± 0.0023) mg CO₂ m^?2 s^?1 and carbon sources in December and March (0.008–0.010 mg CO₂ m^?2 s^?1). In the hot sub-season (May–August), carbon source was observed only in May with (0.011 ± 0.0022), (0.104 ± 0.0029) and (0.036 ± 0.0017) fluxes in LZ, DL and WJ site, respectively, while carbon sinks with (?0.021 ± 0.0041), (?0.213 ± 0.0007) and (?0.110 ± 0.0015) mg CO₂ m^?2 s^?1 fluxes in LZ, DL, and WJ, respectively. Comparing with plain region (WJ), carbon sinks in plateau region (DL and LZ) lasted for a longer time, and the absorption sum was large and up to (–357.718 ± 0.0054) and (?371.111 ± 0.0039) g C m^?2 year^?1, respectively. The LZ site had the weakest carbon sink with (?178.547 ± 0.0070) g C m^?2 year^?1. Multivariate analysis of covariance showed that altitude (AL) as an independent factor explained 39.5 % of F _C (P < 0.026). F _C had a quadratic relationship with Normalized difference vegetation index (NDVI) (R ² ranges from 0.485 to 0.640 for three sites), an exponential relationship with soil temperature at 5-cm depth (ST ₅) at night time and a quadratic relationship with air temperature (T _a) at day time. Path analysis indicated that photosynthetically active radiation (PAR), sensible heat fluxes (H) and other factors all had direct or indirect effects on F _C in all of the three tested sites around the south-eastern TP.     

Field monitoring of groundwater responses to heavy rainfalls and the early warning of the Kualiangzi landslide in Sichuan Basin,southwestern China

The Kualiangzi landslide was triggered by heavy rainfalls in the “red beds” area of Sichuan Basin in southwestern China. Differing from other bedrock landslides, the movement of the Kualiangzi landslide was controlled by the subvertical cracks and a subhorizontal bedding plane (dip angle < 10°). The ingress of rainwater in the cracks formed a unique groundwater environment in the slope. Field measurement for rainfall, groundwater movement, and slope displacement has been made for the Kualiangzi landslide since 2013. The field monitoring system consists of two rainfall gauges, seven piezometers, five water-level gauges, and two GPS data loggers. The equipments are embedded near a longitudinal section of the landslide, where severe deformation has been observed in the past 3 years. The groundwater responses to four heavy rainfall events were analyzed between June 16 and July 24 in 2013 coincided with the flood season in Sichuan. Results showed that both of the water level and the pore-water pressure increased after each rainfall event with delay in the response time with respect to the precipitation. The maximum time lag reached 35 h occurred in a heavy rainfall event with cumulative precipitation of 127 mm; such lag effect was significantly weakened in the subsequent heavy rainfall events. In each presented rainfall event, longer infiltration period in the bedrock in the upper slope increased the response time of groundwater, compared to that of in the gravels in the lower slope. A translational landslide conceptual model was built for the Kualiangzi landslide, and the time lag was attributed to the gradual formation of the uplift pressure on the slip surface and the softening of soils at the slip surface. Another important observation is the effect on the slope movement which was caused by the water level (H _w) in the transverse tension trough developed at the rear edge of the landslide. Significant negative correlation was found for H _w and the slope stability factor (F _s), in particular for the last two heavy rainfall events, of which the drastic increase of water level caused significant deterioration in the slope stability. The rapid drop (Δ?=?22.5 kPa) of pore-water pressure in the deep bedrock within 1 h and the large increase (Δ?=?87.3 mm) of surficial displacement were both monitored in the same period. In the end, a four-level early warning system is established through utilizing H _w and the displacement rate D _r as the warning indicators. When the large deformation occurred in flood season, the habitants at the leading edge of the landslide can be evacuated in time.