^{1}

^{*}

^{2}

In this paper a rainfall threshold and a Bayesian probability model are presented for the landslide occurrence of shallow landslides in Ha Giang city and the surroundings, Vietnam. The model requires the data on daily rainfall combined with the actual dates of landslide occurrences. Careful study on the database is a prerequisite for the paper. For this reason, selecting the input data was carried out carefully to ensure the reliable results of the study. The daily rainfall data covering a time span of 57 years was collected from a unique rain gauge station of National Centre for Hydro-meteorological Forecasting of Vietnam (from 1957 to 2013) and a landslide database with some landslides (37 of total of 245 landslides) that containing dates of occurrence, was prepared from historical records for the period 1989 to 2013. Rainfall thresholds were generated for the study area based on the relationship between daily and antecedent rainfall of the landslide events. The results shows that 3-day antecedent rainfall (with the rainfall threshold was established: R
_{T} = 40.8
− 0.201R
_{3ad}) gives the best fit for the existing landslides in the landslide database. The Bayesian probability model for one-dimensional case was established based on 26 landslides for the period 1989 to 2009, daily rainfall data with the same time and the values of probability varies from 0.03 to 0.44. Next, the Bayesian probability model for two-dimensional case was generated based on 11 landslides, rainfall intensity and duration in three months (May, June and July) of 2013 and the values of probability ranges from 0.08 to 0.67, and computed values of conditional landslide probability
*P*(A|B) from two-dimensional case of Bayesian approach are clearly controlled by rainfall intensity > 40 mm with rainfall duration > 0.3 day.

Landslide occurs every where on over the world, especially in the small and medium river basin, along traffic road, mountain foot area in the mountainous areas in the countries under the influence of climate and storms.

Landslide studies have been carried out from very long time ago until now. The main objectives of these researches are to understand, manage, monitor, control, prevent and mitigate the huge impacts of this natural phenomenon to human’s life, socio-economic development, cultural infrastructure, environmental degradation…

To fulfill these objectives a lot of scientific articles related to this field were published. One of the most important landslide causative factors is rainfall data which are supposed as a triggering factor of landslide. The early studies about the rainfall intensities and duration associated with shallow landslides and debris flows were carried out by Caine [

Extreme rainfall conditions, when intensities remain high for long periods of time, must be accurately examined. In this case, any soil could become unstable. Hence, a better knowledge of the expected rainfall intensity and duration is fundamental for understanding the relationship with soil failure initiation and to establish different warning thresholds [

Since the studies of Caine [

Empirical methods for temporal probability assessment are based on the estimation of rainfall thresholds obtained by studying rainfall conditions that have resulted in landslides. They are usually contained in envelope curves based on variables such as cumulative rainfall, antecedent rainfall, rainfall intensity, and rainfall duration [

Recent studies introduced improvements using statistical methods for a more precise constrain of the ID threshold. This is the case of Bayesian inference method (e.g. [

In this paper, we propose a rainfall threshold for landslide initiation was estimated based on the relationship between the rainfall events and the historical landslide records. Then, the rainfall threshold was validated using the landslide events in 2000, 2007, 2008 and 2013. The method requires complete information on landslides including the dates of occurrence in order to correlate them with rainfall. Besides this method, Bayesian approach will be apply to estimate the probability of landsliding conditional to characteristics of rainfall events. The applications of these methods are based on the quality of landslide information (timing and location). So we used a data set consisting of 37 events that occurred between 1989 and 2013 from a database of 245 landslides occurred in the study area. Daily rainfall data were provided by a unique rain gauge station in the study area and collected rainfall data from 1957 to 2013.

This method is applied to assess landslide currently induced by rainfall in the Ha Giang city and the surroundings, one of the most high risk potential of landslide hazard in mountainous area of Vietnam as well as showing the level of the probability of landslide occurrence in the region. This study area was selected because of the availability of historical landslide records including the dates of occurrence of the landslides, and daily rainfall records from a rain gauge station.

Ha Giang city and the surroundings are located in Ha Giang, one of the northern mountainous provinces in Vietnam (^{2} between longitudes 104˚52'E and 105˚18'E and latitudes 22˚34'N and 22˚56'N. The elevation in the region ranges from 68 to 2000 m above sea level and the mountainous region is strongly dissected and steep. It includes high mountains in the north, south and west but the distribution is not equability and the altitude is gradually decreases to the middle of the region. Moreover, the

karstic area is also representative of the north part of study area. Beside the drawback of terrain, here is also a tectonically active area with the complex of fracture zones, faults system and other important feature is dense drainage network with river and streams that landslides occur frequently as one of the most common natural disasters.

Geologically, the area comprises of limestone, sandstone, silty sandstone, black clay shale, quartz-biotite schist, sericite-chlorit schist, mica schist, green schist, and some kind of rock with small distribution such as: porphyritic biotite granite, granosyenite, olivine gabbro, gabbro diabase, conga diabase… The ages of rocks vary from Proterozoic to Trias. The bedrock of this study area is weathered strongly, the average thickness of weathered layer is 3 - 5 meters, and the most thickness layer is 20 meters as observed in the cut slopes along the transportation road in the region.

The contact between bedrock and weathering soils is often exposed in the cuts slopes, which make them more susceptible to landslides due to the build-up of pore pressure on the contact [

The land use of the study area comprises of forest land, barren land and non-forest rocky mountain, agricultural land, residential areas, bush and grass areas and water surface.

The study area is situated in the monsoonal region, with four seasons. The coldest month is January and the hottest month is June with average temperature of 14˚C and 27˚C, respectively. Seasons in the region can be classified as rainy and dry. Rainfall is main occurs in rainy season. According to the National Centre for Hydro-meteorological Forecasting of Vietnam recorded from 1957 to 2013, the average annual rainfall of this area ranging from 2500 to 3300 mm/year where 90% of the total rainfall occurs in the summer (from May to October every year). Moreover, the frequency and intensity of the rainfall are concentrated over a short period that triggers most of the disastrous hazards in the study area.

Landslide database was collected from historical records of the previous projects. This data is based on landslide inventory from several projects such as: 1) “Investigation and assessment of the types of geological hazard in the territory of Vietnam and recommendation of remedial measures”. Phase II: A study of the northern mountainous province of Vietnam [

In the study area numerous landslides normally occur in the period from May to October. However, almost of landslide points which were collected from field work of previous projects are lacked of exactly dates of landslide occurrence. Some difficulties was mentioned by Hung [^{3} and the maximum size is 9450 m^{3} and median = 134 m^{3} and the total of volume of landslides is approximately 150.000 m^{3} but this type of natural disaster have been increasingly severe in term of magnitude, frequency and volatility, causing vast losses in human life, property, socio-economic and cultural infrastructure as well as environmental degradation. As a mention above, most of the landslides are undefined landslides it means that the surveyor could not collected the exactly date of landslide initiation. Fortunately, from historical landslide records we had list of 37 landslide events that knew exactly when landslide occurred. Inside, 26 landslide-triggering events occurred from 1989 to 2009 and the

remaining 11 events occurred in 2013. Most of landslides occurred in July with 20 events, the next is in June with 10 events, both of May and August are having the same times with 3 events, and then in September with 1 events. Analysis of the landslide inventories shows that most of the slope failures in the study area were caused by rainfall infiltration into the soil causing an increase in soil pore-water pressure [

For the study area, daily rainfall data were collected from a unique rain gauge belonging to the National Centre for Hydro-meteorological Forecasting of Vietnam from the period of time 1957 to 2013. The location of this rain gauge is shown in

The total of yearly rainfall of study area was from 1957 to 2013 ranging from 2500 to 3300. The study of daily rainfall records reveals that the area experiences rainfall in two periods: rainy season, that the rainfall is always concentrate in May to October and the common highest month of rainfall is often July and dry season from November to April. The lowest recorded annual rainfall is 1465 mm and the highest is 3305 mm (

Although in Vietnam, many landslide project have been completed during 20 recent years but we are really lack of some important information of historical landslide data. The first reason of this fact is because all the landslide field work of the previous projects was carry out in April or October. April is at the end of period of time of dry season and October is at the end of period of time of rainy season. This works is not suitable for collecting the landslide information because almost of landslide in mountainous area of Vietnam is always occurs in rainy season from May to October, especially this natural phenomenon is occurred with the most frequently in 3 months from June to August. Due to the fact that, the landslide information collected is very hard to link with exactly date of rainfall episode. Although in the study area collected 245 landslides but most of them could not define the exactly date except 37 landslides triggering events as mention above. Most of these landslides have occurred in areas where the rainfall was relatively high with average daily rainfall occurred landslide is more than 76 mm, even having 9 of 37 landslides occurred the same day with maximum daily rainfall.

One of the most difficulties when using antecedent rainfall measurements to predict landslide occurrence is the definition of the period over which to accumulate the precipitation [

A detailed literature review revealed a complex relationship on the correlation between numbers of days for the antecedent rainfall with the triggering of a landslide [

The thresholds using rainfall intensity-duration are the most widely used method in the literature [

In this study area where daily rainfall (covering a time span of 57 years from 1957 to 2013) and hourly rainfall data (only from May to October in the period of 4 years from 2010 to 2013) are available, so antecedent rainfall plays an important role in estimation of landslide initiation introduced by Jaiswal and Van Westen [

As discussed in Section 2, landslide hazard in the area having strong relationship with certain intensity of rainfall. Hence, we suppose that the probability of occurrence of a landslide is related to the probability of occurrence of the triggering rainfall threshold. It means that landslide has to occur whenever a given rainfall threshold is exceeded may not hold always and everywhere. However, it is also expected that landslides will not occur below the rainfall threshold. Hence, for rainfall-triggered landslides, this assumption can be an acceptable first-approximation to work with and to estimate the frequency of landslides by establishing relations between the landslide trigger, its magnitude and the occurrence of the landslides [

Depending on the type of landslides and their geo-environmental setting, the number of antecedent days can vary from 3 days for shallow landslides to 30 days for deep landslides [

To determine the rainfall threshold (R_{T}), a scatter plot was prepared showing daily rainfall against the corresponding 3-day antecedent rainfall, for each day with one or more triggered shallow landslides. The envelope curve is manually drawn such that it demarcates the lower end of the plotted points. The line can be represented by a linear mathematical equation [

In this study area, R_{T} above which a landslide can occur for the given 3-day antecedent rainfall (R_{3ad}) is represented by the equation R_{T} = 40.8 − 0.201R_{3ad}. Besides, we also established the equations of rainfall threshold for 5, 7, 10 and 15-day antecedent rainfall and the landslides that occurred in study area were used for the threshold validation (

In some cases, landslides were also occurred in the days which no rainfall or with very few of daily rainfall was measured. These were the cases when high antecedent rainfall alone has resulted in landslides and other reason is the failures to pore pressure rising due to water percolating from upslope areas. This holds for landslides associated with cut slopes because during excavation toe of such slopes are removed and the unsupported overburden mass becomes more prone to failure under the given condition. Thus, for all the listed thresholds, the lower boundary of the envelope curve was set to zero daily rainfall [

The rainfall thresholds can be used to predict landslides both spatially and temporally. The temporal aspect is related to the daily variations in rainfall and the spatial aspect can be related to the use of different thresholds for different areas. An indirect way to test the predicting capability of the thresholds is to validate them with the control data sets, which were not used in the model [

In order to validate the model, the historical landslide recorded in some years such as 2000, 2007, 2008, and 2013 were used. A part of the landslide events in these years were also used in the rainfall threshold model creation.

According to the results shows in

Bayesian subjective probability models, forecasting models can be constructed using statistical methods applied to objective data. However, when historical data are unavailable, irrelevant, or there are time and financial restrictions, a subjective approach can be used. The Bayesian subjective probability model relies on querying the opinions and intuitions of a group of experts, obtaining probability estimates, and processing them through the Bayes formula. It should be noted that subjective models are limited in that they are slow in revealing previously unrecognized relationships. Experts can err in making complex forecasts, especially if they have to consider conflicting clues. “Clues” are defined as separate independent facts that impact on the forecast. Because the number and form of clues varies, it is therefore important to aid the experts’ judgments with a methodology to organize them. The Bayes probability model provides a structure for breaking a complex judgment into several smaller components and reducing expert error. Construction of the Bayes probability model requires four steps: 1) Decide on events to forecast; 2) Identify the clues; 3) Assess likelihood estimates; 4) Use Bayes formula to forecast [

Because of only have a unique rain gauge in the study area, so we selected Bayesian probability model applied for one-dimensional case to calculate landslide probability caused by daily rainfall (mm/day). According to Berti [

Bayes’ theorem is a direct application of conditional probabilities. The conditional probability is the probability of some event A (in our case a landslide) given the occurrence of some other event B (a rainfall episode with a certain magnitude, expressed in terms of total rainfall, intensity or any other variable). Conditional probability is written P ( A | B ) and it is read “the probability a landslide (A) occurs given a rainfall episode (B)”. This probability is provided by the Bayes’ theorem:

P ( A | B ) = P ( B | A ) ⋅ P ( A ) P ( B ) (1)

where:

P ( B | A ) = conditional probability of B given A (also called the likelihood), that is the probability of observing a rainfall event of magnitude B when a landslide occurs,

P(A) = prior probability of A (or simply prior), that is the probability a landslide occurs regardless of whether a rainfall event of magnitude B occurs or not,

P(B) = marginal probability of B, that is the probability of observing a rainfall of magnitude B regardless of whether a landslide occurs or not,

P ( A | B ) = conditional probability of A given B (also called posterior probability), that is the probability of observing a landslide when a rainfall event of magnitude B occurs.

Bayesian probability is usually computed in terms of relative frequencies. Thus, if N_{R} is the total number of rainfall events recorded during a given time reference; N_{A} is the total number of landslides occurred during the same period; N_{B} is the number of rainfall events of magnitude B; and N ( B | A ) is the number of rainfall events of magnitude B that resulted in landslides, the probability terms in (1) can be approximated to:

P ( A ) ≈ N A / N R (2a)

P ( B ) ≈ N B / N R (2b)

P ( B | A ) ≈ N ( B | A ) / N A (2c)

And Equation (1) reduces to P ( A | B ) ≈ N ( B | A ) / N B .

Equation (1) can be easily extended to the case of two variables B and C:

P ( A | B , C ) = P ( B , C | A ) ⋅ P ( A ) P ( B , C ) (3)

where: the notation B, C indicates the joint probability of having a certain value (or range of values) of the two variables. If, for example, B ≡ I is rainfall intensity and C ≡ D is rainfall duration, Equation (3) provides the probability of a landslide in response to a rainfall event of given duration and intensity.

In this paper, we selected Bayesian probability model applied for two-dimensional case to calculate landslide probability caused by two factors, rainfall intensity (mm/day) and rainfall duration (days).

According to Driver [

The first step in the evaluation of any rainfall threshold is to identify the rainfall episodes that triggered the historical landslides, here referred as “triggering rainfall”. Ideally a triggering rainfall event should be a well-defined rainfall episode, described by its duration, amount of precipitation, and intensity and clearly related to a given landslide. In some cases the identification is simple (for instance if the landslide occurred after a heavy rainfall preceded by a prolonged dry period) but usually it is not. Landslides may result from complex rainfall sequences made of multiple bursts of variable duration and intensity that make it difficult to detect a well-defined triggering episode. The greatest uncertainties usually derive from the identification of the beginning of the triggering rainfall while the time of landsliding is taken as its end. A certain amount of time without rainfall (or limited rainfall) can be used as criterion to truncate the rainfall sequence that precedes a landslide event [

In this paper, we concentrated on estimating the probability of rainfall threshold for landslide occurrence using daily rainfalls as the triggering rainfall events. For each of the 37 historical landslides (see 2.1), we compared the rainfall data recorded by the reference rain gauge with those recorded. This work was done to find out the relationship between rainfall cumulative data with landslide occurrence in the past events. The triggering rainfall was then defined visually by selecting the rainfall episode closest to the date of occurrence of the landslide. The results were then compared and the discrepancies (usually related to the beginning of the triggering event) were discussed to arrive at a shared definition. We generally agreed to define the triggering event as a period of continuous or nearly continuous precipitation which starts with the onset of the rainfall (or with an abrupt increase of rainfall intensity in a period of light rain) and ends the day of occurrence of the landslide. For those landslides that have occurred after the end of the rainfall, the duration of the triggering event was set equal to the rainfall duration [

Based on the classification of Berti [

Type 1: Well-defined rainfall events can be clearly identified, as shown in the examples of

Type 2: Uncertain events consist of distinct rainfall episodes characterized by uncertain or subjective limits because of the presence of secondary rainfall episodes (

Type 3: Undefined events consist of all those landslides without a significant rainfall event close to the date of occurrence, such as landslides triggered during complex rainfall sequences (

In order to calculate the values of P(A), P(B), P(B|A), P(A|B) we selected only well-defined events (Type 1) for considering in Bayesian analysis. When using this method trigger rainfall must be reliably identified in order to get a reliable likelihood function P(B|A). Landslides triggered by other factors (Type 2 and 3) were then omitted [

Among 37 landslides triggered by rainfall events in the study area, we selected 26 well-defined rainfall events for the input data for Bayesian approach.

The analysis can be applied to the study area data set following the procedure described in Section 4.2.

The marginal rainfall probability P(B) is computed using the N_{R} = 346 rainfall events recorded (we only used rainfall data in the months that occurred landslide triggering rainfall events) in the period of 13 years (May of 1989; May and June of 1993; June of 1995; August of 1996 and 2008; July of 1997, 1998, 2000, 2001, 2004 and 2009; June and July of 2006; June, July and September of 2007). The conditional probability P ( A | B ) using the N_{A} = 26 well-defined rainfall that triggered the historical landslides. Prior landslide probability is therefore P(A) = N_{A}/N_{R} = 26/346 = 0.08. In

Due to lack of historical landslide records in three years, 2010, 2011 and 2012, we decided not using the landslide database of 2013. This data will be combined

Total rainfall events | Number of daily rainfall intensity with 0 < I < 40 mm/day | Number of landslide with 0 < I < 40 mm/day | Number of daily rainfall intensity with 40 < I < 65 mm/day | Number of landslide with 40 < I < 65 mm/day | Number of daily rainfall intensity with I > 65 mm/day | Number of landslide with I > 65 mm/day | Total landslide |
---|---|---|---|---|---|---|---|

346 | 293 | 8 | 21 | 4 | 32 | 14 | 26 |

Rainfall intensity classification | Rainfall events | P(B) | P(B|A) | P(A) | Probability P(A|B) | Total landslide | Total rainfall events |
---|---|---|---|---|---|---|---|

0 < I < 40 | 293 | 0.85 | 0.31 | 0.08 | 0.03 | 26 | 346 |

40 < I < 65 | 21 | 0.06 | 0.15 | 0.08 | 0.19 | ||

I > 65 | 32 | 0.09 | 0.54 | 0.08 | 0.44 |

with two types of rainfall data that are daily rainfall and hourly rainfall for estimating the values of probabilities in two ways one-dimensional and two-dimensional Bayesian analysis for only 2013. The results of one-dimensional Bayesian analysis shows in

Two-dimensional Bayesian analysis evaluates the conditional probability of the event given the joint occurrence of two control variables (Section 4.2).

In order to calculate the values of probability in two-dimensional Bayesian probability model requires the data collection in more detail and based on the rainfall intensity and duration. Because of the hourly rainfall data are not frequently available especially for a long period of time in the past. Therefore, we only had collected this type of data in 4 years, from 2010 to 2013 in 6 months of rainy seasons (from May to October). Although, our landslide database be collected to 2013 but these are unavailable in during three years from 2010 to 2012. So that, we have to be use the landslide data in three months, from May to July of 2013 for calculating the parameters of two-dimensional Bayesian probability model.

The application of probability model follows the Equation (3) described in Section 4.2. The hourly rainfall events recorded in three months of May, June and July in 2013 was 68 events, the total landslide events occurred was 11, total duration happened total rainfall events was 18.46 days (shortest duration was 0.04 day/a rainfall event, and the longest duration was 0.67 day/a rainfall event). To evaluate the conditional probability of the event given the joint occurrence of two control variables are hourly rainfall events and its duration, we classified “rainfall intensity-rainfall duration, I-D” into four categories for calculating easily. Four different combination of rainfall intensity and duration of study area classified basing on the equation of probability of occurrence of the triggering

Total rainfall events | Number of daily rainfall intensity with 0 < I < 40 mm/day | Number of landslide with 0 < I < 40 mm/day | Number of daily rainfall intensity with 40 < I < 65 mm/day | Number of landslide with 40 < I < 65 mm/day | Number of daily rainfall intensity with I > 65 mm/day | Number of landslide with I > 65 mm/day | Total landslide |
---|---|---|---|---|---|---|---|

51 | 40 | 3 | 5 | 3 | 6 | 4 | 10 |

Rainfall intensity classification | Rainfall events | P(B) | P(B|A) | P(A) | Probability P(A|B) | Total landslide | Total rainfall events |
---|---|---|---|---|---|---|---|

0 < I < 40 | 40 | 0.78 | 0.30 | 0.20 | 0.08 | 10 | 51 |

40 < I < 65 | 5 | 0.10 | 0.30 | 0.20 | 0.60 | ||

I > 65 | 6 | 0.12 | 0.40 | 0.20 | 0.67 |

rainfall threshold of the study area which established in Section 3.2 with R_{T} = 40.8 − 0.201R_{3ad}.

1) Rainfall events has I ≤ 40 mm/day and D ≤ 0.3 day,

2) Rainfall events has I ≤ 40 mm/day and D > 0.3 day,

3) Rainfall events has I >40 mm/day and D ≤ 0.3 day, and

4) Rainfall events has I >40 mm/day and D > 0.3 day.

Based on

Classification of rainfall intensity-duration (RID) | RID did result in landslides | RID did not result in landslides | Total of RID events of each class | Total of landslide events | Total of RID events |
---|---|---|---|---|---|

I ≤ 40 & D ≤ 0.3 | 4 | 39 | 43 | 11 | 68 |

I ≤ 40 & D > 0.3 | 2 | 12 | 14 | ||

I > 40 & D ≤ 0.3 | 0 | 0 | 0 | ||

I > 40 & D > 0.3 | 6 | 5 | 11 |

Classification of rainfall intensity-duration (RID) | P(I|D) | P(I,D|A) | P(A) | P(A|I,D) | Total of landslide events | Total of RID events |
---|---|---|---|---|---|---|

I ≤ 40 & D ≤ 0.3 | 0.63 | 0.36 | 0.16 | 0.09 | 11 | 68 |

I ≤ 40 & D > 0.3 | 0.21 | 0.18 | 0.16 | 0.14 | ||

I > 40 & D ≤ 0.3 | 0.00 | 0.00 | 0.16 | 0.00 | ||

I > 40 & D > 0.3 | 0.16 | 0.55 | 0.16 | 0.55 |

According to the main objectives of the paper, we proposed a rainfall threshold model for landslide initiation based on the relationship between the rainfall events and the historical landslide records. Besides this method, a Bayesian approach to estimate the probability of landsliding conditional to characteristics of rainfall events to shows the possible risk of landslide occurrence.

In this paper, a threshold model based on antecedent rainfall have selected, because only daily rainfall data are available for a long time period and easy to implement. Data required for the model are derived from the landslide events in the past and daily rainfall records. The selection of the most suitable rainfall threshold model for the study area is based on historical landslide records and daily rainfall records in the period of 13 years, and the works to establish rainfall threshold equations from the relationships of landslides triggered by the rainfall events with antecedent rainfall in 3 days, 5 days, 7 days, 10 days and 15 days. Moreover, a validation of the rainfall threshold was also implemented and the results of this work has supported our selection is appropriate for the study area.

Jaiswal and Van Westen [

Selection a rainfall threshold method that building based on rainfall intensity and duration for this study area is not suitable because of the limitation of historical landslide records, rainfall records as well as rain gauge station network.

Bayesian probability model applied to calculate the probabilities of landslide occurrence in the study area follow two cases, one-dimensional case and two-dimensional case. In the first case, we determined that the selection of the input data is the most important steps. Therefore, only 26 landslide events in the past selected to take part in the calculation of the model. Beside, the ideal to classify the rainfall intensity threshold applied for the study area was also evaluated carefully basing on the fact of data. According to the rainfall threshold equation that we suggested for the study area in Section 3.2, the classification of rainfall intensity threshold to compute values of conditional landslide probability was applied from this rainfall threshold equation. In the second case, due to the unavailable of hourly rainfall database in the period of 13 years. The application of Bayesian probability model in the case of two-dimensional cannot carry out in this period. Hence, only the hourly rainfall data in May, June and July of 2013 (a unique hourly rainfall data that we collected) was took part in calculating the values of probability of this case. In this area, two variables that used for evaluating of two-dimensional Bayesian probability were rainfall intensity and rainfall duration. The classification of rainfall intensity threshold was also based on the same way in one-dimensional Bayesian probability to compute the values of the conditional probability.

Computed values of conditional landslide probability-P(A|B) in one-dimensional case in the period of 13 years in comparison with the values of its in May, June and July of 2013 have the differential levels relative clearly. These differences can explain as a simple way by the unavailable of the database.

Landslides are generally isolated processes which individually may not be very large in size but which can occur with a high frequency in a region. Unlike hazardous events that affect large areas, such as earthquakes or flooding, the generation of landslide inventory maps and databases is a tedious procedure. Landslides have to be mapped and described one by one, and each one might have different characteristics. In most countries there is no single agency that has the responsibility for maintaining a landslide database. At best, several organizations, such as the public works department or road department, will deal only with those landslides that have affected their area of interest (e.g., road network), resulting in incomplete and biased databases. Newspaper and other historical records only record those events that caused substantial damage. Universities and research institutes do work on landslide inventory maps but mostly as part of a research project, with a limited duration, after which a landslide database is no longer kept up-to-date. Therefore, it is very difficult to obtain landslide inventory maps that are complete, both with respect to the area covered and to the time period investigated [

Determining temporal probability is done either by correlating the data on landslide occurrences with those of the triggering factors (provided that the historical records are sufficient for this) or through dynamic modelling [

With the limitations of database and a lack of rain gauge station network that can be distributed regularly in the study area, so we cannot produce a map on spatial landslide probability distribution as in the research results of Berti [

Finally, we still believe that the results of this paper can helps local authorities and local people to be aware of the negative impacts of landslides triggered by rainfall events as well as help them in the selection the most suitable development planning methods for the study area to avoid badly effects, risk mitigation and damages in the future.

In the study area, landslide triggering is significantly controlled by an empirical minimum rainfall threshold: R_{T} = 40.8 − 0.201R_{3ad}. The Bayesian probability model for one-dimensional case was established based on 26 landslides for the period 1989 to 2009; daily rainfall data with the same time and the values of probability varies from 0.03 to 0.44. Next, the Bayesian probability model for two-dimensional case was generated based on 10 landslides; rainfall intensity and duration in three months (May, June and July) of 2013 and the values of probability ranges from 0.08 to 0.67, and computed values of conditional landslide probability P(A|B) from two-dimensional case of Bayesian approach are clearly controlled by rainfall intensity > 40 mm with rainfall duration > 0.3 day.

We have proposed an approach by combining an empirical minimum rainfall threshold model with a probability model for determining probability of landslide occurrence over an area by using Bayesian probability method. In this model, the value of minimum rainfall threshold of the study area is an important referential parameter related to rainfall intensity and rainfall duration for the application of Bayesian probability method in calculating the probability of landslide occurrence.

The model may not be applicable if the exact dates of landslide events in the past are not known. The results of the model will be better if the input data on daily rainfall, hourly rainfall from a well-distributed network of rain gauges is combined with the actual dates of landslide occurrences.

The computed landslide probability is dependent strongly on a well-distributed network of rain gauges in a well-defined area.

From the results of this research, we also want to continue our studies in other mountainous area of Vietnam as future works. To fulfill this objective, the installation of a well-distributed network of rain gauges in mountainous area of Vietnam is a strong recommendation of us in order to test the practical effectiveness of the proposed method in a wide variety of geological conditions.

We acknowledge the help of the Vietnam Institute of Geosciences and Mineral Resources (VIGMR) for providing relevant data and support. This paper was carried out as a part of the first author’s PhD studies at the Faculty of Engineering, China University of Geosciences (Wuhan), China. The authors also thank the three anonymous reviewers of Open Journal of Geology for relevant comments and Joanna Cheng, the Managing Editor, for assistance.

Do, H.M. and Yin, K.L. (2018) Rainfall Threshold Analysis and Bayesian Probability Method for Landslide Initiation Based on Landslides and Rainfall Events in the Past. Open Journal of Geology, 8, 674-696. https://doi.org/10.4236/ojg.2018.87040