Research Article  Open Access
Shen Zhang, Junwei Ma, Huiming Tang, "Estimation of Risk Thresholds for a Landslide in the Three Gorges Reservoir Based on a KDECopulaVaR Approach", Geofluids, vol. 2020, Article ID 8030264, 13 pages, 2020. https://doi.org/10.1155/2020/8030264
Estimation of Risk Thresholds for a Landslide in the Three Gorges Reservoir Based on a KDECopulaVaR Approach
Abstract
The Three Gorges Reservoir area, one of the most landslideprone areas in China, is characterized by widely distributed deepseated landslides exhibiting creep deformation due to rainfall and reservoir fluctuation. Thresholds, which are a key component for a reliable landslide early warning system, are still lacking for the prediction of movements of deepseated reservoir landslides with creep deformation because information about reservoir fluctuation indicators is lacking, uncertainty is ignored, and binary output is provided. The risk threshold, defined as the tolerance criteria for risks that will lead to action, is an effective measure of the degree of uncertainty. In the present study, a hybrid approach utilizing kernel density estimation, a copula function, and the value at risk is proposed for the estimation of the risk threshold for the Baishuihe Landslide, a typical deepseated landslide in the Three Gorges Reservoir area. Historical observations over approximately 15 years including rainfall, reservoir fluctuation, and landslide velocity were used to extract the risk threshold. A threelevel risk threshold describing the minimum magnitudes of rainfall and reservoir fluctuation for changing the landslide movement state under three confidence levels was developed. A threelevel risk response procedure, including risk responses in yellow alert, orange alert, and red alert, is proposed for risk management. Given its successful application, the present approach can be used to estimate the risk threshold for deepseated landslides.
1. Introduction
The movement and failure of landslides can cause significant societal and economic losses. According to statistics, landslides have caused 10,338 deaths during the period of 20082017 [1]. Rainfall is considered to be the primary causal factor of landslide movement and failure [2, 3]. A threshold is defined as the minimum or maximum level of a quantity needed for a process to take place or a state to change [4], and this concept is a key component for landslide early warning systems. For rainfallinduced landslides, the rainfall threshold defines the lower bound of known hydrological conditions (e.g., rainfall intensity, rainfall duration, or soil moisture) that has resulted in landslides [5–7], and this concept has been widely used to forecast the possible occurrence of a landslide in a given study area. In the 1980s, a global rainfall threshold for shallow landslides and debris flows was proposed based on 73 cases [8]. Following this pioneering work, a series of rainfall thresholds ranging from global to regional to local scales [6, 7, 9–21] have been proposed based on landslide inventories and precipitation. The available rainfall thresholds can also be classified according to both the estimation approaches that are used (i.e., physically based approaches or empirical approaches) and the rainfall variables that express the threshold (e.g., rainfall duration, cumulative event rainfall, rainfall intensity, or antecedent rainfall) [22]. These rainfall thresholds have been widely applied due to the advantage of determining the timing for the occurrence of shallow landslides.
The Three Gorges Reservoir area, one of the most landslideprone areas in China, is characterized by deepseated landslides. The deepseated landslides there exhibit creep deformation which can be described as a series of steps consisting of rapid movements at certain times and suspended activities during other periods due to rainfall and reservoir fluctuation [23, 24]. And changes of landslide movement state from suspended activities to rapid movements pose a significant threat to lives and properties. For example, the movement of the Qianjiangping Landslide, which was triggered by the initial filling of the Three Gorges Reservoir, caused the deaths of 24 people and damage to 129 houses [25]. However, only a few efforts have attempted to estimate the thresholds in the Three Gorges Reservoir area [26–28]. The available thresholds are classifiers that provide a binary output (landslide or no landslide), which are more applicable for sudden and rapid rainfallinduced failures. Moreover, hydrological indicators (such as reservoir fluctuation) that have been acknowledged to be domain causal factors for the movement of deepseated landslides [29, 30] are not usually included in the available thresholds. Furthermore, little attention has been paid to the uncertainty of the models. In fact, a certain degree of uncertainty exists in the identification of rainfall thresholds [31], given that only the major variable, rainfall, is included, and less important variables are excluded. Poor quality observations or insufficient data sets may also increase the involved uncertainty.
Therefore, it is of great value to define thresholds to describe the minimum magnitudes of rainfall and reservoir fluctuation for changing the landslide movement state and to thus implement an effective tool for landslide risk management in the Three Gorges Reservoir area. The risk threshold, defined as the tolerance criteria for risks that will lead to action, is an effective measure of the degree of uncertainty [32, 33]. It has been widely applied as a popular risk management measure in the fields of finance [34], hydrology [35, 36], and energy storage operation [37] due to the advantage of conceptual simplicity.
The purpose of the present study is to estimate the risk threshold for deepseated reservoir landslides based on kernel density estimation (KDE), a copula function, and the value at risk (VaR). The Baishuihe Landslide, a typical deepseated landslide in the Three Gorges Reservoir, is selected as a case study. Based on approximately fifteen years of observations, a threelevel threshold under three confidence levels is developed for risk management.
2. Materials and Methods
2.1. Features of the Baishuihe Landslide
The Baishuihe Landslide (latitude and longitude: N31^{°}0134 and E110^{°}3209, Figure 1(a)) is located on the south bank of the Yangtze River approximately 56 km upstream from the Three Gorges Dam. The main sliding direction is 20°. The toe elevation of the landslide is approximately 130 m, but the toe is submerged by the Yangtze River, while the crown varies from 350 to 400 m (Figures 1(b) and 1(c)). The landslide is approximately 600 m long and 700 m wide. The average thickness is 30 m, and the landslide volume is estimated to be 12.6 million cubic meters. The mean inclination of the landslide surface is 30°.
(a)
(b)
(c)
Based on the deformation characteristics, the Baishuihe Landslide can be divided into two blocks [38]: the warning zone and the relatively stable block (Figures 1(b) and 1(c)). The warning zone is located in the front part, with elevations ranging from 130 to 270 m. The warning zone is approximately 450 m long and 500 m wide. The planar area of the warning zone is 0.225 million square meters, and the volume is 6.75 million cubic meters (Figures 1(b) and 1(c)). Borehole analysis shows that the landslide materials are Quaternary deposits consisting of gravel clasts and silty clay. The underlying bedrock is sandstone and mudstone of the Triassic Shazhenxi Formation, with an average dip direction of 15° and a dip angle of 36°.
The Baishuihe Landslide is an ancient landslide. Since the 1990s, continuous deformation has been recorded and has caused damage to roads and retaining walls (Figure 2(a)). In August 25, 1993, a local slump occurred in the rear part, and fifteen residents had to be relocated. In June 2003, after the initial impoundment of the Three Gorges Reservoir, transverse cracks were noticed in the eastern part and reached 300 meters long. During August 2005 and August 2008, multiple slump events (Figure 2(b)), with a maximum volume of , occurred along the road at an elevation of 220 m [38, 39].
(a)
(b)
2.2. Input Data
With the aims of realtime monitoring and early warning, an integrated monitoring system consisting of highaccuracy double frequency Global Positioning Systems (GPSs), a Leica TCA2003 total station, inclinometers, water level gauges, earth pressure cells, PulsePre Pump Brillouin Optical Time Domain Analysis (PPPBOTDA) optical fiber, and an automatic weather station was installed in the landslide area. To ensure the quality of the data, multiagency monitoring was performed routinely by local residents, a local meteorological station, the National Field Observation and Research Station of Landslide in Three Gorges Reservoir Area of Yangtze River, and the China Three Gorges Corporation. More than fifteen years of regular monitoring data have been collected since June 2003. The selected observations from regular monitoring (Table 1) show that the displacement time history can be described as a series of steps consisting of rapid movements at certain times (June and July of each year) and suspended activities during other periods. And significant deformation with an average annual velocity up to 300 mm/year occurred in the warning zone; however, small deformations with an annual velocity of millimeters in magnitude were observed in the relatively stable block. At the end of the monitoring period, the observed displacements for monitoring points ZG118 and ZG93 located in the warning zone were 3148.4 and 3206.7 mm, respectively. The displacement direction of monitoring point ZG118 is 20°, which corresponds to the main sliding direction. Therefore, the observations from monitoring point ZG118 are selected for the estimation of the risk threshold. The observations of rainfall intensity, reservoir fluctuation, cumulative displacement, and velocity at ZG118 from June 2003 to December 2017 are shown in Figure 3. The available data reveal that the Baishuihe Landslide exhibits steplike deformation and that rapid movements occur in the rainy season (May to September) and then essentially cease during the dry season (October to April).

Figure 4 displays histograms and empirical cumulative distribution functions (CDFs) for rainfall, reservoir fluctuation, and landslide velocity. Rainfall, reservoir fluctuation, and landslide velocity clearly follow different distribution patterns. Therefore, a universal framework capable of handling different distribution patterns and being further extended to more complex variables is required.
2.3. Methodology
2.3.1. Risk Threshold and Value at Risk (VaR)
It is generally accepted that risk arises from uncertainty [40]. According to the definition from the Project Management Institute, a risk is an uncertain event or condition that, if it occurs, has a positive or a time cost, span, or quality, which implies an uncertainty in the identified events and conditions [32]. A risk threshold refers to measures along the level of uncertainty at which a stakeholder may have a specific interest. Below that risk threshold, the organization will accept the risk. Above that risk threshold, the organization will not tolerate the risk [32, 33]. In the present study, VaR, an effective risk management tool in the field of finance, is utilized to estimate the risk threshold.
In finance, VaR defines as a threshold value to measure the potential financial loss within a given time horizon and at a confidence level. It has been widely adopted by practitioners and regulators as the standard method of measurement of the financial loss. Define as the confidence level ranging from 0 to 1 and as the loss measured as a positive value. is the threshold at which the probability of a loss exceeding the threshold is equal to (schematically illustrated in Figure 5) and can be expressed as follows:
Take, for instance, a 99% confidence level (). is the loss such that the probability of experiencing a greater loss is less than 1%.
Various methods, including the variancecovariance method, historical simulations, and Monte Carlo simulations, are available for the estimation of the VaR [41]. However, these methods are based on the assumption of probability distributions, such as normal distributions, lognormal distributions, or any of a number of other distributions [42, 43]. Moreover, a large number of samplings are required in a Monte Carlo simulation. Additionally, when addressing practical applications with multiple factors, VaR estimation can become very difficult due to the complexity of joint distribution estimation. In the present study, the copula, an efficient mathematical tool capable of constructing multivariate distributions from their univariate marginal functions without any assumptions of normality and linear correlation, is utilized for VaR estimation.
2.3.2. Copula
The copula theory, first developed by Sklar in the 1960s, states that any multivariate distribution function can be written in terms of univariate marginal distributions and a copula function [44]. The flexibility of specifying the marginal distributions separately from the copula function that links these distributions to form the joint distribution has resulted in its wide application by researchers for multivariate analysis [45]. Taking the bivariate copula as an example, let be a twodimensional joint distributional function of random variables and with univariate margins and . Then, there exists a copula such that
Various copulas exist, including Gaussian, Clayton, and Frank. A more detailed description of the copula function is listed in Table 2.

The parameter of the copula can be estimated from the original data through the maximum likelihood method based on the loglikelihood function as follows:
In the present study, the Euclidean distance between the empirical copula and the optional copula is chosen to determine the optimal copula as follows: where is the data size and is the empirical copula. The empirical copula [46, 47] is defined as follows: where denotes the indicator function and represent order statistics from the sample .
The construction of bivariate or multivariate copulas usually starts by specifying marginal univariate distributions. In this study, nonparametric KDE is applied for marginal univariate distribution estimation.
2.3.3. Kernel Density Estimation
As one of the most wellknown nonparametric density estimators, KDE is capable of learning the shape of the density of a data set automatically without any distribution hypotheses. Let be a sample from a distribution with density . Nonparametric density estimation refers to the estimation of with as few assumptions as possible. Various approaches, including the nearest neighbor method, Parzen windows, and KDE, are available for nonparametric density estimation. However, the performances of the nearest neighbors and Parzen windows methods are generally poor [48]. The output of the Parzen windows method is usually characterized by discontinuous curves with steplike features. Fortunately, these drawbacks can easily be eliminated by adopting KDE.
The estimator of is given by where is the kernel and is the bandwidth. Technically, KDE smooths out each data sample into a bump, the shape of which is determined by the kernel. Then, KDE estimates the shape of the data set by summing over all these smoothed bumps.
By adoption of the most efficient Epanechnikov kernel [48], the estimator can be written as where is the indicator function, that is, for and for .
2.3.4. Risk Threshold Estimation Based on KDECopulaVaR
The overall framework of estimating the risk threshold based on a KDECopulaVaR approach mainly consists of three steps: (1) marginal distribution estimation by KDE, (2) joint distribution construction by the copula, and (3) risk threshold extraction by VaR.
First, the marginal distribution functions of rainfall intensity (), reservoir fluctuation (), and landslide velocity () were estimated by Epanechnikov KDE. The bandwidths of the Epanechnikov KDE for variables related to rainfall, reservoir fluctuation, and velocity were set to 0.05, 0.05, and 0.5, respectively. Figure 6 shows the empirical CDFs and estimated CDFs based on KDE for rainfall intensity (), reservoir fluctuation (), and landslide velocity (). The CDF estimations based on KDE match well with the empirical CDFs. This perfect performance is based on the nonparametric nature of KDE.
Second, joint probability distribution functions were constructed with a copula. The popular copula functions listed in Table 2, including the Gaussian, , Gumbel, Clayton, and Frank functions, were set as candidate copula functions. The Euclidean distance between the empirical copula and the optional copula was calculated to determine the optimal copula function. The Euclidean distances between the empirical copula and the optional copula of (, ) and (, ) are listed in Table 3. The Euclidean distances indicate that the Frank copula was the optimal copula function with a minimum Euclidean distance. Therefore, the Frank copula was selected to form a multivariate distribution function.
 
The bestfit copula is in italic font. 
Third, a threelevel risk threshold at three popular confidence levels (, 0.90, and 0.85) [49] was extracted from quantiles of the multivariate joint probability distribution functions.
3. Results
The joint Frank copula probability density functions (PDFs) and CDFs of (), (), and () are shown in Figures 7 and 8. Moreover, the correlation coefficients between rainfall intensity and landslide velocity and between reservoir fluctuation and landslide velocity are measured by two indexes including Kendall’s tau and Spearman’s rho and listed in Table 4. Kendall’s tau and Spearman’s rho have values ranging from 1 to +1. Zero indicates no correlation, and +1 and 1 indicate complete positive and negative correlations, respectively. The copula analysis shows a moderate correlation (Kendall’s , Spearman’s ) between rainfall and landslide velocity. The constructed joint PDF between the rainfall and landslide velocity (Figure 7) is characterized by a symmetric tail that indicates a moderate correlation between the lower tail (0,0) and the upper tail (1,1). That is, if the landslide area experienced a heavier rainfall event, the landslide mass underwent more deformation. The effect of rainfall infiltration on movement of the Baishuihe Landslide can be described as follows: under prolonged rainfalls, most of the rain infiltrates into the landslide mass with permeability coefficients up to according to in situ tests. The infiltrated water may cause the increases of landslide mass weight and porewater pressures and the reduction of shear strength along the sliding zone, thus triggering the landslide movements. Meanwhile, infiltrated rainfall may reduce the mechanical strength of the sliding zone, accelerating movements. A negative moderate correlation (Kendall’s , Spearman’s ) is found between the reservoir fluctuation and landslide velocity, suggesting that a higher landslide velocity was more likely to occur during a sharp drawdown of the reservoir. The movements of the Baishuihe Landslide during the sharp drawdown of the reservoir were caused by the reduction of support provided by the reservoir at the toe and high seepage pressures toward the toe.

The obtained results also indicate that the correlation of reservoir fluctuation and landslide velocity is stronger than that of rainfall and landslide velocity. This finding agrees well with the conclusion derived from the monitoring. The fast landslide movements mainly occur between May and June, when the water in the reservoir rapidly decreases to its lowest. For example, during June 2015, the reservoir rapidly decreased to 145 m at a daily fluctuation speed ranging from 0.6 to 0.99 m per day, and fast movement occurred with displacement increment up to 497 mm.
A threelevel threshold relating the rainfall intensity, reservoir fluctuation, and landslide velocity under three popular levels of confidence is shown in Table 5. To illustrate the essence of the obtained threshold models, three models are described in the following paragraphs.

The threshold at the 0.95 confidence level states that under prolonged rainfall (exceeding 314 mm per month) and rapid drawdown with a magnitude of 11 m per month, the landslide was certain to change into rapid movement state with a probability up to 95%. When those thresholds are fulfilled, the copulabased velocity for landslide movements will reach to 320 mm per month. The correlation relationships between rainfall intensity and rapid movement (Figure 9(a)) show that rapid landslide movements positively correlated to the rainfall intensity. The heavier the rainfall is, the more likely it is that significant landslide movements will occur. A negative correlation exists between the reservoir fluctuation and landslide velocity (Figure 9(b)). The sharper the reservoir drawdown is, the more likely it is for significant landslide movements to occur. Almost all rapid movement events fall below the obtained thresholds. These results correspond with the obtained values from the KDECopulaVaR approach.
(a)
(b)
The threshold at the 0.90 confidence level shows that under heavy rainfall (exceeding 206 mm per month) and water level drawdown at a moderate speed, namely, 8 m per month or 0.27 m per day, the landslide would be extremely likely to change into rapid movement state with probability up to 90%. When those thresholds are fulfilled, the copulabased velocity for landslide movements will reach to 51 mm per month. The reservoir fluctuation thresholds correspond to the work of He et al. [28], in which the reservoir drawdown criterion for landslide instability is 0.2 m per day.
The threshold at the 0.85 confidence level shows that under heavy rainfall (exceeding 174 mm per month) and water level drawdown at a speed of 6 m per month, the landslide would be highly likely to change into rapid movement state with probability up to 85%. When those thresholds are fulfilled, the copulabased velocity for landslide movements will reach to 42 mm per month.
The obtained thresholds also indicate that the copulabased velocities rapidly increase in exponential growth rate with confidence levels. This is because with higher confidence levels, people will be more risk tolerant and therefore willing to take more risk to achieve higher expected return.
The obtained thresholds can be valuable for risk mitigation and the implementation of an effective tool for landslide early warning systems. Decisions on risk response can be made based on the comparison between measured or forecasted values and the obtained thresholds (Figure 10). For example, the rainfall thresholds can be compared with the rainfall intensity recorded from the beginning of the rainfall event or the forecasted rainfall intensity. When the rainfall is close to or exceeds a preestablished rainfall threshold, the appropriate risk response can be raised. For the purpose of landslide early warning, a threelevel risk response procedure is proposed for risk mitigation (Figure 10 and Table 5). The threelevel risk response procedure mainly consists of the following steps: risk identification and risk quantification, comparison with the risk threshold, and response to the risk (Figure 10).
A red alert is raised, indicating a high potential risk of landslide failure, when the risk value exceeds the risk threshold at the 0.95 confidence level. Mitigation measures including 24 h continuous and comprehensive monitoring, general inspection, deployment of warning signs, closure of roads and river channels, and issuance of evacuation orders should be employed.
An orange alert is raised when the risk value exceeds the risk threshold at the 0.90 confidence level but is not higher than the threshold at the 0.95 confidence level. The following mitigation measures should be implemented: 24 h continuous and comprehensive monitoring and general inspection; development of emergency mitigation, preparedness, response, and evacuation plans; and a consultation meeting between experts and government decision makers.
A yellow alert is raised when the risk value exceeds the risk threshold at the 0.85 confidence level but is not higher than the threshold at the 0.90 confidence level. Implementation of the following measures is recommended to minimize the landslide risk: more frequent and comprehensive monitoring and issuing the information to expert and decisionmaker groups.
Otherwise, the risk events are withheld or ignored.
4. Discussions
In order to estimate early warning criteria of rainfall and reservoir fluctuation and quantify uncertainty associated with the estimation, a threelevel risk threshold under three confidence levels is extracted using a KDECopulaVaR approach. In fact, the hybrid approach utilizing the KDE, a copula function, and the VaR is a datadriven approach. The major advantages of the proposed approach consist of the following.
The modeling procedure is quite straightforward, and less detailed physical mechanism information involved in landslides is required, which will save both time and cost for threshold estimation. The process of data acquisition and processing is relatively convenient, which means the threshold model can be easily established. The hybrid datadriven approach provides satisfactory performance once the threshold model is well trained.
However, limitations relating to data size and representativeness still exist. Firstly, sufficiently large size of observations is needed for training models with satisfactory performance. Generally, the trained models only extrapolate the past pattern of landslide movement. Since the future is never exactly like the past, a mere extrapolation of past patterns cannot provide accurate predictions, whichever sophisticated model is used to do so. To ensure better model performance, one should consider a sample period that is sufficiently large, so that the landslide movement patterns can be modelled more appropriately. Another limitation concerns the representativeness. Natural systems including landslide may change rapidly due to changing material supply and topography, restoration of vegetation, and human actions. The conditions that triggered landslide movements in the past may not be representative for the future. Therefore, the trained threshold models should be retrained and updated frequently when changes are made or more data is available.
In spite of the limitations described above, we still believe that the proposed hybrid approach can be useful for end users. From the perspective of practice, our research provides an analytical procedure for end users, enabling them to minimize landslide risk, especially when the associated safety factors are already marginal. However, end users still face essential challenges. Two of these challenges related to uncertainty and confidence level are discussed as follows.
The first challenge is the attitude of end users towards uncertainty. For some, uncertainty is perceived as a curse, with risk arising from uncertainty. For others, uncertainty is perceived as a blessing that offers possible benefits if addressed properly. The ways in which end users interpret uncertainty also depend on their attitude. When users face uncertainty, they may choose to either downplay uncertainty for the sake of efficiency or prepare for the worstcase scenario because of highrisk aversion [50]. In fact, as Kahneman [50] has noted, an unbiased appreciation of uncertainty is a basis for rationality. However, this viewpoint is at odds with the will of the people and organizations. In general, expressions of uncertainty by researchers are typically perceived as an indication of weakness and vulnerability. Researchers who display more confidence will easily gain the trust of end users and the public and will take the place of researchers who acknowledge the full extent of their ignorance.
The second challenge is how to evaluate the degree of confidence level. This factor is particularly important for those situations in which estimations are established to aid decisionmaking. In fact, the level of confidence will have an impact on the application of the obtained results.
To address the abovementioned challenges in practical applications, we should begin by changing the discourse about uncertainty, moving from a lack of knowledge to a more realistic and nuanced view. To this end, uncertainty could be treated as confidence and additional knowledge instead of as a lack of knowledge. Moreover, making uncertainty explicit and transparent should also be valued as a virtue similar to honesty, humility, and trust.
Based on experience and the expertjudgment framework within the Intergovernmental Panel on Climate Change (IPCC) [51], the following process (Figure 11) is provided to address the second challenge, evaluating the degree of confidence. Evidence and agreement are two metrics for evaluating the degree of confidence. The first step in the process is the identification of existing evidence (see step 1 in Figure 11). Types of evidence include experimental results, mechanistic understanding, observations, and models. The next step in the process is to evaluate evidence and agreement (see step 2 in Figure 11), particularly the type, amount, quality, and consistency of evidence and the degree of agreement. The degree of agreement is a measure of the consensus across the scientific community on a given topic. Qualitative language is used to characterize the amount of evidence and the degree of agreement among author teams. Evidence ranges from limited to robust, and the consistency of evidence ranges from low to high. In general, evidence is most robust when there are multiple independent and consistent sources of highquality evidence. The next step in the process is to evaluate the level of confidence based on evidence and agreement by integrating the evaluation of evidence and agreement into one metric (see step 3 in Figure 11). Increasing confidence is correlated with an increasing level of evidence and degree of agreement. The obtained confidence is expressed qualitatively using five qualities: very low (less than 1 in 10 chance of being subsequently corrected), low (approximately 2 in 10 chance of being corrected), medium (approximately 5 in 10 chance of being corrected), high (approximately 8 in 10 chance of being corrected), and very high (at least 9 in 10 chance of being corrected).
5. Conclusions
Due to rainfall and reservoir fluctuation, the Three Gorges Reservoir area is highly affected by the creep deformation of deepseated landslides. However, few attempts have been made to estimate early warning criteria that describe the magnitudes of rainfall and reservoir fluctuation for changing the landslide movement state. The risk threshold is a measurement along the level of uncertainty at which a stakeholder may have a specific interest, below which the organization will accept the risk, and above which the organization will not tolerate the risk. In the present study, a hybrid approach utilizing the KDE, a copula function, and the VaR is proposed to estimate the risk threshold describing the minimum magnitudes of rainfall and reservoir fluctuation for changing the landslide movement state for the Baishuihe Landslide, a typical deepseated landslide in the Three Gorges Reservoir. The following conclusions were obtained from this study: (1)A threelevel risk threshold describing the minimum magnitudes of rainfall and reservoir fluctuation for rapid movement under three confidence levels is extracted. The thresholds at the 0.95 confidence level depict a minimum monthly rainfall of 314 mm and a monthly reservoir decrease of 11 m. The thresholds at the 0.90 confidence level result in a minimum monthly rainfall of 206 mm and a monthly reservoir decrease of 8 m. The thresholds at the 0.85 confidence level yield a minimum monthly rainfall of 174 mm and a monthly reservoir decrease of 6 m(2)A threelevel risk response procedure, including risk responses in yellow alert, orange alert, and red alert, is proposed for risk management(3)Given the successful use of the KDECopulaVaR approach for estimating the risk threshold for the Baishuihe Landslide, this methodology is likely to be useful to estimate the thresholds for deepseated reservoir landslides. More specifically, the threelevel thresholds found in this paper relating rainfall and reservoir fluctuation and the threelevel risk response procedure could be applicable to other deepseated landslides in the Three Gorges Reservoir
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This study was financially supported by the National Key R&D Program of China (Grant No. 2018YFC1507200), the National Natural Science Foundation of China (Grant No. 41702328), and the Hubei Provincial Natural Science Foundation of China (Grant No. 2019CFB585). All support is gratefully acknowledged.
References
 D. Tiranti and R. Cremonini, “Editorial: landslide hazard in a changing environment,” Frontiers in Earth Science, vol. 7, no. 3, 2019. View at: Publisher Site  Google Scholar
 A. Dikshit and N. Satyam, “Probabilistic rainfall thresholds in Chibo, India: estimation and validation using monitoring system,” Journal of Mountain Science, vol. 16, no. 4, pp. 870–883, 2019. View at: Publisher Site  Google Scholar
 S. L. Gariano, R. Sarkar, A. Dikshit et al., “Automatic calculation of rainfall thresholds for landslide occurrence in Chukha Dzongkhag, Bhutan,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 6, pp. 4325–4332, 2019. View at: Publisher Site  Google Scholar
 I. D. White, D. N. Mottershead, and S. J. Harrison, Environmental Systems, Chapman & Hall, London, 2nd edition, 1996.
 R. K. Dahal and S. Hasegawa, “Representative rainfall thresholds for landslides in the Nepal Himalaya,” Geomorphology, vol. 100, no. 34, pp. 429–443, 2008. View at: Publisher Site  Google Scholar
 S. Segoni, L. Piciullo, and S. L. Gariano, “A review of the recent literature on rainfall thresholds for landslide occurrence,” Landslides, vol. 15, no. 8, pp. 1483–1501, 2018. View at: Publisher Site  Google Scholar
 A. Sengupta, S. Gupta, and K. Anbarasu, “Rainfall thresholds for the initiation of landslide at Lanta Khola in North Sikkim, India,” Natural Hazards, vol. 52, no. 1, pp. 31–42, 2010. View at: Publisher Site  Google Scholar
 N. Caine, “The rainfall intensity  duration control of shallow landslides and debris flows,” Geografiska Annaler: Series A, Physical Geography, vol. 62, no. 12, pp. 23–27, 1980. View at: Google Scholar
 F. Guzzetti, S. Peruccacci, M. Rossi, and C. P. Stark, “The rainfall intensity–duration control of shallow landslides and debris flows: an update,” Landslides, vol. 5, no. 1, pp. 3–17, 2008. View at: Publisher Site  Google Scholar
 F. Guzzetti, S. Peruccacci, M. Rossi, and C. P. Stark, “Rainfall thresholds for the initiation of landslides in central and southern Europe,” Meteorology and Atmospheric Physics, vol. 98, no. 34, pp. 239–267, 2007. View at: Publisher Site  Google Scholar
 H. Yang and R. F. Adler, “Predicting global landslide spatiotemporal distribution: integrating landslide susceptibility zoning techniques and realtime satellite rainfall estimates,” International Journal of Sediment Research, vol. 23, no. 3, pp. 249–257, 2008. View at: Publisher Site  Google Scholar
 P. Reichenbach, M. Cardinali, P. De Vita, and F. Guzzetti, “Regional hydrological thresholds for landslides and floods in the Tiber River Basin (Central Italy),” Environmental Geology, vol. 35, no. 23, pp. 146–159, 1998. View at: Publisher Site  Google Scholar
 G. Martelloni, S. Segoni, R. Fanti, and F. Catani, “Rainfall thresholds for the forecasting of landslide occurrence at regional scale,” Landslides, vol. 9, no. 4, pp. 485–495, 2012. View at: Publisher Site  Google Scholar
 S. Segoni, A. Rosi, G. Rossi, F. Catani, and N. Casagli, “Analysing the relationship between rainfalls and landslides to define a mosaic of triggering thresholds for regionalscale warning systems,” Natural Hazards and Earth System Sciences, vol. 14, no. 9, pp. 2637–2648, 2014. View at: Publisher Site  Google Scholar
 S. Lainas, N. Sabatakakis, and G. Koukis, “Rainfall thresholds for possible landslide initiation in wildfireaffected areas of western Greece,” Bulletin of Engineering Geology and the Environment, vol. 75, no. 3, pp. 883–896, 2016. View at: Publisher Site  Google Scholar
 G. Wieczorek, “Effect of rainfall intensity and duration on debris flows in central Santa Cruz Mountains, California,” Debris Flows/Avalanches: Processes, Recognition and Mitigation, Geological Society of America, Reviews in Engineering Geology, pp. 93–104, 1987. View at: Google Scholar
 M. L. Lee, K. Y. Ng, Y. F. Huang, and W. C. Li, “Rainfallinduced landslides in Hulu Kelang area, Malaysia,” Natural Hazards, vol. 70, no. 1, pp. 353–375, 2014. View at: Publisher Site  Google Scholar
 M.C. Chung, C.H. Tan, and C.H. Chen, “Local rainfall thresholds for forecasting landslide occurrence: Taipingshan landslide triggered by Typhoon Saola,” Landslides, vol. 14, no. 1, pp. 19–33, 2017. View at: Publisher Site  Google Scholar
 M. I. Sameen, R. Sarkar, B. Pradhan, D. Drukpa, A. M. Alamri, and H.J. Park, “Landslide spatial modelling using unsupervised factor optimisation and regularised greedy forests,” Computers & Geosciences, vol. 134, p. 104336, 2020. View at: Publisher Site  Google Scholar
 A. Dikshit, R. Sarkar, B. Pradhan, R. Jena, D. Drukpa, and M. A. Alamri, “Temporal probability assessment and its use in landslide susceptibility mapping for Eastern Bhutan,” Water, vol. 12, no. 1, p. 267, 2020. View at: Publisher Site  Google Scholar
 A. Dikshit, R. Sarkar, B. Pradhan, S. Acharya, and K. Dorji, “Estimating rainfall thresholds for landslide occurrence in the Bhutan Himalayas,” Water, vol. 11, no. 8, p. 1616, 2019. View at: Publisher Site  Google Scholar
 Y. Galanti, M. Barsanti, A. Cevasco, G. D. Avanzi, and R. Giannecchini, “Comparison of statistical methods and multitime validation for the determination of the shallow landslide rainfall thresholds,” Landslides, vol. 15, no. 5, pp. 937–952, 2018. View at: Publisher Site  Google Scholar
 B. Chai, K. Yin, J. Du, and L. Xiao, “Correlation between incompetent beds and slope deformation at Badong town in the Three Gorges reservoir, China,” Environmental Earth Sciences, vol. 69, no. 1, pp. 209–223, 2013. View at: Publisher Site  Google Scholar
 J. W. Ma, H. M. Tang, X. Liu, X. L. Hu, M. J. Sun, and Y. J. Song, “Establishment of a deformation forecasting model for a steplike landslide based on decision tree C5.0 and twostep cluster algorithms: a case study in the Three Gorges Reservoir area, China,” Landslides, vol. 14, no. 3, pp. 1275–1281, 2017. View at: Publisher Site  Google Scholar
 F. W. Wang, Y. M. Zhang, Z. T. Huo, T. Matsumoto, and B. L. Huang, “The July 14, 2003 Qianjiangping landslide, Three Gorges Reservoir, China,” Landslides, vol. 1, no. 2, pp. 157–162, 2004. View at: Google Scholar
 H. K. Chen, L. Wei, and L. Tan, “Review of research on empirical rainfall threshold of rainfallinduced landslide,” Journal of Chongqing Jiaotong University (Natural Science), vol. 5, 2012. View at: Google Scholar
 J. Chen, Z. F. Yang, and X. Li, “Relationship between landslide probability and rainfall in Three Gorges Reservoir area,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 17, pp. 3052–3056, 2005. View at: Google Scholar
 K. Q. He, X. R. Li, X. Q. Yan, and D. Guo, “The landslides in the Three Gorges Reservoir Region, China and the effects of water storage and rain on their stability,” Environmental Geology, vol. 55, no. 1, pp. 55–63, 2008. View at: Google Scholar
 J. W. Ma, H. M. Tang, X. L. Hu et al., “Identification of causal factors for the Majiagou landslide using modern data mining methods,” Landslides, vol. 14, no. 1, pp. 311–322, 2017. View at: Publisher Site  Google Scholar
 Z.l. Wei, Q. Lü, H.y. Sun, and Y.q. Shang, “Estimating the rainfall threshold of a deepseated landslide by integrating models for predicting the groundwater level and stability analysis of the slope,” Engineering Geology, vol. 253, pp. 14–26, 2019. View at: Publisher Site  Google Scholar
 E. I. Nikolopoulos, S. Crema, L. Marchi, F. Marra, F. Guzzetti, and M. Borga, “Impact of uncertainty in rainfall estimation on the identification of rainfall thresholds for debris flow occurrence,” Geomorphology, vol. 221, pp. 286–297, 2014. View at: Publisher Site  Google Scholar
 K. H. Rose, A Guide to the Project Management Body of Knowledge (PMBOK® Guide)Fifth Edition, Pennsylvania, USA, Newtown Square, 2015.
 S.G. Kim, J.J. Kim, and K.R. Kim, “A risk threshold calculation methodology for the construction projects applying value at risk,” Construction Research Congress 2005, pp. 1–9, 2005. View at: Google Scholar
 J. Cotter, “Downside risk for European equity markets,” Applied Financial Economics, vol. 14, no. 10, pp. 707–716, 2004. View at: Publisher Site  Google Scholar
 J. Piantadosi, A. V. Metcalfe, and P. G. Howlett, “Stochastic dynamic programming (SDP) with a conditional valueatrisk (CVaR) criterion for management of stormwater,” Journal of Hydrology, vol. 348, no. 34, pp. 320–329, 2008. View at: Publisher Site  Google Scholar
 L. Wang, B. Wang, P. Zhang, M. Liu, and C. Li, “Study on optimization of the shortterm operation of cascade hydropower stations by considering output error,” Journal of Hydrology, vol. 549, pp. 326–339, 2017. View at: Publisher Site  Google Scholar
 S. Moazeni, W. B. Powell, and A. H. Hajimiragha, “Meanconditional valueatrisk optimal energy storage operation in the presence of transaction costs,” IEEE Transactions on Power Systems, vol. 30, no. 3, pp. 1222–1232, 2015. View at: Publisher Site  Google Scholar
 F. S. Miao, Y. P. Wu, Y. H. Xie, and Y. N. Li, “Prediction of landslide displacement with steplike behavior based on multialgorithm optimization and a support vector regression model,” Landslides, vol. 15, no. 3, pp. 475–488, 2017. View at: Google Scholar
 X. L. Yi, H. M. Tang, Y. P. Wu, Y. Ge, X. Fan, and S. Zhang, “Application of the PPPBOTDA distributed optical fiber sensor technology in the monitoring of the Baishuihe landslide,” Chinese Journal of Rock Mechanics and Engineering, vol. 35, no. S1, pp. 3084–3091, 2016. View at: Google Scholar
 R. Olsson, “In search of opportunity management: is the risk management process enough?” International Journal of Project Management, vol. 25, no. 8, pp. 745–752, 2007. View at: Publisher Site  Google Scholar
 M. G. Sampid and H. M. Hasim, “Estimating valueatrisk using a multivariate copulabased volatility model: evidence from european banks,” International Economics, vol. 156, pp. 175–192, 2018. View at: Google Scholar
 D. L. Olson and D. Wu, “The impact of distribution on valueatrisk measures,” Mathematical and Computer Modelling, vol. 58, no. 910, pp. 1670–1676, 2013. View at: Publisher Site  Google Scholar
 P. J. G. Vlaar, “Value at risk models for Dutch bond portfolios,” Journal of Banking & Finance, vol. 24, no. 7, pp. 1131–1154, 2000. View at: Publisher Site  Google Scholar
 A. Sklar, “Fonctions de repartition an dimensions et leurs marges,” Publications de l'Institut de statistique de l'Université de Paris, vol. 8, pp. 229–231, 1959. View at: Google Scholar
 J. W. Ma, X. X. Niu, H. M. Tang, Y. K. Wang, T. Wen, and J. R. Zhang, “Displacement prediction of a complex landslide in the Three Gorges Reservoir Area (China) using a hybrid computational intelligence approach,” Complexity, vol. 2020, Article ID 2624547, 15 pages, 2020. View at: Publisher Site  Google Scholar
 P. Deheuvels, “La fonction de dépendance empirique et ses propriétés — un test non paramétrique d'indépendance,” Académie Royale de Belgique  Bulletin de la Classe des Sciences 5e Série, vol. 65, pp. 274–292, 1979. View at: Google Scholar
 P. Deheuvels, “A non parametric test for independence,” Publications de l’Institut de Statistique de l’Université de Paris, vol. 26, pp. 29–50, 1981. View at: Google Scholar
 A. Gramacki, Nonparametric Kernel Density Estimation and its Computational Aspects, Springer, New York, NY, USA, 2017.
 C. Filippín, F. Ricard, and S. Flores Larsen, “Evaluation of heating energy consumption patterns in the residential building sector using stepwise selection and multivariate analysis,” Energy and Buildings, vol. 66, pp. 571–581, 2013. View at: Publisher Site  Google Scholar
 D. Kahneman, Thinking Fast and Slow, Farrar, Straus and Giroux, New York, 2011.
 M. D. Mastrandrea, K. J. Mach, G. K. Plattner et al., “The IPCC AR5 guidance note on consistent treatment of uncertainties: a common approach across the working groups,” Climatic Change, vol. 108, no. 4, pp. 675–691, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Shen Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.