Inventory and Distribution Characteristics of Large-Scale Landslides in Baoji City, Shaanxi Province, China

: Inventories of historical landslides play an important role in the assessment of natural hazards. In this study, we used high-resolution satellite imagery from Google Earth to interpret large landslides in Baoji city, Shaanxi Province on the southwestern edge of the Loess Plateau. Then, a comprehensive and detailed map of the landslide distribution in this area was prepared in conjunction with the historical literature, which includes 3440 landslides. On this basis, eight variables, including elevation, slope, aspect, slope position, distance to the fault, land cover, lithology and distance to the stream were selected to examine their inﬂuence on the landslides in the study area. Landslide number density (LND) and landslide area percentage (LAP) were used as evaluation indicators to analyze the spatial distribution characteristics of the landslides. The results show that most of the landslides are situated at elevations from 500 to 1400 m. The LND and LAP reach their peaks at slopes of 10–20 ◦ . Slopes facing WNW and NW directions, and middle and lower slopes are more prone to sliding with higher LND and LAP. LND and LAP show a decreasing trend as the distance to the fault or stream increases, followed by a slow rise. Landslides occur primarily in the areas covered by crops. Regarding lithology, the regions covered by the Quaternary loess and Cretaceous gravels are the main areas where landslides occur. The results would be helpful for further understanding the developmental characteristics and spatial distribution of landslides on the Loess Plateau, and also provide a support to subsequent landslide susceptibility mapping in this region.


Introduction
Landslides are a common geological hazard, causing serious casualties and economic losses every year [1,2]. Regional historical landslide inventories play an important role in the assessment of natural hazards. For example, the Oso landslide occurred in northwest Washington State on 22 March 2014, leading to devastating loss of life and destruction of property. The study of Kim et al. [3] on this event showed that a landslide had already occurred there in 2006. Compiling a complete historical landslide database is useful for analyzing the spatial distribution of landslides and their impact on the evolution of geomorphology and rivers [4]. Under the same conditions, the areas affected by previous landslides may become unstable in the future, thus moving from past records to future projections [5,6]. In many countries or regions around the world, the inventory of regional historical landslides has been completed, which is gradually maturing in terms of content, structure, accuracy and completeness. For example, the EM-DAT (International Disaster Database) was collated by KU Leuven [7]. Kir et al. [8] established a new publicly available global landslide catalogue, covering 124 countries and territories. Prior to the widespread  The northeastern part of Baoji is located in the Ordos Plate, which is mainly covered by a thick layer of Quaternary loess, usually between 0.075 and 0.25 mm in particle size [29]. Figure 1b shows the main strata and lithology of the study area, which can be divided into 12 types of rock groups according to stratigraphic ages (Table 1), including basement of Pre-Cambrian, limestone and sandstone from Cambrian to Silurian, biotite-quartz schist and granitic sandstone from Devonian to Permian, sandstone and shale from Triassic to Cretaceous, shale and marl of Tertiary, and loess of Quaternary. There are also five types of intrusive rocks: biotite granite, quartz diorite, hornblende pyroxenite (bedrock) and hornblende gneiss. Table 2 summarizes the historical earthquakes in Baoji City and neighboring areas [30], and Figure 2 shows the geographical location of the earthquake sites. In terms of climate, Baoji belongs to a semi-humid climate zone in the warm temperate zone of mid-latitudes, a continental monsoon climate type, with a mean annual temperature of 12.9°C and an average annual rainfall of~600 mm; we collected data from six meteorological stations (Figure 2) in and around the study area for the last 10 years, and then plotted the annual average temperature and annual average precipitation in Figure 3.

Visual Interpretation
In recent years, the Google Earth platform has been widely used in landslide investigations [31,32]. It contains available Landsat imagery (30 m or 15 m pan-sharpened), orthophotos (0.5-2 m), and high-resolution commercial satellite imagery (SPOT, FOR-MOSAT-2: 0.5-8 m; World View-1 and World View-2, 0.5-2.5 m) [33,34]. The satellite images are continuous and 100% coverage of the whole study area, and their resolution also meets the requirements for identifying large-scale landslides. This work performed a

Visual Interpretation
In recent years, the Google Earth platform has been widely used in landslide investi gations [31,32]. It contains available Landsat imagery (30 m or 15 m pan-sharpened), or thophotos (0.5-2 m), and high-resolution commercial satellite imagery (SPOT, FOR MOSAT-2: 0.5-8 m; World View-1 and World View-2, 0.5-2.5 m) [33,34]. The satellite im ages are continuous and 100% coverage of the whole study area, and their resolution als meets the requirements for identifying large-scale landslides. This work performed

Visual Interpretation
In recent years, the Google Earth platform has been widely used in landslide investigations [31,32]. It contains available Landsat imagery (30 m or 15 m pan-sharpened), orthophotos (0.5-2 m), and high-resolution commercial satellite imagery (SPOT, FORMOSAT-2: 0.5-8 m; World View-1 and World View-2, 0.5-2.5 m) [33,34]. The satellite images are continuous and 100% coverage of the whole study area, and their resolution also meets the requirements for identifying large-scale landslides. This work performed a visual interpretation of large landslides in the study area based on Google Earth, which was proven to be an objective and accurate approach in many previous studies [35,36].
Cruden and Varnes [37] suggested that some landslides in the past remain visible in the landscape for thousands of years after they initially moved and then stabilized. Landslides caused by earthquakes or rainfall, with the passage of time, may have been reformed to varying degrees. The vegetation restoration in the study area, however, has little influence on the outlines of historical landslides, which remain clearly visible. To identify the historical landslides, this work relied on geomorphology and morphological characteristics of the landslide body: the back wall of a landslide often shows a lap-chair shape, with obvious curved ridge lines; platform depressions may develop in the middle and rear parts of landslide body; and cracks are distributed in the middle and front edge. On the remote sensing images, they usually present tongue-shaped or circle-chair-shaped features, with irregular step-like distribution on the landslide body [38]. Figure 4 shows four typical large landslides identified in the study area. visual interpretation of large landslides in the study area based on Google Earth, which was proven to be an objective and accurate approach in many previous studies [35,36]. Cruden and Varnes [37] suggested that some landslides in the past remain visible in the landscape for thousands of years after they initially moved and then stabilized. Land slides caused by earthquakes or rainfall, with the passage of time, may have been re formed to varying degrees. The vegetation restoration in the study area, however, has little influence on the outlines of historical landslides, which remain clearly visible. To identify the historical landslides, this work relied on geomorphology and morphologica characteristics of the landslide body: the back wall of a landslide often shows a lap-chair shape, with obvious curved ridge lines; platform depressions may develop in the middle and rear parts of landslide body; and cracks are distributed in the middle and front edge On the remote sensing images, they usually present tongue-shaped or circle-chair-shaped features, with irregular step-like distribution on the landslide body [38]. Figure 4 shows four typical large landslides identified in the study area.

Inventory of Large Landslides
A total of 3440 landslides were identified in the study area, with a cumulative exten of 359.5 km 2 . From the landslide area, 12 landslides, 1021 landslides and 2407 landslides have areas ≥ 10 6 m 2 , 10 6 m 2 > area ≥ 10 5 m 2 and 10 5 m 2 > area, respectively, accounting for 0.3%, 29.7% and 70%, respectively, of the total number of landslides. The average

Inventory of Large Landslides
A total of 3440 landslides were identified in the study area, with a cumulative extent of 359.5 km 2 . From the landslide area, 12 landslides, 1021 landslides and 2407 landslides have areas ≥ 10 6 m 2 , 10 6 m 2 > area ≥ 10 5 m 2 and 10 5 m 2 > area, respectively, accounting for 0.3%, 29.7% and 70%, respectively, of the total number of landslides. The average landslide area is 104,198 m 2 , and the largest landslide area is approximately 1.71 × 10 6 m 2 , which is located in Feng County in the south of Baoji City. Figure 5 shows the geographic location of the landslide sites with different sizes; it can be seen that the number of landslides in the northwestern Long County and around the central Jintai District is the highest, and the overall size of the landslide in the northeastern corner is large. Meanwhile, we collected many typical landslides from historical articles and disaster records (Table 3), which further confirmed the accuracy of the interpretation. The landslides recorded in the articles have more or less the typical characteristics of landslides mentioned in Section 3.1 above.
ISPRS Int. J. Geo-Inf. 2021, 10, x FOR PEER REVIEW 7 of 22 landslide area is 104,198 m 2 , and the largest landslide area is approximately 1.71 × 10 6 m 2 , which is located in Feng County in the south of Baoji City. Figure 5 shows the geographic location of the landslide sites with different sizes; it can be seen that the number of landslides in the northwestern Long County and around the central Jintai District is the highest, and the overall size of the landslide in the northeastern corner is large. Meanwhile, we collected many typical landslides from historical articles and disaster records (Table  3), which further confirmed the accuracy of the interpretation. The landslides recorded in the articles have more or less the typical characteristics of landslides mentioned in Section 3.1 above.

Influencing Factors
Referring to the statistics of landslide influencing factors in previous work [44], and considering the purpose of this study, we selected eight influencing factors of landslides,

Influencing Factors
Referring to the statistics of landslide influencing factors in previous work [44], and considering the purpose of this study, we selected eight influencing factors of landslides, including elevation, slope, aspect, slope position, distance to the fault, land cover, lithology and distance to the stream. The DEM data used were ALOS PALSAR DEM (resolution 12.5 m), from which the stream, slope and aspect information were extracted. The relative position index (RPI) was introduced to define the slope position [45]. The RPI of a landslide on the slope is defined as the ratio of the shortest distance from its centroid to the valley and the sum of the shortest distances from the point to the valley and ridge (Figure 6a), which was divided into 5 levels ( Figure 6b). The fault and lithology data were derived from Deng et al. [46] and a 1:200,000 geological map, separately. The land cover data (resolution 12.5 m, after resampling) were from European Space Agency (ESA). including elevation, slope, aspect, slope position, distance to the fault, land cover, lithology and distance to the stream. The DEM data used were ALOS PALSAR DEM (resolution 12.5 m), from which the stream, slope and aspect information were extracted. The relative position index (RPI) was introduced to define the slope position [45]. The RPI of a landslide on the slope is defined as the ratio of the shortest distance from its centroid to the valley and the sum of the shortest distances from the point to the valley and ridge ( Figure  6a), which was divided into 5 levels ( Figure 6b). The fault and lithology data were derived from Deng et al. [46] and a 1:200,000 geological map, separately. The land cover data (resolution 12.5 m, after resampling) were from European Space Agency (ESA).

Statistics and Analysis
Besides the RPI, the location of a landslide is determined in more detail by using its elevation, distance from the river (Figures 7 and 8), slope and aspect (Figures 9 and 10). Zhuang et al. [24] found that the overall distribution of landslides in the Loess Plateau has a certain clustering characteristic, i.e., landslides may appear in clusters or bands. In the process of spatial analysis, a similar phenomenon was noted, which was inferred that the zonal distribution of landslides is influenced by faults (Figure 11), and the aggregated distribution may be influenced by land-use types ( Figure 12).

Statistics and Analysis
Besides the RPI, the location of a landslide is determined in more detail by using its elevation, distance from the river (Figures 7 and 8), slope and aspect (Figures 9 and 10). Zhuang et al. [24] found that the overall distribution of landslides in the Loess Plateau has a certain clustering characteristic, i.e., landslides may appear in clusters or bands. In the process of spatial analysis, a similar phenomenon was noted, which was inferred that the zonal distribution of landslides is influenced by faults (Figure 11), and the aggregated distribution may be influenced by land-use types ( Figure 12). including elevation, slope, aspect, slope position, distance to the fault, land cover, lithology and distance to the stream. The DEM data used were ALOS PALSAR DEM (resolution 12.5 m), from which the stream, slope and aspect information were extracted. The relative position index (RPI) was introduced to define the slope position [45]. The RPI of a landslide on the slope is defined as the ratio of the shortest distance from its centroid to the valley and the sum of the shortest distances from the point to the valley and ridge ( Figure  6a), which was divided into 5 levels (Figure 6b). The fault and lithology data were derived from Deng et al. [46] and a 1:200,000 geological map, separately. The land cover data (resolution 12.5 m, after resampling) were from European Space Agency (ESA).

Statistics and Analysis
Besides the RPI, the location of a landslide is determined in more detail by using its elevation, distance from the river (Figures 7 and 8), slope and aspect (Figures 9 and 10). Zhuang et al. [24] found that the overall distribution of landslides in the Loess Plateau has a certain clustering characteristic, i.e., landslides may appear in clusters or bands. In the process of spatial analysis, a similar phenomenon was noted, which was inferred that the zonal distribution of landslides is influenced by faults (Figure 11), and the aggregated distribution may be influenced by land-use types ( Figure 12).          Referring to previous work [47], we chose the landslide number density (LND) and landslide area percentage (LAP) as two indicators of landslide distribution to measure the spatial distribution of landslides in the study area. Of them, LND describes the concentration of landslides, defined as landslide number per km 2 (Equation (1)); LAP means the scale of landslides, which is defined as the percentage of landslide area in each factor category (Equation (2)). We can gain initial insight into the subsequent susceptibility map [48] from the statistical results. The extraction of all attributes was done in ArcGIS. The results show that the LND and LAP of the whole study area are 0.19/km 2 and 1.98%, respectively, in the study area, the LND distribution map is shown in Figure 13. Table 4 shows the results of LND and LAP for all variables.

=
(1) Referring to previous work [47], we chose the landslide number density (LND) and landslide area percentage (LAP) as two indicators of landslide distribution to measure the spatial distribution of landslides in the study area. Of them, LND describes the concentration of landslides, defined as landslide number per km 2 (Equation (1)); LAP means the scale of landslides, which is defined as the percentage of landslide area in each factor category (Equation (2)). We can gain initial insight into the subsequent susceptibility map [48] from the statistical results. The extraction of all attributes was done in ArcGIS. The results show that the LND and LAP of the whole study area are 0.19/km 2 and 1.98%, respectively, in the study area, the LND distribution map is shown in Figure 13. Table 4 shows the results of LND and LAP for all variables.

LND =
Landslide number Classi f ied area o f variable I (1) scale of landslides, which is defined as the percentage of landslide area in each factor egory (Equation (2)). We can gain initial insight into the subsequent susceptibility [48] from the statistical results. The extraction of all attributes was done in ArcGIS. results show that the LND and LAP of the whole study area are 0.19/km 2 and 1.98% spectively, in the study area, the LND distribution map is shown in Figure 13. Tab shows the results of LND and LAP for all variables. = = Figure 13. The LND distribution map derived from this study.

Overall Distribution Characteristics of Large Landslides
The distribution of large landslides in Baoji presents two characteristics: unevenness (or a strong spatially non-uniform pattern) and aggregation. The former is reflected in that the number of the landslides is much higher in the north than in the south throughout the study area, while there is no such a feature in the horizontal direction. Landslides are sparsely distributed in the southeast of this area, with several dense landslide areas in the northwestern part. In areas where a large number of landslides occurred, there is also an obvious zonal and aggregated distribution of slope failures, of which the former (Figure 14a) is located near the Longxian-Baoji fault in the Jintai District on the north bank of the Wei River, and the latter (Figure 15a) in the affected area of the Longxian earthquake ( Figure 2).

Correlation between Landslides and Influencing Factors
This work analyzed the relationship between landslides and influence factors such as topography and geology. The statistical results from this analysis are shown in Figure  16.

Correlation between Landslides and Influencing Factors
This work analyzed the relationship between landslides and influence factors such as topography and geology. The statistical results from this analysis are shown in Figure  16.

Correlation between Landslides and Influencing Factors
This work analyzed the relationship between landslides and influence factors such as topography and geology. The statistical results from this analysis are shown in Figure 16.

Topographic Factors
The elevation of the study area ranges from 404 to 3719 m, divided at 100 m intervals in this work. Statistics show that most of the landslides are located at elevations from 500 to 1400 m, accounting for 94.2% of the total. In the range of elevations greater than 1800 m, there are only 45 landslides, although the land area occupies 16.7% of the study area. The peak values of LND and LAP are 0.37/km 2 and 3.9% at 600-700 m and 1100-1200 m, respectively (Figure 16a).
Slope means the degree of steepness of the surface unit, which determines the exposed face when a landslide occurs to some extent. The slope of the study area ranges from 0.8 to 62.1 • and was divided into 10 • intervals in this work. The area of the interval less than 40 • occupies 94.6% of the total, among which the number of landslide development in the 10-20 • interval is the largest, 1555, where the LND and LAP also reach their peak at 0.35/km 2 and 3.8%, respectively (Figure 16e). The statistics show no case that the higher the slope is, the more likely landslides are to occur.
The slope aspect (facing direction) has a great effect on the mountain ecology, thus affecting the occurrence of landslides. The area of 16 aspect divisions in the study area is approximately equal. The peaks of both LND and LAP appear in the WNW direction with 0.29/km 2 and 2.7%, respectively. In general, the LAP values of west-facing slopes are significantly larger than those of east-facing slopes, and the average landslide size is also slightly larger than that of east-facing slopes.
Slope position links to the landform feature may also affect the occurrence of landslides. There is a positive correlation between the area of landslides and the area of slope location division in the study area, generally showing a trend of first increasing and then decreasing. There are 1716 landslides on the middle slope, with LND and LAP of 0.34/km 2 and 3.5%, respectively.
The influence of land-use types on the occurrence of landslides is also obvious. Four main land types exist in the study area: trees, crops, shrub and built area. Among them, the area covered by trees occupies 62.6% of the total, with 1153 landslides. There are 1383 landslides in the areas covered by crops, where both LND and LAP reach the peaks of 0.39/km 2 and 4.8%, respectively, implying landslide-prone areas (Figure 16c).

Geological Factors
In the study district, the area of loess-dominated Quaternary and granite-dominated intrusive rocks is the largest, occupying 60.8% of the total, where the landslides mainly developed in loess sediments, with 1995 landslides identified. From the Cretaceous to the Quaternary strata, large LND and LAP values are present, with peaks of 0.42/km 2 and 3.7% in the Cretaceous, indicating a major influence of loess and sandy conglomerates on the occurrence of landslides.
Active faults are closely related to genesis of earthquakes, thus affecting seismic landslides. Figure 16h shows that faults have obvious influence on the occurrence of landslides. There are 3046 landslides developed within 20 km from the fault, accounting for 88.3% of the total, and the LND shows a decreasing trend as the distance to the fault increases. The peaks of both LND and LAP appear at the closest location to the fault, with 0.34/km 2 and 2.9%, respectively. Although LND and LAP fluctuate somewhat at 20 km away from the fault, the total number of landslides occupies only 11.7%, meaning that the degree of influence is relatively weak.

Hydrological Factors
The study area hosts several tributaries of the Wei River that is a major stream in Shaanxi. There are 531 landslides developed within 400 m from the stream. LND and LAP show a gradual decrease in the range of 200-1200 m to the stream. LND reaches the maximum value of 0.33/km 2 in the distance range of 200-400 m to the stream, and LAP reaches the maximum value of 3.2% in the range of 600-800 m (Figure 16d).

Landslide Inventory
For the identification of large landslides in Baoji, we mainly referred to the existing standards for manual interpretation of remote sensing images [4,49,50]. The methods used in this study include the following: (1) Based on the basic outline of typical landslides, including steep walls at the back edge of landslides, sidewalls, undulating topography at the front edge, irregular stepped distribution on landslides, and the formation of weirs by landslide blockage [51], we established the interpretation marks of large landslides in the study area. (2) Using multi-period high-resolution satellite imagery from Google Earth, acquired from 2000 to 2021, we identified key parts of landslides based on interpretation experiences, and used polygons to circle the spatial extent of single landslides. (3) By comparing the landslide locations recorded in the existing literature (Table 3), the landslides we interpreted earlier were further checked so as to produce a relatively complete database of large landslides. Considering that the force source mechanisms that generate landslides include pushing and traction [52], we used the centroid of landslides (Figure 6a) to count the spatial distribution pattern of landslides. Xu et al. [53] adopted the same method to count the co-earthquake landslides by the 2008 Wenchuan earthquake, indicating that the results are objective and accurate.

Spatial Distribution of Large Landslides
As Xu et al. [54] mentioned, there is some clustering of landslides in the study area, but our results do not conform to the characteristic of dense landslides in the south and gradual sparseness toward the north. Le et al. [55] suggested that the basin edge, river bank, terrace edge, fault zone, and town surroundings tend to be zones of aggregation of landslides in loose soils, which can be seen in Figures 6 and 11. The intense neo-tectonic movements and accompanying paleoclimatic environmental changes are the primary driving force for the formation of river systems in the Loess Plateau, where rivers flow from northwest to southeast along the terrain [56]. The Wei River is the main stream in the study area, with numerous tributaries (Figure 8). Landslides are concentrated on both sides of this river valley and the intersection of several rivers, where the LND is as high as 0.33/km 2 within 200 to 400 m from the stream, which can be explained by differential lifting and lowering movements of the crust in the Loess Plateau after the Tertiary time that led to the instability of slopes caused by erosion [57]. In addition, earthquakes as a typical representative of the neotectonic movement, might have resulted in the "flaky" dense distribution of landslides in loose soils, which is also reflected in this study. For instance, an earthquake with Ms6.0 in 1704 was recorded in Long County, northwest corner of Baoji, where many old landslides are present around the earthquake area (Figure 14a). With reference Peng et al. [24], the study area can also be zoned according to the location of large landslides, and such subsequent work needs to be done based on further collection of detailed geological data.

Correlation between Landslides and Influencing Factors
The geological structure of Baoji is complex, with many fault zones in the territory, such as the Qinling North Edge fault zone. The Wei River is the main stream in the study area, with numerous tributaries, and the land around the twelve districts and counties in the territory has been greatly reformed ( Figure 12).
The study area is dominated by low-middle elevations; most landslides occurred in the regions with elevations of 500 m to 1400 m, accounting for 94.2% of the total, which is in good agreement with Hu [58]. The reason may be that these areas have accumulated many loose deposits, which are easy to slide under external disturbance. From this work, the peak value of LND is 0.37/km 2 at elevations 600-700 m, likely due to the effect of streams there.
Every indicator reaches their peaks in the range of slope 10 • to 20 • , and the number of landslides shows a decreasing trend as the slope increases, consistent with the conclusion obtained by Xu et al. [59].
For the aspect (Figure 10), the largest number of landslides is in the WNW and NW directions, with the peaks of both LND and LAP occurring in the WNW direction. We note that the orientation of several faults within the northern margin of the Qinling North Edge fault zone is close to the dominant sliding direction of the landslides, which leads to the conclusion that the slope aspect has a strong controlling on the occurrence of landslides. Qiu et al. [60] obtained similar results by coupling the relationship between the occurrence of landslides in loose soils and slope aspect. They also found that the aspect effect varied with the type of landforms: the more mature the development of geomorphology, the more obvious the slope aspect effect on the landslides. Hu [58] analyzed 6627 landslides on the Loess Plateau; however, he found that there is not much correlation between landslides and the aspect, probably because he did not ensure the completeness of the landslides in the area involved. Therefore, it is necessary to construct a complete regional landslide database.
With the position of a slope from low to high, LND and LAP show a trend of first increasing and then decreasing; meanwhile, the number of landslides reaches the maximum in the middle slope, which is likely the result of the joint influence of the elevation and river (Figure 17b). The Quaternary strata are the largest area exposed, occupying 38.4% of the total area of the study area, mainly distributed in the plains and basins on the north bank of the Wei River, followed by granite-based intrusive rocks occupying 25.2% of the area. The number of landslides in the Quaternary loess-covered strata is the largest, with 1995, consistent with Xu et al. [23]. Every indicator reaches their peaks in the range of slope 10° to 20°, and the number of landslides shows a decreasing trend as the slope increases, consistent with the conclusion obtained by Xu et al. [59].
For the aspect (Figure 10), the largest number of landslides is in the WNW and NW directions, with the peaks of both LND and LAP occurring in the WNW direction. We note that the orientation of several faults within the northern margin of the Qinling North Edge fault zone is close to the dominant sliding direction of the landslides, which leads to the conclusion that the slope aspect has a strong controlling on the occurrence of landslides. Qiu et al. [60] obtained similar results by coupling the relationship between the occurrence of landslides in loose soils and slope aspect. They also found that the aspect effect varied with the type of landforms: the more mature the development of geomorphology, the more obvious the slope aspect effect on the landslides. Hu [58] analyzed 6627 landslides on the Loess Plateau; however, he found that there is not much correlation between landslides and the aspect, probably because he did not ensure the completeness of the landslides in the area involved. Therefore, it is necessary to construct a complete regional landslide database.
With the position of a slope from low to high, LND and LAP show a trend of first increasing and then decreasing; meanwhile, the number of landslides reaches the maximum in the middle slope, which is likely the result of the joint influence of the elevation and river (Figure 17b). The Quaternary strata are the largest area exposed, occupying 38.4% of the total area of the study area, mainly distributed in the plains and basins on the north bank of the Wei River, followed by granite-based intrusive rocks occupying 25.2% of the area. The number of landslides in the Quaternary loess-covered strata is the largest, with 1995, consistent with Xu et al. [23]. From the Cretaceous to the Quaternary, large LND and LAP values occur in each stratum for the following reasons: (1) Sandstone within the Cretaceous strata, which was the product of alluvial fan and fluvial deposits, is unstable. (2) Shales and marls in Tertiary have small strength, prone to softening and swelling after soaking in water. (3) The loess layer in Quaternary is loose and permeable, and the underlying bedrock is relatively watertight, so the water is collected in the weak rock layer, softening the beds and forming a sliding zone, most likely to trigger a "dual structure" (loess-bedrock) landslide [61].
Active faults have obvious influence on the occurrence of landslides: 3046 landslides developed within 20 km from the faults, occupying 88.3% of the total, and the LND shows a decreasing trend as the distance to the faults increases. Peng et al. [62] analyzed the Liujiabao landslide, which occurred under the influence of loess-mudstone stratigraphy, faults, seismic forces and heavy rainfall. Moreover, their extensive investigations found From the Cretaceous to the Quaternary, large LND and LAP values occur in each stratum for the following reasons: (1) Sandstone within the Cretaceous strata, which was the product of alluvial fan and fluvial deposits, is unstable. (2) Shales and marls in Tertiary have small strength, prone to softening and swelling after soaking in water. (3) The loess layer in Quaternary is loose and permeable, and the underlying bedrock is relatively watertight, so the water is collected in the weak rock layer, softening the beds and forming a sliding zone, most likely to trigger a "dual structure" (loess-bedrock) landslide [61].
Active faults have obvious influence on the occurrence of landslides: 3046 landslides developed within 20 km from the faults, occupying 88.3% of the total, and the LND shows a decreasing trend as the distance to the faults increases. Peng et al. [62] analyzed the Liujiabao landslide, which occurred under the influence of loess-mudstone stratigraphy, faults, seismic forces and heavy rainfall. Moreover, their extensive investigations found that such a geological background and climatic condition are common in northwest China, which provide favorable conditions for large-scale landslides in this region.
The study area is dominated by crops in addition to the area covered by trees, which occupied 30.2% of the total area, and the building area accounts for 9.7% of the total. There are 1383 landslides developed in the area covered by crops, where both LND and LAP reach their maximum values. Figure 12 shows the location where landslides are concentrated, indicating that landslides are directly related to human activities, especially agricultural irrigation. The residents of Loess Plateau have long been engaged in agricultural irrigation for, which poses an influence on the groundwater in the loess plateau area, resulting in a continuous rise of the water level, which gradually deforms the soil at the foot of the slope and leads to the instability of the slope [63]. The most typical affected area of such landslides is the Hefangtai area in Lanzhou City [64].
The number of landslides is mostly within the distance of 200 to 400 m to the stream, which tends to decrease as the distance increases. As the largest tributary of the Yellow River, the Wei River has well-developed drainage that has created a heavily eroded valley landscape in the study area. The sand content of the Wei River is gradually decreasing in the study area, for which the reasons are as follows: after the Wei River entered the Guanzhong Plain from the Baoji Gorge, the valley suddenly released, and the bed slope decreased greatly, which led to the decrease in the sediment-carrying capacity of the water flow. This further intensified the scouring on both sides of the river, leading to changes in the composition and structure of the slope body and causing landslides.

Comparison with Existing Landslide Databases
The study area, Baoji, is located at the southwest edge of the Loess Plateau, and we attempted to make some qualitative and quantitative comparative analyses based on this paper and the published literature involving the spatial distribution of landslides in Baoji. Firstly, Peng et al. [24] mapped the distribution of landslides on the Loess Plateau, including 14,544 landslides, and divided the Loess Plateau into eight landslide zones. In his study, the maximum LND, 0.18/km 2 , was obtained in the Haiyuan seismic area, compared to the average LND of 0.19/km 2 in our study; thus, it seemed that although each landslide mentioned by Peng possessed more detailed geological information, his database was clearly incomplete. A similar study was conducted by Hu [58], but he only selected 6627 typical landslides of the Loess Plateau for mapping landslide distribution, and the LND of Baoji was 0.0141/km 2 , which is far lower than the results of our study. In addition, although Xu et al. [23] compiled a database of landslides on the Loess Plateau containing more than 80,000 landslides, the regions he covered were mainly distributed in the four extremely seismic areas of Hongdong, Huaxian, Tongwei, and Haiyuan, and the area affected by the 2010 Tianshui rainstorm event. His study did not include the distribution of landslides within the entire city of Baoji; the 4122 landslides in Shaanxi Province mentioned by Zhuang et al. [24] were similar.
Secondly, the literature with Baoji as the study area was of great help to this study and is listed as follows: The 392 large loess landslides in Baoji mentioned by Xu et al. [29], with an average area of 47.18 million m 2 and the approximate geographical locations of the landslides, could be obtained from his paper, which provided a good reference for our interpretation as well as cross-checking. Shi et al. [65] analyzed the causes of more than 70 landslides in a small area of Jintai District, Baoji city; Figure 13a shows a part of his study. Xu et al. [39] analyzed the distribution pattern of the geological hazard in the area from Long County to Qianyang in the northwest of Baoji city, and proposed five concentrated distribution areas and six concentrated distribution zones. The detailed locations of the hazards were in good agreement with our interpretation, and the LND was 0.181/km 2 , which was close to our study. It could be seen that the database of geological hazards in small areas built by Xu et al. [39] was very complete. Most of the remaining literature deals only with the detailed analysis of single landslides within Baoji, as shown in Table 3.

Conclusions
Based on Google Earth and remote sensing imagery, this work prepared an inventory of large-scale landslides in the Baoji city, Shaanxi Province, China. In total, 3440 landslides were identified, of which the largest individual is 1,709,700 m 2 , the smallest is 634.5 m 2 , and the average area is 104,198 m 2 . We chose the landslide number density (LND) and landslide area percentage (LAP) as evaluation indicators and conducted statistical analysis on the correlation between the landslides and influencing factors. The results show that LND and LAP are 0.19/km 2 and 1.98%, respectively. The landslides are more likely to occur in the elevation range of 500-1400 m. The LND and LAP reach their peaks at slopes of 10-20 • . Slopes facing WNW and NW are close to most of the orientation of faults, which are the dominant area of landslides. Middle and lower slopes exhibit high landslide susceptibility. LND and LAP show a general trend of decreasing with increasing distance to the fault and stream, followed by a slow rebound. There are more landslides in the areas covered by crops, compared to other regions. For lithology, areas covered by Quaternary loess and Cretaceous gravels are the main areas where landslides occurred. Compared with relevant previous work, the database of large-scale landslides within Baoji obtained in this study is more complete. However, our current work is very superficial, and we need to conduct more in-depth research in some aspects. We believe that the following aspects should be the focus of the next research, which are subdivided as follows: (1) the landslide susceptibility mapping in Baoji; (2) to conduct more detailed investigations, including the types of landslides [66], the morphology of the landslide bodies, etc.; and (3) to analyze the occurrence mechanism and triggering factors of typical landslides based on the previous information. Overall, our results provide necessary support for various subsequent studies.