Study on the Spatially Variable Saturated Hydraulic Conductivity and Deformation Behavior of Accumulation Reservoir Landslide Based on Surface Nuclear Magnetic Resonance Survey

Saturated hydraulic conductivity (Ks) is spatially variable in accumulation landslide sites that exert significant effort onto landslide seepage and deformation behavior. To better understand spatial variability and the effect of Ks on the slide mass of an accumulation landslide, this study introduced the surface nuclear magnetic resonance (SNMR) technology to study a representative reservoir accumulation landslide field in the -ree Gorges Reservoir area (TGRA), the Baishuihe landslide, to obtain a series of relative reliable spatial measurements of Ks effectively on the basis of calibration in terms of the field tests measurements. -e estimated Ks values were distributed log-normally for the overall landslide mass site with a wide range of 3.00 × 10∼7.80 × 10 cm/s, which reaches about 3 orders of magnitude. Variogram analysis indicated that the Ks values have the range (A) of 295.89m and 65.56m for the overall site and major cross-sectional analysis, respectively. A finite-element seepage-stress analysis associated with a Kriging-interpolated spatial Ks variable calculation model based on the best-fitted theoretical variogram was subsequently performed to study the seepage and deformation behavior of the landslide. -e available monitored data and simulated results of the finite-element seepage-stress analysis indicated that the Baishuihe landslide is a progressive landslide, and the main factor influencing the deformation is rainfall and reservoir water fluctuation. -is study provides an unconventional framework for studying the heterogeneous geomaterial and contributes to a better understanding of the spatial variation of the hydraulic property of accumulation reservoir landslides at a field scale.


Introduction
Accompanied by the reservoir operation of hydropower projects, reservoir landslides have become a typical geohazard in reservoir areas in the world [1,2].e instability of the reservoir landslide may arouse devastating economic and social costs, even the loss of human life.Under the combined effect of reservoir water-level fluctuation and rainfall, saturated hydraulic conductivity (K s ) is a key soil hydraulic parameter, which controls the water seepage process, determines water distribution, and thus exerts a dominant impact on landslide deformation and instability [3][4][5][6].
It is generally accepted that K s is characterized by spatial variability at different scales due to the parent material and environmental process [7][8][9][10].However, K s values are usually assumed to be identical and are generally derived from experiments of a few samples.e variations of the values derived from the characteristics of discontinuity, heterogeneity, and anisotropy have not been taken into account in this assumption [11].In this sense, the understanding of spatial variability of hydrological properties

Basics of SNMR.
e magnetic spin is an inherent property of an atom that possesses an angular momentum, with an associated magnetic moment.Under the action of the stable earth's stable magnetic field, the hydrogen protons spin continuously in a precess motion about the magnetic field axis at the characteristic proton Larmor frequency, which is proportional to the magnetogyric ratio of the intensity of the earth's magnetic field and atomic nucleus.e orientation of the static field and the plane of precession perpendicular to it are commonly referred to as longitudinal direction and transverse direction, respectively [18].During the process of the SNMR survey, an alternating current with the resonance of Larmor frequency is supplied to the large surface-based wire-loop transmitter, which is typically a circular or square loop of 50-100 m in diameter/edge length.A slight macroscopic spin magnetization moment is produced at the same time by the excitation pulse of the oscillating current.When the current pulse is terminated, magnetic resonance signals produced by the excitation of different pulse moments can be explored by the same surface-based transmitting coil.After an instrument delay known as "dead time" (τ dead ), the magnetic resonance signal can be measured and expressed as follows: where E 0 is the initial amplitude of the signal, which can be obtained by the measured amplitude E 0d and the relaxation time T * 2 as E 0 � E 0d exp(τ dead /T * 2 ), and f 0 and φ 0 denote the initial frequency and initial phase, respectively.e phase φ 0d correlates with the electrical conductivity of the geomaterials, but it is not currently used for SNMR data interpretation.With regard to the relaxation time, it is believed that T * 2 is susceptible to the disturbance of electromagnetic noise; therefore, T * 1 is adopted for hydrological parameter interpretation [18].e strength of the signal and rate of decay are related to the magnitude of hydrogen protons in groundwater molecules.In other words, the amplitude of the magnetic resonance signal is proportional to the free water content of the surveyed area [24].
e initial signal amplitude E 0d is proportional to the number of protons in the measured sample, which is equivalent to the volume of water in the measured sample [25], so the SNMRestimated effective porosity (∅) can be determined from E 0d using an instrument-specific calibration factor.
With regard to hydraulic conductivity estimation, Schlumberger-Doll Research (SDR) models give an equation to estimate the permeability k SDR [26]: 2 Advances in Civil Engineering where b, n, and m are the empirical constants and T 2(1)ML is used to present the average surface-area-to-value ratio of the pore of the entire saturated pore space, which can be replaced by T 2(1) or T * 2 (1) .Under the condition of the unconsolidated near-surface aquifer material as we measured, n � 2 and m � 1, b is regarded as a turning parameter to obtain the best possible match between the measurement and estimation [27][28][29].erefore, hydraulic conductivity K can be related to k by the equation K � ((ρg)/μ)k, where g � 9.8 m/s 2 is the gravitational acceleration, μ is the temperature-dependent dynamic viscosity, and ρ is the temperature-dependent density, which can be determined in terms of table lookup according to the temperature.
For all of the analyses above, the physical constants can be absorbed into the constant b ′ in the SDR equation as follows: where b ′ should be calibrated based on the temperature as well as the drilling data and other forms of measurement on core samples for the requirement of accuracy [29].

Background.
e Baishuihe landslide is a typical accumulation reservoir landslide in TGRA.It is located on the south bank of the Yangtze River, 56 km west of the ree Gorges Dam (Figures 1(a) and 1(b)), with the main sliding direction of 20 °with respect to the north.It produced remarkable deformation since the first impoundment of the ree Gorges Reservoir in 2003 and has a great influence on the shipping safety of Yangtz River.
e multiyear timedisplacement data of the landslide surface obtained from GPS monitoring exhibited a step-like pattern, which aroused much research attention [30][31][32][33].However, most previous research studies aimed at predicting the deformation behavior and paid little attention to the engineering geological model and landslide mechanism.
e sophisticated engineering surveys involving 23 boreholes with a total length of 699.20 m and 14 exploratory trenches with a total volume of 101.24 m 3 (shown in Figure 1(c)) have been carried out by Sanxia University in 2012 and China university of Geosciences, Wuhan, in 2013 [34], respectively.Based on the revealed observations, the boundary of the Baishuihe landslide had been determined as the predetermined warning area of the landslide.As shown in Figure 1(b), the morphology of the landslide exhibited an irregular gentle concave shape with the elevation from 120 m to 297 m and the volume of 645 × 10 4 m 3 .To date, multiple monitoring devices had been arranged in this typical landslide, involving GPS, borehole clinometer, SNMR, and distributed optical fiber sensor deformation monitoring, as shown in Figure 1(c).e GPS therein was arranged in 2007 by the government, but the monitoring data after 2012 are still unpublicized till now.

Geological Setting.
e Baishuihe landslide is developed in the clastic rocks of the Lower Jurassic Xiangxi Formation (J 1x ), with the specific lithology of dark gray thin-bedded and medium-bedded siltstone mingling with thin stratified mudstone and coalbed.e inclination direction is 15 °∼20 °, coinciding with the main sliding direction.erefore, the landslide is a monoclinic consequent bedding landslide along the weak intercalated layer initially.However, GPS monitoring data and in situ observations indicated that the deformation is mainly along the formatted sliding zone or the local area of sliding mass rather than along the weak intercalated layer.
According to the drilled holes logs, the slide mass was composed of a soil-rock mixture with the thickness of 7.75∼38.50m.It comprised various dimensions of subangular siltstone gravels together with silty clay and clay; the gravel content was approximately 10∼70% (Figure 2(a)).
ere was no unified and continuous layered boundary among the geomaterials.As revealed, the slide zone was composed of subrounded or subangular gravels mixed with silty clay, with an average thickness of 0.7 m. e gravel with the dimension of 1∼6 cm accounted for a proportion of 20∼30% in the sliding-zone soil (Figure 2(b)).
According to the drilled holes logs and engineering geological investigation, the representative engineering geological cross section can be plotted, as shown in Figure 3.

Interpretation of SNMR Data.
ere were 15 SNMR survey points arranged in the landslide site in total, which are distributed in 4 survey lines, as shown in Figure 1(c).e SNMR device NUMIS Poly from IRIS Instruments was adopted in this study.e data were measured and interpreted in September 2013 for the spatial variability study and landslide deformation behavior analysis.e built-in procedures of the NUMIS system were applied to address the first-hand surveyed data.After backward extension and noise filtering of the signal, the data can be interpreted using a classic inversion method of least squares based on Tikhonov regularization.
e measured and estimated data for each vertical survey line were not continuous.erefore, the estimated data with a presupposed interval of 2 m were set for data extraction.
e points of discontinuity were taken into account as well at the same time.As has been measured by the comprehensive field/laboratory tests, the double-ring infiltrometer tests indicated that the K s values of the landslide surface were within the range of 4.2 × 10 −5 ∼5.0 × 10 −4 cm/s; the borehole water injection tests indicated that the K s values of the slide mass were within the range of 2.1 × 10 −6 ∼1.4 × 10 −4 cm/s.e laboratory tests indicated that the K s values of silty clay were within 3.05 × 10 −5 ∼7.12 × 10 −5 cm/s, while for clay, the result was 5.77 × 10 −6 cm/s.According to the aforementioned field measurements, the constant b ′ was calibrated as 4.32 herein.Figure 4 presents the comparison of the lithology indicated by drilled cores together with the surveyed T * 1 values as well as the estimated ∅ and K s values.Notably, we discarded the measured data of the bedrock and slide zone based on the consideration that the research object was the slide mass of the Baishuihe landslide.Advances in Civil Engineering of the estimated data, as shown in Table 1.Notably, the estimated K s results for the overall landslide site and major cross section had the range between 3.00 × 10 −6 ∼7.80 × 10 −3 cm/s and 4.70 × 10 −6 ∼7.00 × 10 −3 cm/s, which reaches about 3 orders of magnitude.e magnitude of the mean and median values corresponds with the field measurements.e median values were smaller than the mean value for the two analyzed research locations, which reflected the left skewness of the analyzed data.e median value is believed to be more representative generally because it rises superior to the effect of the extreme values compared to the mean value.e CV, a dimensionless parameter of the standard deviation divided by the mean value, representing the degree of variation, was 46% and 53% for the overall site and major cross section, respectively.ese results implied moderate variability for both the overall site and the major cross section.

Advances in Civil Engineering
With regard to the probability distribution, the frequency histograms of the estimated parameters are shown in Figure 5.As has been exhibited, both of the histograms for K s for the overall site and the major cross section presented the pattern of left skewness, indicating that the data had a nonnormal distribution.A logarithmic transformation (natural base) was performed on the estimated data for further statistical analysis.As shown in Table 1 and Figure 5, the log-transformed values were normally distributed.Lilliefors tests also indicated that the K s values for both the overall site and major profile were distributed log normally.

Spatial Variability Analysis.
Geostatistics is a useful tool for analyzing spatial variability and interpolating between point observations [35].e variogram which describes the spatial variability of K s is a function of the sampling distance under the intrinsic hypothesis of the variable, as shown in the following equation:  ) are employed to t the calculated experimental variograms.e best-tted theoretical models indicated by the greatest coe cient of correlation (R) and the least residual sum of squares (RSS) can be undertaken to represent the experimental variograms.
In this study, we investigated the omnidirectional spatial variable patterns for ln K s values of the overall site and major pro le because of the insu cient data pairs along some characteristic directions.e scale e ect caused by the di erent data intervals along vertical and horizontal directions was also taken into consideration by calculating di erent experimental variograms, which was characterized by small-scale variograms and large-scale variograms, respectively.
e calculated experimental variograms and the modeled theoretical variograms tted by the exponential model, Frequency (%) e corresponding fitting parameters for the 3 models are listed in Table 2.It is noticeable that a was the fitting parameter which equaled to A for the spherical model, while it equaled to 3a and � 3 √ a for exponential and Gaussian models, respectively.A greater A represents a greater range of variable regular changes and greater spatial correlation [37].
e greatest R and the least RSS indicated the best fit model for ln K s , i.e., K s .For the overall site study, the exponential model provided the best fit for the small-scale study because of the least RSS and the slight difference among the three models, while the spherical model provided the best fit for the large-scale study.For the major cross section, the spatial variation for a small-scale study was insignificant and the spherical model provided the best fit for a large-scale study.According to the fitting parameters, the (nested) variogram functions can be deduced for representing the variogram structure [34].
e variogram function for ln K s of the overall site and major profile can be deduced as follows.
Overall site: Major profile: Equations ( 5) and ( 6) indicated that the studied variables (ln K s values) were spatially correlated within A of the largescale study.
e A values reached 295.89 m and 65.56 m, respectively, for the overall site and major cross section, which suggested that the investigated spatial variation was remarkably dependent on the research direction.
e directional variogram structure can be investigated in the further studies in terms of more comprehensive measurements and estimations.Meanwhile, the effect of small-scale spatial variability contributed significantly to the nugget and sill, which represents the randomness derived from the microstructure and errors as well as the variable variation degree within the studied range.
e nugget-to-sill ratio (NSR) was employed to define the degree of spatial dependency.e greater the NSR value, the greater the variation degree caused by randomness.A low NSR value (<25%) indicates a strong spatial dependence of the variables under the strong impact of structural factors.A higher NSR value (25%∼75%) indicates a moderately spatial variable dependence.If it is higher than 75%, the variable is of weak spatial dependence, and randomness plays an important role [11,38].As shown in Table 2, all NSR values (except for the NSR of ln K s for the overall site with a small scale) were lower than 25%, suggesting that the spatial variation of K s is mainly dependent on the deposit structure.
e K s for the overall site with a small scale depended on structural factors and randomness.

Deformation Behavior Analysis Based on
Spatial Variable K s

Analyzing Method and Calculation Model.
A 2D deformation behavior analysis of the major profile associated with the spatial variable K s value under the combined effect of reservoir water fluctuation and rainfall has been carried out in this study.e commercial finite-element code SEEP/W (by Geoslope International Ltd.) [39], which is capable of solving finite-element equations for the steadystate flow governed by Darcy's law and for transient flow under the hypothesis of Terzaghi-Rendulic uncoupled consolidation, was used herein.
An ordinary Kriging estimation method was then applied to the orthogonal structural grids (4 m × 4 m) subjected to the geometry of the slide mass on the basis of conditional spatial data (SNMR-measured K s ) and a covariance matrix derived from the variogram structure analyzed in the aforementioned section.e method is believed to provide the best unbiased linear estimation for the unknown points [36,37] using a system equation: where K is the covariance matrix; I is a column vector of all points with a length of p; β is the vector of weights of the known data with  p i�1 β i � 1; k is the vector of covariance between the estimated points and known points; and λ is a Lagrange multiplier.As shown in Figure 7, the variation range of the K s values transformed from the interpolated ln K s values of the research domain was 1.37 × 10 −6 ∼1.11 × 10 −2 cm/s, which coincided with the laboratory and field measurements.
e calculation model was generated in SEEP/W.As shown in Figure 8, the mesh element of the slide mass had the same specification with the estimation grid.ere were 1968 nodes in total and 1970 elements available in this calculation model.e K s value of each element of the slide mass can be assigned individually by rewriting the XML file 8 Advances in Civil Engineering of SEEP/W.e front of the slope related to the reservoir water level, while the surface above the water level was assigned as a rainfall boundary with variable precipitation.
e back edge was assigned a constant water head of 240 m, while the bottom was assigned as a zero-ux boundary.e SIGMA/W code was subsequently utilized for analyzing the  Advances in Civil Engineering deformation behavior associated with the seepage eld results.
According to the laboratory experimental measurements and engineering geologic analogy method, the adopted calculation parameters are shown in Table 3. Apart from K s values of the slide mass, the other parameters for the landslide were identi ed as a single variable.e volumetric water content function for the slide zone was estimated on the basis of laboratory experimental results, of which α 0.029 kPa −1 , n 1.24, θ r 0.103, and θ s 0.315, while the volumetric water content function for the slide mass was estimated on the basis of the grain size distribution and the estimated K s according to the modi ed Kovacs method [40].
ree inclinometer holes were arranged in this major pro le corresponding to ZK01, ZK02, and ZK03, as shown in Figure 8.However, two of the three inclinometer tubes had been destroyed because of undetermined reasons.erefore, ZK01 was the remaining inclinometer hole for in-depth displacement monitoring.
e actual daily rainfall precipitation and reservoir water level of the exact monitoring period, from December 28, 2013, to October 20, 2014, were selected as the calculation condition for seepage and deformation behavior analysis.

Results and Comparison.
e simulated seepage and deformation behavior for 4 sections (as shown in Figure 9) were investigated in this study, including the three sections (S1, S2, and S3) corresponding to the drilled holes with inclinometers of ZK03, ZK02, and ZK01, respectively, as well as S4 representing the unavailable bottom section beneath the water level.e horizontal cumulative displacement time history of the surfaces of the 4 studied sections are presented in Figure 10 to compare with the outer condition of rainfall and uctuating water level.A signi cant deformation was detected at the bottom of the landslide (S4) at the very beginning under the condition of reduction of the reservoir water level.However, the displacement sustained at approximately 200 mm from January to July 2014 until a heavy rainstorm occurred in August 7 and facilitated the deformation under the combinational e ect of reservoir water storage.e time history of S4 indicates that the deformation of the front of the landslide is sensitive to both uctuating water level and rainfall.However, the upper sections of the landslide (S1-S3) produced negligible deformation during the entire simulation period relatively, which indicates that the deformation of the landslide is progressive.e onset of the deformation of S1-S3 was initially facilitated by an expedited drawdown of the water level and relatively intense rainfall.
e step-like displacement was re-onset by the intensive rainfall in July and August of that year.e deformation for S2 and S3 increased drastically after a remarkable heavy rainfall on August 7, 2014, and lasted for two months after the remarkable rainfall was associated with the increasing water level.Miao et al. also proposed that the displacement of the Baishuihe landslide changed stepwise at the rainfall-concentrated months from May to October of every year through a multiyear displacement analysis [33].e multiyear displacement monitoring data also indicated that the landslide has cumulated a spectacular displacement (approximately 1∼2 m) at the intermediate section by stages since the first impoundment of the ree Gorges Reservoir in 2007, but it remains stable [30][31][32][33] because of the toe resistance result from the gentle slope of the interface of the sliding mass and the bedrock.
During the simulation period, five representative points in time were selected to evaluate the seepage and deformation processes, as illustrated in Table 4. e simulated pore water pressure and displacement profile of the investigated points in time along the 4 sections were extracted from the numerical simulation, as presented in Figure 11.As indicated by Figures 11(a)-11(c), the stable depth of the zero-pore water pressure implying the embedded water table during the entire simulation period suggests that the water level and rainfall have little impact on the underlying water table at the upper and intermediate sections.However, the water table at the toe section varies remarkably in close relation to the fluctuating water level and rainfall (Figure 11(d)).e pore water pressure profiles above the bedrock for S1, S2, and S3 present a convoluted pattern with a negative value, indicating that the suction is greatly dependent on the daily rainfall at the upper section of the landslide.
As suction dissipated in the slide mass above the bedrock because of the rainfall, the kinematics of the landslide deformation changed with progressive movement in the geomaterials.e deformation of the top section (S1, as in Figure 11(e)) presents the pattern of greatest displacement at the surface dissipating along with the depth, without drastically changing at the sliding zone, which is implicative of a superficial deformation.e simulated displacement of the surface of S1 reached the peak value (58.40 mm) on October 7 2014 after a heavy storm, and the subsequent deformation was negligible.e deformation for S2 and S3 followed a translational mechanism, represented by relative uniform displacements above the sliding zone (Figures 11(f ) and 11(g)).
e in-field-monitored depth-displacement profile for S3 (ZK01) is also exhibited in Figure 11(g).
e monitored displacement deflected at the depth of 28 m beneath the ground surface, corresponding to the depth of the sliding zone.Such commensurate patterns with the simulated displacement curve demonstrate the reliability of the numerical results.Relatively, the displacement of the toe (S4) exhibits a distinct pattern that is tortuous with a great rotation at the depth of 8 m for the entire investigated period, which is implicative of a secondary sliding surface located in this depth at the toe.e deformation of S4 was significant from the very beginning and maintained its stability until October 2014, which is implicative of the great effect of the heavy rainfall and fluctuating water level.
According to the field investigation on October 20 2014, the Baishuihe landslide has produced remarkable crack in the middle platform (Figure 12(a)) and shallow collapse of the back edge (Figure 12(b)).e inclinometer PVC tubes located in ZK02 and ZK03 were discovered to have been destroyed at the depth of less than 1 m. e distributed optical fiber sensor arranged at the front of the landslide was discovered to be pulled apart as well.However, the exact time of occurrence of the great deformation is not verifiable.e newly outcropped large crack and collapse suggested that the deformation was not occurring for a long time, no more than 2 months.It is noticeable that the maximum simulated displacement for ZK02 and ZK03 (59.84 mm and 31.39 mm) was not large enough for the complete damage of the inclinometer pipes in general.e reason lies on the shallow deformation along the weaker zone (not sliding zone) due to the complex effects of vegetation cover, development of tension cracks, erosion, etc. e shallow deformation is also difficult to be captured during the numerical simulation.
All of the field observations and simulated results provided insights into the complex deformation behavior of the

Advances in Civil Engineering
Baishuihe landslide.
e deformation behavior is progressive, occurring initially at the front of the landslide and progressing backwards.
e deformation is greatly influenced by rainfall and reservoir water fluctuation.When the reservoir water level drops, a hydraulic gradient is produced due to the quick decrease of the water level compared to that of the groundwater table and it exerts an adverse effect on the landslide stability.When the reservoir stores water, the groundwater table level lifts gradually, thus increasing the weight of the geomaterial beneath the water table and softening the material.e pore water pressure in the slide mass and the sliding zone increases simultaneously; thus, the effective stress of the geomaterial beneath the water table decreases and produces deformation.In the case of heavy rain, the infiltration of the rain decreases the suction of the slide mass and thus decreases the effective stress.Besides, the heavy rainfall produces surface runoff and causes the erosion of the landslide surface, which exacerbates the deformation of the surface.Consequently, the reservoir water impoundment and heavy rainfall are the most important factors that affect the Baishuihe landslide deformation.12 Advances in Civil Engineering the landslide associated with the in-eld observations.e conclusions from the study are as follows:

Conclusions
(1) e SNMR technology provides a strategy for obtaining a series of relatively reliable spatial measurements of hydrological parameters of landslide elds using a noninvasive method.e estimation should be calibrated according to the eld tests measurements.
(2) e estimated K s values for the overall landslide mass site are log-normally distributed within a wide range between 3.00 × 10 −6 cm/s and 7.80 × 10 −3 cm/s and are associated with the CV of 34%.e variogram    is study provides an unconventional framework for studying the heterogeneous geomaterial and contributes to a better understanding of the spatial variability of the hydraulic property of accumulation reservoir landslides at a field scale.Advances in Civil Engineering

Figure 5 :
Figure 5: Histograms of the estimated K s values: (a) K s for the overall site; (b) K s for the major cross section; (c) ln K s for the overall site; (d) ln K s for the major cross section.

Figure 6 :
Figure 6: Omnidirectional variograms of the estimated ln K s values: (a) the small-scale variogram for ln K s for the overall site; (b) the largescale variogram for ln K s for the overall site; (c) the large-scale variogram for ln K s for the major cross section.

Figure 7 :Figure 8 :
Figure 7: Estimated K s value of the slide mass based on the ordinary Kriging interpolation method.

e
presented case study described the spatial variability of K s values of a representative accumulation reservoir landslide in TGRA, China, estimated by the SNMR survey work arranged in the landslide field.e finite-element seepagestress simulation results based on the spatially variable K s values provided insights into the deformation behavior of

Figure 11 :
Figure 11: Simulated pore water pressure and displacement over time along sections of the Baishuihe landslide.
Advances in Civil Engineeringanalysis considering the scale effect caused by the different data intervals along different directions shows that the nested variogram has A of 298.89 m and 65.56 m, respectively, for the overall site and major cross-sectional analysis.e spatial variation of K s is mainly dependent on the structure of the deposit.(3) e displacement time history as well as the pore water pressure and displacement of the studied sections based on the finite-element seepage-stress analysis of a Kriging-interpolated spatial variable K s value calculation model provide important insights into the deformation behavior of the landslide associated with the in-field observations and monitoring data.All of the results indicate that the Baishuihe landslide is a progressive landslide with significant secondary sliding at the front section and superficial deformation at the back section.e available information indicates the main factor influencing the landslide deformation is rainfall and reservoir water fluctuation.

Figure 12 :
Figure 12: Large deformation observations in the Baishuihe landslide site.(a) Great fissure at the intermediate platform.(b) Small collapse at the top section of the landslide.

Table 1 :
Descriptive statistics for the K s value of the slide mass of the Baishuihe landslide.K s , saturated hydraulic conductivity (cm/s); SD, standard deviation; CV, coe cient of variation (%). Note:

Table 2 :
Parameters of tting theoretical variogram models for the ln K s value.
Note: NSR, nugget-to-sill ratio; R, coe cient of correlation; RSS, residual sum of squares; the best-t theoretical model is marked in bold.

Table 3 :
Parameters of the calculation model.

Table 4 :
Representative points in time and description.
Note: √ means having inclinometer monitoring data for comparison.