Application of the Multitype Strauss Point Model for Characterizing the Spatial Distribution of Landslides

Landslides are common but complex natural hazards. They occur on the Earth’s surface following a mass movement process. This study applies the multitype Strauss point process model to analyze the spatial distributions of small and large landslides along with geoenvironmental covariates. It addresses landslides as a set of irregularly distributed point-type locations within a spatial region. Their intensity and spatial interactions are analyzed by means of the distance correlation functions, model fitting, and simulation. We use as a dataset the landslide occurrences for 28 years from a landslide prone road corridor in the Indian Himalayas. The landslides are investigated for their spatial character, that is, whether they show inhibition or occur as a regular or a clustered point pattern, and for their interaction with landslides in the neighbourhood. Results show that the covariates lithology, land cover, road buffer, drainage density, and terrain units significantly improved model fitting. A comparison of the output made with logistic regression model output showed a superior prediction performance for the multitype Strauss model. We compared results of this model with the multitype/hard core Strauss point process model that further improved the modeling. Results from the study can be used to generate landslide susceptibility scenarios. The paper concludes that a multitype Strauss point process model enriches the set of statistical tools that can comprehensively analyze landslide data.


Introduction
Landslides are defined as the movement of a mass of rock, debris, or soil along a downward slope, due to gravitational pull. The inherent properties of the Earth material, encompassing various geoenvironmental factors, can make a particular area susceptible to landslides. Landslides are among the most common natural hazards. They exhibit themselves in different mass movement processes and are considered as complex natural hazards occurring on the Earth's surface [1]. Although individual landslides are not as spectacular or damaging as earthquakes, floods, and hurricanes, they are widespread and frequently occurring. Over the years they have caused great loss of life and property and their effects on the economy are growing at a rapid pace in many countries.
Spatial zonation of landslide occurrences, also known as landslide susceptibility mapping, aims to differentiate a land surface into homogeneous areas according to their degree of failure caused by mass movement at specific locations [2]. It relies on understanding complex mass movement processes and their controlling factors [3]. Approaches to the spatial modeling of landslides can broadly be divided into two groups [4]. The first approach consists of deterministic, dynamic modeling of the physical mechanisms that control slope failure, using mathematical methods. This approach is highly localized because of the detailed data requirements. The second approach uses the relation between the locations of previous landslides and geoenvironmental variables, to predict areas of different landslide susceptibility, using heuristic or statistical methods. The statistical methods used successfully in landslide susceptibility mapping to date include discriminant analysis [5,6], multivariate statistics [7], likelihood ratio [8], information value method [9], and logistic regression methods [10]. These methods 2 Mathematical Problems in Engineering allow the analysis of geoenvironmental variables controlling landslide occurrence with respect to previous landslides without looking at the mutual interactions of landslides and their distribution patterns. Commonly applied generalized linear modeling uses a maximum likelihood estimation that results into point parameter estimates with standard errors [11]. As individual landslides cover only a small fraction of an unstable area, landsliding can be considered as a spatial point process that is controlled by a number of surface and subsurface spatial variables.
A spatial point process underlies a pattern of spatial point data within a region. Spatial point patterns are characterized by the 1st-and 2nd-order effects of a point process [12], specifying the intensity and point interactions, respectively. Nearest neighbor functions between pairs of points commonly model the 2nd-order effects as a function of positions of points and their distances. Such functions usually consider the relative position of two points in a bounded region and are taken as a function of distance only [13]. Spatial point processes play a fundamental role in spatial statistics and exhibit an active area of research. Disciplinary applications occur in forestry addressing positions of trees in a forest and log-landing sites [13][14][15], ecology addressing locations of bird's nests [13], seismology addressing earthquake epicentres [16,17], astrophysics addressing locations of stars in a constellation [18], and environmental modeling for peak concentrations of a pollutant in a geographical region [19]. In this study we use a spatial point process model to analyze the landslide data for inferring properties of the spatial distribution pattern.
The objective of this study was to identify significant factors for landslide susceptibility by applying a multitype Strauss point model. In this way the landslide occurrence patterns could be better understood. The point process model used information at the level of detail offered by the landslide data. That information was combined with geoenvironmental covariates to explain the underlying process for determining landslide susceptibility within an area. The model output was compared with the output of logistic regression model. The study was applied to a landslide prone area in the northern Himalayas, India.

Methods
A landslide distribution pattern can be considered as a collection of point data spread irregularly in an essentially planar region. A basic assumption for the analysis is that the data can be regarded as a realization of a stochastic point process [13]. This process is characterized by the intensity at location , ( ), defined as where | | is the area of a small region , is expectation operator, and ( ) is number of points in the region .
Landslide occurrences, considered as a spatial point pattern, show variation in the relative frequency as a function of the distance between positions Assuming stationarity, (2) depends only upon the relative position of the two landslides between positions and , that is, on their distance [13]. A landslide process is second-order stationary if its intensity is independent upon translation over , so that ( ) = and the second-order intensity depends only upon the distance vector between and and not on their locations.

Conditional Intensity and the Gibbs Model.
For a landslide process that exhibits inter-landslide interaction the pairwise interaction models define the intensity in the form of probability densities: where is a normalizing constant, ( ), ∈ is the intensity or the first-order term, and ( , ), , ∈ , is the pairwise interaction or second-order term in a bounded window . Pairwise interaction models are a special case of Markov point process, called Gibbs point process models. Analysis of a Gibbs point process model is based on its conditional intensity [20]. The probability of the occurrence of a landslide ( ) at location is determined by the conditional intensity defined by ( , ). For a landslide process in a bounded area the conditional intensity is related to the probability density by For the general pairwise interaction process the conditional intensity is Landslides of a particular zone with a specified radius of influence, however, might have a homogeneity condition with respect to each other. A Strauss process emphasizes such homogeneity conditions for deriving the relationship among the events. For the Strauss process, a simple model of dependence between landslides has the conditional intensity where ( , ) is the number of points of the landslides that lie within a distance of the location . The conditional Mathematical Problems in Engineering 3 intensity is useful in modelling because the two distinct components in the functional form represent the interaction of landslides that can be interpreted in a straightforward way, following Baddeley [20]. A major restriction of Gibbs models is that the parameters cannot be estimated using maximum likelihood estimation, and hence a maximum pseudolikelihood is returned for each model [20]. The unknown scaling factor is intractable; hence the calculated pseudolikelihood does not involve any unknown factor and it is easier to use when estimating the parameters.

Nearest Neighbour -Function.
Methods based upon the distances between landslides can be used for investigating inter-landslide interactions, for example, to identify secondorder effects of the landslide pattern data. These secondorder properties are specified by the pair correlation function that is assessed using the inter-landslide interaction methods like the -and -functions [20]. In this way, the nature of the departure from complete spatial randomness (CSR) can be identified. This, in turn, is useful in determining the kind of interaction and interaction distances between the landslides [21]. Alternatively, Ripley's -function can be used for detecting deviations from spatial homogeneity. The shape of this function indicates the specific type of pattern displayed by the data, that is, indicating whether the landslides show inhibition or occur as a regular or a clustered pattern.
The -function quantifies the distance distribution of a landslide to the nearest other landslide. It is expressed aŝ where is the radius of a disk centered at the location of the th landslide , ( , ) is an edge correction weight such that ( ) is approximately unbiased, = min ̸ = ‖ − ‖ is the distance of each landslide location to its nearest neighbor, and { ≤ } is the indicator function equal to 1 if ≤ and it is 0 otherwise.
An estimate of derived from a spatial landslide pattern dataset can be used in exploratory data analysis and formal inference about the pattern [22,23]. The shape of this function provides information about the way the landslides are clustered in a particular area. If the landslides are clustered, increases rapidly at short distances, whereas for landslides that are evenly spaced, increases slowly up to the distance at which most events are spaced, and only then it increases rapidly [24]. For a homogeneous Poisson point process of intensity , the nearest neighbor distribution distance function for landslide distribution equals If ( ) > pois ( ) then the landslide pattern is clustered, whereas if ( ) < pois ( ) the landslide pattern is regular.

Strauss Point Process for Marked Point Pattern Analysis of Landslides.
A multitype pairwise interaction process is a Gibbs process which assumes symmetric interactions of landslides have the probability density of the form where is a normalizing constant, ( ) is a function determining the first-order trend for landslides ( ) of the th type of patterns, indicating whether the landslides show inhibition or occur as a regular or a clustered pattern, and , ( , ) are symmetric functions that describe the interaction between a pair of landslides and of given types and ; that is, , ( , ) = , ( , ) and , ≡ , [20]. Thus, the conditional intensity of two categories of landslides Small and Large as defined in (11) for a multitype Strauss process is given by The multitype Strauss process has pairwise interaction terms where , > 0 are interaction radii as above and , ≥ 0 are interaction parameters [20].
To fit a multitype Strauss process model to landslides, the matrix of interaction radii , between individual landslides is specified based on field conditions and distribution patterns. Model fitting generates values for the interaction parameter , and the model coefficients describing spatial inhomogeneity and inter-landslide interactions.

The Multitype/Hard Core Strauss Point Process Model.
The multitype/hard core Strauss point process is a hybrid of the multitype Strauss process and the hard core process, for example, the case that , = 0, , = 0, or , = 0 of the Strauss process. A pair of points of types (Small) and (Large) must not lie closer than ℎ , units apart; if the pair lies more than ℎ , and less than , units apart, it contributes a factor , to the probability density similar to the , pairs and the , pairs. For landslides, this extension makes sense, as the large landslides are usually well separated and small distances do not occur; the same applies to the smaller landslides, although to a lesser degree as they can be closer. Moreover, a small landslide does not occur below a large one, whereas, if the location of a landslide is indicated as a point, the landslides cannot occur immediately close to each other.

Goodness of Fit.
Akaike's Information Criterion (AIC) is a measure of the goodness of fit of an estimated statistical model. The AIC is defined as AIC = −2ln ( ) + 2 , where is the maximum likelihood value for the estimated model and is the number of parameters in the model. It can be interpreted as the trade-off between bias and variance in model construction indicating that of accuracy and complexity of the model [25].
The study system is based on the point process modeling of landslide data along with the geoenvironmental covariates influencing landslide such as lithology, topography, or geology. The models are selected based on the AIC. The study makes use of the spatial distribution of landslide to make the exploratory data analysis, their interactions using -and -functions, and the possible susceptibility intensity in the area using the Strauss point model. The AIC is a test between models and hence it may serve as a tool for model selection. Given a dataset, several competing models may be ranked according to their AIC, with the one having the lowest AIC being the best [25].

Study Site.
The study area lies between 30 ∘ 47 29 N and 30 ∘ 54 45 N latitude and 78 ∘ 37 41 E and 78 ∘ 44 03 E longitude in the northern Himalayas, India, in the catchment of the river Bhagirathi, a tributary of the river Ganges ( Figure 1). This study area of a 12 km long road corridor with a total area of 8.88 km 2 was selected judiciously with corroboration that any landslide that occurs in the area affects the road. In the Himalayan terrain rock strength and geological structures play a major role in the landslide activity. The dominant rock types in the area include low grade metamorphic rock such as chlorite schist, schistose quartzite, and quartz mica schist along with high grade migmatites and gneisses. Rock mass properties, such as intact rock strength (IRS) computed for the area, varies between 50 and 200 MPa and corresponding cohesion of rock mass varies between 9 and 29 KPa [26]. Detailed assessment showed that the IRS varies due to compositional changes, the spacing and orientation of the joints present in these rocks, and the degree of weathering in each rock type. Elevation in the area ranges between 1550 and 2100 m with a high relative relief; average elevation of the area is around 1900 m.
The last three decades of rainfall information between 1982 and 2009 showed that the highest (1900 mm) and lowest (600 mm) annual rainfall occurred in years 2003 and 1991, respectively, with an annual average of approx. 1200 mm [27]. The area receives heavy precipitation during the summer months starting from mid of June to mid of October and moderate rainfall during the winter months from January to March. However, the rainfall is uniform in the area and no variation of rainfall was observed spatially in this road corridor. In the Himalayan region, landslides are recurring annually and are prominent during the summer months between June and October when the seasonal monsoon occurs. Landslides in this area were the result of a combination of geotectonics, adverse natural topography, such as steep slopes, weathered rocks and soils, human influences on the topography, and high rainfall [28,29].
3.2. Spatial Data Quality. Spatial data quality plays an important role for the precise landslide identification and accurate landslide mapping. Landslide data were thus collected from reliable sources. The major organizations which keep the updated record of landslides in the Indian Himalayan terrain are Border Road Organization (BRO) and Geological survey of India (GSI). The historical landslide records of BRO during 1982-2009 was used in this study for preparing the inventory. A total of 178 active landslides were mapped at the 1 : 10,000 scales. Areas affected by these landslides were clearly recognizable from the remote sensing images. They were correlated with BRO records for the road corridor occurring along the cut slopes as well as in the natural slopes of the road corridor. The mapped landslides covered an area of 0.45 km 2 , corresponding to 5.6% of the total area (min. 125 m 2 , max. 41,000 m 2 , median 1884 m 2 , and mean 3967 m 2 ). For the point pattern analysis the centroids of landslide points were marked. Each landslide was attached with the attribute, that is, the area of the landslide. The area was a proxy for the area of influence of a landslide, directly reflecting its intensity. The landslide data were converted into the multitype marked point pattern by classifying the landslides into Small and Large according to their size. Categorization was done by considering the mean area of the landslides. Therefore, for our study any landslide with size less than 4000 m 2 (mean of the distribution) was Small and more than 4000 m 2 was Large. To investigate the nature of distribution of the landslides the nearest neighbor -function was calculated for the landslide spatial patterns whether they follow a random, normal, or clustered distribution. A multitype Strauss model was fitted and the AIC was determined by applying the optimum interaction radii. To improve model fitting, various geoenvironmental factors like lithology, slope, aspect, lineament density, drainage density, weathering, soil depth, terrain units, road buffer, and land cover were considered as covariates [30]. Lithology, road buffer, terrain units, land cover, and drainage density were helpful in significantly reducing the AIC. The final susceptibility map was created by including these as explanatory variables for the landslide spatial pattern. In this study the Spatstat module was used in the R software for spatial statistical modeling [20].

Distance Curves. The temporal distribution of Large and
Small landslides in the area during 1982-2009 is shown in Figure 2(a). The area apparently contained more Small landslides in the early years than in recent years. This indicated that major landslides were subsequent to the occurrence of Small landslides as could be confirmed during the field survey. However, the occurrence of Large landslides could also be due to the regional geomorphological settings and modification of slopes in later years in this area. Figure 2(b) indicated the spatial distribution of the multitype point pattern of Small (represented by circle) and Large (represented by triangle) landslides in the road corridor resembling a clustered pattern.
The intensity was inhomogeneous throughout the road corridor as shown in Figure 3. The estimated average intensity of the landslide was 2.02 × 10 −5 landslides m −2 with a    curves for different estimates of with edge corrections were plotted well above the theoretical line indicating clustering of landslides ( Figure 4). In addition, it was clear from Figure 4 that the Small landslides were spatially correlated to large landslides within the distance of 1 km.
Next, the -functions were used for the multitype marked point patterns of landslides. Figure 5 presented a graphical interaction of the Small and Large landslides. Figure 5  type. The estimated curve for the nearest neighbor distances laid well above the curve of the Poisson process, indicating clustering. The maximum interaction radius as shown in the figure was 80 m, whereas a sharp rise of the -function occurred at distances above 30 m. These indicated more landslides in a close vicinity. Figure 5(b) revealed the observed cumulative distribution pattern of the distances from Small landslide to its nearest neighbor landslide of type Large. Here also clustered pattern was observed as all the curves are above the curve of the Poisson process. Figure 5(c) showed clustering of the Large landslides around the Small landslide within distances of 80 m. Clustering of Small around Large landslides was different from the clustering pattern of Large around Small landslides. Apparently, the occurrence of Small landslides around and following large landslides differed from that of large landslides around and following Small landslides. Large landslides were tightly clustered around Small landslides. Figure 5(d) indicated the clustered pattern of Large landslide as the plotted -function is above the curve of the Poisson process. Many large landslides however were tightly clustered within radii of 80 m. Therefore, the -function plotted in Figure 5 clearly indicated the clustering of the landslides along the road corridor.

Model Fitting.
It is a common practice to determine the interaction radii by observing -function values between the different pairs of points. There is no single optimal way to determine the interaction radii for the multitype Strauss point pattern that fits to a dataset [20]. The following matrix was determined interactively for the interaction radii that satisfy the best model fit: 19 105 ] .
The values in (13) suggested that the distance for the optimum interaction of Small and Large landslides in the study area equaled 19 m. Estimation of the interaction radii is one of the requirements for fitting multitype Strauss point models. After trying the different Strauss model for different interaction terms, the model in which the intensity was log linear function of the Cartesian coordinate was selected as the appropriate one. Since landslide is a natural phenomenon of slope failure the underlying cofactors practically follow a log linear relationship with the landslide occurrence data. The model output shown in Figures 6(a) and 6(b) represented a log linear model with a multitype Strauss process of Small and Large landslides with AIC equal to 1416 (Table 1) Next, we fitted the model for various numbers of covariates adding sequentially and monitored the AIC. Two models are significantly different if the difference of their AIC is at least 2 [31]. Various covariates were explored during model fitting in order to reduce the AIC to identify the best model. The covariates that contributed to a significant decrease were maintained. The covariates that were finally selected in the model were the causal factors for the landslides in the area. Lithology, land cover, road buffer, drainage, and terrain significantly reduced the AIC. The resultant model including these covariates along with the landslide density reflected the susceptibility of the area to landslides. The model was implemented as Fitting this model we obtained the corresponding parameter values listed in Table 1 with AIC value equal to 1395. The estimated interaction parameters for the models were 6.16, 0.56, and 1.86 for the Small-Small, Small-Large, and Large-Large interactions, respectively. Interaction parameters of the Small-Small and Large-Large landslides had values larger than 1, thus showing clustering. Small-Large landslides had an intensity value below 1 suggesting inhibition. The detail estimated parameters of the fitted intensity are presented in Table 1  influence in the final model. Further, the overall intensity was increasing because of the influence of the covariates that act as causal factors for the occurrence of landslides.

Comparison of the Model Output.
The evaluation of the model was done by comparing the Receiver operation characteristics (ROC) curves generated for the logistic regression model and the Strauss point process model. A ROC curve, also called success rate curve, was generated to analyze the success of the developed model for landslide susceptibility. By using the ROC curve, success was assessed by comparing the calculated probability values vis-à-vis their actual present condition. The area under the ROC curve (AUC) characterizes the quality of a prediction system by describing the system's ability to anticipate correctly the occurrence or nonoccurrence of predefined "events" [32]. True positive rates (sensitivity) are compared against false positive rates (1−specificity) to assess the prediction accuracy by the models. The ROC curve is shown in Figure 7.
The area under the curve (85.0% for Strauss model and 79.6% for LR model) implied a higher success rate for Strauss model, even with same sets of variables. Further, the ROC curve for the Strauss model showed a stable increase in the sensitivity as compared to the LR model finally achieving the higher accuracy. This indicated the higher sensitivity of the Strauss model to the output probability. Table 2 shows the fitted coefficients and the AIC values for the multitype/hard core Strauss point process model. As Table 1: Coefficients and AIC for the fitted models multitype Strauss model.

Intercept
Covariates AIC Lithology Land use Road buffer Drainage density Terrain units   in the multitype Straus point process model, initially the inclusion of covariates improved the model, as the AIC takes significantly lower values at each step. The last two steps, however, show that the model with three covariates is to be preferred above the model with four and five covariates. The best fitting model thus has geomorphology, drainage density, and land use as the significant covariates. As both geomorphology and land use are categorical variables, not much value can be assigned at present to the positive or negative values observed.

Discussion
Landsliding is a geomorphic slope failure process triggered by natural as well as anthropogenic factors and is controlled by unfavorable terrain conditions that act as causal factors.
To understand the landslide mechanism in an area and to identify the unknown factors affecting their occurrences, several geoenvironmental variables are included in the analysis [32]. Therefore, the problem of spatial zonation of landslides lies in the landslide inventorization, as well as in their integration with causal factors in a conceptual framework. Various statistical methods have been applied successfully to model the landslide susceptibility mapping [26,[33][34][35][36][37].
All these methods handle the landslide data as spatially independent locational variables, whereas a point pattern analysis allows preserving the level of detail offered by the landslide data itself. This is achieved through inter-landslide interactions and variation in the relative frequency of the pairs of landslides as a function of their position. A multitype Strauss point process model was used in this study with ten explanatory variables, including morphological, lithological, and structural covariates. Thefunction was able to classify landslides as a clustered pattern. The AIC of the fitted model was calculated for each of the covariates and significant ones were recorded [38]. Covariates slope, aspect, soil, weathering, and lineament density did not reduce the AIC significantly. Covariates lithology, land cover, road buffer, drainage density, and terrain units significantly reduced AIC of the fitted model and were included in the model to demonstrate the landslide susceptibility of the area. This study showed that the Strauss point process is capable of modeling the clustered pattern of landslide data for identifying causal factors for susceptibility zonation. The landslide database typically used by experts for landslide susceptibility mapping includes a comprehensively prepared landslide inventory map supported by geoenvironmental variables that cause landsliding [30,39]. The Strauss point The Strauss point process model discussed here can be applied to any set of landslide database for characterizing its inherent property for susceptibility zonation. This study presented a case for analyzing landslide data as well as identifying significant variables for generating a landslide susceptibility scenario with available information. However, the method is data driven and, therefore, the reliability of the results of modeling is always associated with the quality of the input dataset used in the model development. Thus, accuracy of the outputs invariably depends on the accuracy of the input dataset. For this study the centroid of landslide points were marked which in a way helped to reduce the positional uncertainty of landslides. The landslide data were converted into a point pattern with size (area) of landslide as mark of each landslide making it a marked point pattern. The model used the first-order (intensities) and second-order (interactions) effects of point pattern sequentially for characterizing the landslides. The required parameters like interaction radii were derived iteratively by running the model several times. After model fitting, the corresponding AIC values have been calculated for the fitted model. The model fitting with AIC showed that the covariates like lithology, land cover, road buffer, drainage, and terrain units were significant to the fitted model. All the significant covariates were combined with the landslide data in the multitype Strauss model to generate the map showing areas susceptible to landslides. The nearest neighbor distance analysis carried out for point pattern data of landslides opened an opportunity to include the spatial distribution of landslides in an alternate fashion to the existing methods. Model fitting was carried out by means of GLM functions including a combination of distance correlation functions, model fitting, and simulation which helped in describing and understanding the landsliding pattern in detail.
Several physical and terrain parameters influence the landslide process. To understand the landslide occurrence mechanism, these parameters were analyzed systematically. All the factors invariably had control on the landslides occurring along natural slopes. Small landslides occurring exclusively along the cut slopes of road corridor, however, might be controlled mainly by anthropogenic factors rather than by the natural terrain factors. Model fitting through AIC was a good way to address the sensitivity of the fitted model to the landslides data as well as to the significant covariates data. For the Strauss model being a Gibbs model, AIC acted as a goodness-of-fit of the model to the data for global sensitivity assessment of the susceptibility model.
Landslides are spatially discrete events and are controlled by number of geoenvironmental factors that are not straightforward to be easily modeled using statistical methods. Fitted models of the multitype Strauss process to landslide occurrences reflected the nature of model fit to the data. This mathematical model is the model that best fit the data; however, it is not necessarily the best model to serve in situations of practical use. Therefore, it is essential that the fitted model keeps pace with a priori knowledge for consistency. For those reasons we also fitted the multitype Strauss hardcore model. This showed slight differences with the multitype Strauss model, as the geomorphology now more stands out as the most explaining covariate, besides the -and -coordinates. The model provided a slightly better fit, that is, lower AIC values.
A point pattern analysis of the landslide data was helpful in understanding the intensity of the landsliding spatially as well as relating the interaction of the geoenvironmental variables in each landslide location with the processes generating the landslide distribution patterns. In addition, an automated module created by considering landslides as point data reduced the bias of sampling errors. In landslide studies more emphasis is given on field work and field data collection and research. Little has been done so far towards dissemination of information and facility for nonexperts. The needs of the disaster mitigating agencies working in remote hazardous front where no geomorphology expert can be made available are critical as also in case of adventure travelers or groups of pilgrims in a hilly terrain. Hence a data driven model can be useful to provide information on susceptibility of landslide in a particular region.
The factors that control landsliding in a particular region largely depend upon the geoenvironmental setup of the area. Point pattern analysis of the landslide can be helpful to urban planners in understanding the distinctive geological and geophysical characteristics at individual landslide location. These provide insight into the processes behind the landslide distribution patterns. The interaction between the landslides and their causal factors derived by means of the point process model as shown in this paper can be useful to address the management practices in hazard prone areas. It was observed from the area in this study that geoenvironmental factors related to rock types and their weathering conditions, land use patterns, and flow of water, occurring both at the surface and at the subsurface, influenced the occurrence of landslides. The marked point pattern analysis technique had the advantage of including the size of landslides into the analysis, as well as their interactions amongst themselves and with the covariates controlling the landslides. This made a unique combination to address the locations and time of occurrence of landslides in particular area leading to adopt suitable mitigation measures like slope modifications, construction of retaining walls and water channels, afforestation and alternate route alignments, and so forth.

Conclusions
A multitype Strauss point process model applied in this study successfully addressed the intensity and spatial interactions of landslides by means of the distance correlation functions, model fitting, and simulation. We concluded that this model enriches the set of statistical methods that can comprehensively analyze landslide data. It expresses properties inherently associated with each landslide and their interactions with landslides in the neighbourhood. As such, point pattern models are advantageous over GLM based models that perform a model fitting and simulation for landslide susceptibility analysis. Thus, our study adds to already existing statistical methods for landslide susceptibility mapping. Model fitting also demonstrated the significance of covariates. In this way, significant causal factors could be extracted that help to better understand the pattern of landsliding in an area.