Available online at www.sciencedirect.com
Procedia Environmental Sciences 12 (2012) 598 – 603
2011 International Conference on Environmental Science and Engineering (ICESE 2011)
Using Logistic Regression and GIS to Analyze the Relationship between Precipitation and Debris Flow in Sichuan, China Wenjuan Yua,Wenbo Xua,Guoping Zhangb a
Institute of Geo-Spatial Science and Technology University of Electronic Science and Technology of China Chengdu, Sichuan Province, 611731, China
[email protected] [email protected] a Public Weather Service Center of CMA Beijing, China
[email protected]
Abstract Debris flow is a serious disaster that frequently happens in Sichuan province China. The relationship between precipitation and debris flow is one of the important components of prediction research. We chose the geographic information system (GIS) as a tool to estimate the precipitation of hazard point and used statistical techniques to calculate attenuation coefficient of effective antecedent precipitation. With such methodologies, the logistic regression model was used to comparatively analyze the contribution rate of two types precipitation combination: 1) intraday rainfall and 10-day previous rainfall, 2) intraday rainfall and effective antecedent rainfall. Our results indicate that the contribution rate of intraday rainfall is the highest and that the category coincidence rate of second type precipitation is 3% higher than the first. © by Elsevier B.V. Selection and/or peer-review under responsibility of National University of Singapore. © 2011 2011Published Published by Elsevier Ltd. Selection and/or peer-review under responsibility of [name organizer] Keywords: debris flow; logistic regression; attenuation coefficient; effective antecedent rainfall
1. Introduction Debris flow is a sudden natural hazard which is caused by precipitation and a large volume of loose debris torrent solid in valley. The precipitation is an important factor for debris flow. It provides a trigger condition and acts as transport media.
1878-0296 © 2011 Published by Elsevier B.V. Selection and/or peer-review under responsibility of National University of Singapore. doi:10.1016/j.proenv.2012.01.323
Wenjuan Yu et al. / Procedia Environmental Sciences 12 (2012) 598 – 603
Previously, researches about relationship between precipitation and debris flow have been conducted, but the study areas mainly focused on small watershed and the prefecture-level cities, such as Ganzi, Liangshan, and Ya’an. The time span for such studies was as short as 10 years or below. Therefore, there have been insufficient information for analysis. Furthermore, the relation studies were conducted basing on sample statistical and qualitative descriptions [1-3] and researches method focused mostly on critical rainfall calculating. The conclusion of these researches is that debris flow was triggered by critical rainfall when it reached a certain value [4-6]. However, except for intraday rainfall, the occurrence of a hazard still depends on the antecedent precipitation. And because of the surface runoff, wind and other factors, antecedent precipitation would reduce along with time. So far, some researchers have calculated the attenuation coefficient by measuring the moisture content of soil [7]. Some used the empirical statistic value 0.84 as attenuation coefficient [8-9]. This method, however, was always employed in small watershed researches. While the third method was to calculate it by analyzing the relationship between occurrence time of debris flow and rainfall time [10]. In view of insufficiencies of the previous researches, e.g. short time span and limited study area, the mainly idea of this article is to analyze the relationship of regional rainfall and debris flow by a set of long time span precipitation data and debris flow data in large scale region. 2. Material And Method 2.1 Study Area Sichuan Province is located in southwest of China, and is bounded by longitudes of 97°21'E and 108°31'E, latitudes of 26°03'N and 34°19'N. Topography is complex and diversified, including mountainous area, basin and plateau. The climate of study area belongs to warm temperature climate zone, characterized by abundant and concentrated rainfall. The debris flow historical records from the China Geological Environmental Monitoring Institute shows that there are a total of 467 disaster points with precise position in the region from 1951 to 2004 2.2 Data In this paper, daily precipitation data that comes from 54 meteorological stations is used to research the relationship between rainfall and debris flow. For sufficient observational data for our study, we used the daily precipitation records during the period of 1981-2004 extracted from the China Meteorological Data Sharing Service System and the corresponding disaster data from the China Institute of Geo-Environmental Monitoring. The analysis factors include intraday rainfall, 10-day previous rainfall and effective antecedent rainfall. 2.3 Methodolody GIS technology and Logistic regression model were used to research. In reality, it is difficult to measure the rainfall of disaster points when the debris flows happen. Therefore, a useful method spatial interpolation is use to estimate the precipitation of disaster points. Our method is: 1) estimate the precipitation of disaster points, 2) calculate the attenuation coefficient of effective antecedent rainfall by statistical method, 3) use logistic regression model to analyze the relationship between debris flow and two types of comparative rainfall which are intraday rainfall and 10-day previous rainfall, intraday rainfall and effective antecedent rainfall. Fig. 1 demonstrates the flow of research process.
599
600
Wenjuan Yu et al. / Procedia Environmental Sciences 12 (2012) 598 – 603
Daily Precipitation Data Selection
Kriging Spatial Interpolation
Attenuation Coefficient and Effective Antecedent Rainfall
Intraday Rainfall and Antecedent Rainfall
Logistic Regression Model
Build the Equation
Discuss the Results
Fig. 1 Flow of Research Process
1αKriging Interpolation Kriging is an interpolation based approximate meta-model. It is a widespread and high-accuracy spatial interpolation method used to estimate rainfall of disaster points [11]. This method is used here to estimate the intraday rainfall and 10-day previous rainfall of points at which the disaster has happened. The calculation is carried out by the geostatistical analyst functional model in software ArcGIS Map. 2αMethod of Calculating Effective Precipitation The effective antecedent rainfall is calculated as a summation of daily rainfall, multiplied by a decrease coefficient as presented as the following equation: n
Ra
¦k R i
i
(1)
i
Where Ra is effective antecedent Rainfall and Ri (i=0,1,2…n) is the previous “ n day ” of rainfall. And k is attenuation coefficient which calculated by a statistical method. 3αLogistic Regression Logistic regression model is a logarithmic model that describes the relationship between one or more variables and binary dependent variables. The binary logistic model is as follows: P eE0 E1x1 !Emxm /1 eE0 E1x1!Emxm (2) Where P is the probability of an events occurring, ȕ0 is the constant of the equation, x1, x2,Ăxm are the values of the independent variables and ȕ1, ȕ2…ȕm are the coefficient calculated for independent variables. Actually, the following simpler format is always used:
P 1/(1 e B ) B E 0 E1 x1 ! E m xm
(3)
Here, the presence or absence of debris flow is used as dependent variable, which takes a value of 0 and 1 (for) to represent the absence and presence of disaster respectively. The intraday rainfall, 10-day previous rainfall and effective antecedent rainfall are used as independent variables to establish the equation.
Wenjuan Yu et al. / Procedia Environmental Sciences 12 (2012) 598 – 603
3. Results 3.1 Spatial and Temporal Distribution
Frequency
Rainfall (0.1mm)
At the Temporal distribution, 93.8% of debris flows occur in summer from June to September. In this flood season of Sichuan Province, the rainfall accounts for 70% of annual precipitation. Fig.2 shows the relationship among occurrence of time, frequency and rainfall which from 1981 to 2004 in Sichuan Province. It can be seen that rainfall increases the risk and frequency of the disasters.
Time (month) Fig. 2 Relationship among time, frequency and rainfall
At the Spatial distribution, hazards are all located in the heavy rain areas. The distribution structure is along the south-north direction, with obvious regionalization (see Fig. 3). In the eastern part, disasters occurred in the biggest rainfall area with low frequency. Heavy rainfall is the dominant factor in this part. The average annual rainfall around the Sichuan Basin and south is lower than eastern part, while the occurrence frequency is the highest in the study area. Except for the precipitation, debris flow is also affected by other factors such as geological features and other factors.
Fig. 3 Relationship Between space, Frequency and rainfall
3.2
Analysis of Antecedent Precipitation
As for the occurrence of debris flow, precipitation is a dynamic factor. Continuous rainfall makes the soil moisture, and decrease the shear strength of soil. Intraday rainfall and 10-day previous rainfall are used as independent variables. Regression result is as follows:
601
602
Wenjuan Yu et al. / Procedia Environmental Sciences 12 (2012) 598 – 603
P 1/1 e B .
(3) B 0.528 0.024r0 0.002r1 0.001r2 0.001r3 0 .0 0 2 r4 0 .0 0 1r6 0 .0 0 2 r8 0 .0 0 1r9 Where r0 represents the intraday rainfall which debris flow occurred, r1-r10 respectively 10 days before the debris flow rainfall. Category coincidence rate was 70%. It can be seen from the regression equation, the coefficient of intraday rainfall is the highest, 0.024. That means the contribution rate of intraday rainfall to debris flow is the biggest. 3.3
Analysis of Effective Antecedent Precipitation
The frequency of debris flow after the day with heavy rain has been counted. For example, “one day before” means the rainy event happened is one day earlier than the disaster (see Table I). Table I Rainy Time and Frequency Rainy time Total Frequency Intraday 275 83.586 One day before 29 8.814 Two days before 10 3.039 Three days before 10 3.039 Four days before Five days before Six days before Seven days before Eight days before Nine days before Ten days before
5 0 0 0 0 0 0
1.522 0 0 0 0 0 0
Frequency
As can be seen from Table I, there are 83.6% of the debris flow occurred after intraday rainfall. The following exponential relationship can be built by fitting a curve with the total frequency. As shown in Fig. 4 and (4), the fitting accuracy is as high as 99%.
Time (day) Fig. 4 Attenuation Coefficient Curve
f ( x) 83.5526 exp(2.1869 x) . (4) Where x is the day before debris flow occurred. For example, when x=0, it means the intraday when disaster happened, x=1, it means the previous one day. f(x) is the attenuation coefficient. The intraday rainfall and effective antecedent rainfall are used as independent variables. Result equation is as follows: P 1/1 e B
(5)
Wenjuan Yu et al. / Procedia Environmental Sciences 12 (2012) 598 – 603
B
0.411 0.044r0 0.024ra
Where, r0 is the intraday rainfall, ra is effective antecedent rainfall. Category coincidence rate of this group is 73%. And the contribution of intraday rainfall, which is still larger than other rainfall’s, is 0.044. 4. Conclusion The location of debris flow disasters and the distribution of rainfall are factors interrelated. At the Temporal distribution, debris flow occurs mainly in summer, the rainy season. At the Spatial distribution, the locations of debris flows are also within the rainy region. As a result, rainfall is the dominant factor of debris flows. The expression of the relationship between rainfall and the two types of comparative rainfall, i.e. intraday rainfall and 10-day previous rainfall, intraday rainfall and effective antecedent rainfall, comes out by Logistic regression analysis. The results indicate the coefficient of intraday rainfall is the highest. It means that intraday rainfall contribute most to the debris flow. According to comparison of the category coincidence rate between two types of rainfall, the second type is 3% higher than the first. It indicates that effective antecedent rainfall is more significantly related to debris flow. Acknowledgment This research work is funded by National Natural Science Foundation of China, project number: 40971016. References [1]
Su Pengcheng, Liu Xilin, Guo Jie, “Primary analysis of relationship between debris flows and rainfall in Sichuan
Province,” Journal of Natural Disasters, vol. 15, no. 4, pp. 19-23, August 2006. [2]
Yu Shuhua, “Analyses on spatial- temporal characteristics of mud-rock flow and landslip in Sichuan Basin and its
meteorological cause”, Plateau Meteorology, vol. 22, suppl, pp. 83-88, October 2003. [3]
B.F.Di, N.S.Chen, P.Cui, “GIS-based risk analysis of debris flow: an application in Sichuan, southwest China”,
International Journal of Sediment Research, vol. 23, no. 2, pp. 138-148, 2008. [4]
Cong Weiqing, Pan Mao, Li Tiefeng, “Quantitative analysis of critical rainfall-triggered debris flows”, China
Journal of Rock Mechanics and Engineering, vol. 25, supp. 1, pp. 2808-2812, February 2006. [5]
Servel Miller, Tim Brewer, Norman Harris, “Rainfall threshholding and susceptibility assessment of rainfall-
induced landslides: application to landslide management in St Thomas, Jamaica”, Geol Environ, vol. 68, pp. 539-550, 2009. [6]
Li Tiefeng, Cong Weiqing, “A method for rainfall-induced landslides prediction based on Logistic regression and
effective antecedent rainfall,” The Chinese Journal of Geological Hazard and Control, vol. 17, no. 1, January 2006 [7]
Wei Fangqiang, Hu Kaiheng, Chen Jie, “Determination of effective antecedent rainfall for debris flow forecast,”
Journal of Mountain Science, vol. 23, no. 4, pp. 453-457, July 2005. [8]
Cui Peng, Yang Kun, Chen Jie, “Relationship between occurrence of debris flow and antecedent precipitation:
taking the Jiangjia Gully as an example,” Science of Soil and Water Conservation, vol. 1, no. 1, pp. 11-14, March 2003. [9]
C.L.Shieh, Y.S.hien, Y.J.Tsai.etc, “Variability in rainfall threshold for debris flow after the Chi-Chi earthquake in
central Taiwan, China,” International Journal of Sediment Research, vol. 24, pp. 177-188, 2009. [10]
Ma Li, You Yangsheng, Miao Qilong, “The research of landslip forecast reduced by strong precipitation”, Journal
of Mountain Science, vol.26, no.5, pp.583-589, September 2008. [11]
Gao Ge, Gong Lebing, Zhao Shanshan, “Spatial interpolation method of daily precipitation”, Journal of Applied
Meteorological Science, vol. 18, no. 5, pp. 732-736, October 2007
603