Geological Hazard Assessment of Secondary Collapses Due to Volcanic Earthquakes on Changbai Mountain in China

: In recent years, the volcanic activity of Changbai Mountain has been accompanied by several earthquakes, and the frequent human engineering activities have led to a gradual increase in the number of collapses in the region, which severely impacts residents’ lives and property safety. In northeastern China, the Changbai Mountain area in the southeastern Jilin Province is a typical mountain environment. This paper selects 12 evaluation indicators to build a hazard assessment system, including slope, aspect, elevation, curvature, lithology, NDVI, land use type, distance from the fault, the river from the road, volcanic earthquake, and annual average precipitation. Using emotional weight (G1 method) and objective weight (WOE-CV method), the hazard due to collapses in the study area is evaluated too. Among them, the transcendence probability of volcanic earthquakes greater than VI degree represents the relationship between Changbai Mountain volcanic earthquakes and the assessment of geological collapse hazard. The results show that high-and very high-hazard areas are mainly distributed in densely populated areas and national and provincial trunk lines, with apparent spatial agglomeration characteristics. The low-hazard area, medium-hazard area, high-hazard area and very high-hazard area accounted for 19.33%, 44.19%, 33.85% and 2.63% of the total area of the study area, respectively. By comparing the previous geological hazard survey reports in the area with the collapse hazard zoning map in this paper, 87.72% of the known collapse hazard areas are distributed within high and very high hazard zones, indicating that the conclusions of the article are more accurate and in line with the actual situation. Results from collapse-related hazards can provide relevant guidance for preventing and controlling geological risks posed by volcanic earthquakes aﬀecting Changbai Mountain.


Introduction
In the past few years, with the continuous development of cities and the increasing number of human activities, the disruption of ecological environments has led to the increasing frequency of geological related hazards, and the economic losses caused by hazards are also increasing [1,2]. As the world's second largest economy, China's economic growth has always been in the forefront of the world. How to reduce the losses caused by disasters while maintaining economic development is an important part of disaster prevention and reduction, Therefore, it is of great significance to carry out hazard assessments concerning geological hazards [3].
Changbai Mountain is a famous tourist a raction and an important part of the region's economic development. With a nearby resident population in the tens of millions, in the event of eruptive activity from the Changbai Mountain Volcano, the resulting secondary geological hazards could cause huge economic losses and human casualties.
From 2002 to 2005, the monitoring of Changbai Mountain Tianchi V showed a significant increase in earthquake magnitude, topographic variation, frequency and greenhouse gas emissions associated with the mountain [4,5]. In May 2010, the elevation of the northern slope of the Changbai Mountain Tianchi volcanic cone reversed and decreased by 12.72 mm, breaking the previous pa ern of an annual increase of 4 mm. Based on the above monitoring results, the Changbai Mountain Tianchi volcano may be entering an active phase. Therefore, it is necessary to evaluate the hazard of secondary collapse geological hazards caused by volcanic activity in Tianchi, Changbai Mountain.
The development characteristics, distribution areas and inducing factors of different types of geological hazard are different, so the selection of evaluation indicators should equally differ [6,7]. The study area is located in the southeast of the Jilin Province (NE China), and its population is mainly distributed near cities and major transportation roads, while other areas are mostly forested and sparsely populated. According to the analysis of the previous field investigation data, the number of hazards has a positive correlation with human engineering activities. At the same time, the study area is close to the Changbai Mountain crater. Through the analysis of seismic monitoring data acquired over the last 20 years, the frequency of Changbai Mountain volcanic earthquakes has been on the rise, which has a positive impact on the collapse geological hazard in the study area. Therefore, this article selects indicators such as slope, aspect, elevation, curvature, lithology, annual average rainfall, NDVI, distance from faults, and distance from rivers. Additionally, indicators related to human engineering activities such as distance from road and land use type were selected [8]. At the same time, due to the study area being located near the crater of Changbai Mountain, volcanic earthquakes are selected as the inducing factors to assess the collapse geological hazard caused by volcanic earthquakes in the study area.
Volcanic earthquakes and their induced secondary geological disasters not only have a certain impact on ecosystems, but also cause damage to human life and property safety. Many domestic and foreign scholars have been commi ed to this research. Italian scholar De Vita Sandro [9] found that volcanic activity is an undeniable triggering factor in his research on landslide hazard. American scholar Keefer [10,11] analyzed the data of collapse and landslide hazard caused by earthquakes worldwide, and obtained the distribution pa erns of earthquake parameters and collapse and landslide hazard. Japanese scholar Tamura [12,13] analyzed case data of earthquake-induced landslide hazards over the past 100 years in 1978, and obtained the relationship between earthquake parameters such as magnitude and epicenter distance and landslide boundary distance. Huang Yidan [14], a Chinese scholar, and others found that with the increase in earthquake intensity, the development density of collapse geological disasters also increased by interpreting the remote sensing images of the earthquake area and combining with the field survey data. Zhang Jianqiang, Fan Jianrong [15] and others conducted risk analyses on earthquake-induced landslides in Wenchuan County, Sichuan Province by selecting lithology, slope, distance from fault and other indicators.
Previous evaluation methods mostly resort to a single evaluation model. In order to improve the accuracy of the evaluation results of collapse geological hazards in the study area, this article first uses the objective weight method [16,17] to analyze known data, and then combines the subjective weight method [18] to conduct hazard assessment of collapse geological hazards in the study area. At the same time, in the selection of evaluation indicators, comprehensive consideration was given to the terrain and geomorphic characteristics, human engineering activities, special triggering factors, and other aspects of the study area. Indicators such as slope, aspect, elevation, curvature, lithology, annual average precipitation, NDVI, distance from faults, and distance from rivers were selected. At the same time, indicators related to human engineering activities such as distance from roads and land use type were also selected. Combined with the special triggering factor of volcanic earthquakes, the hazard assessment of the secondary geological hazard due to vol-canic earthquakes was carried out. The receiver operating characteristic curve (ROC) [19] is used to compare and evaluate the accuracy of the model before and after the combination.

Study Area
Changbai Mountain is located in the southeast of Jilin Province, NE China, bordering North Korea to the south and Russia to the east. The study area includes the Changbai Mountain Scenic Area and three surrounding counties, namely, Fusong County, Antu County and Changbai County (Figure 1), with an area of 16,108.96 km 2 . Figure 2 shows that during the field investigation, most of the collapse geological hazard points recorded are distributed along the highway.
The study area belongs to the continental humid monsoon climate in the northern temperate zone. Summer is warm and rainy with concentrated precipitation, while winter is cold and long. The annual average temperature is 2-4 °C, and the annual average rainfall is 630~800 mm. The drainage system in the study area is relatively organized, and drainage density is large, with the Yalu, Songhua and Tumen Rivers system. According to the genesis, the geomorphological types can be divided into three categories: volcanic landforms, tectonic denudation landforms and erosion accumulation landforms. The geological structures such as fault zone and fold are mainly distributed in the northeast and southwest of the study area.
The population is mainly distributed in the south of Changbai County, the northeast of Antu County and the central and western parts of Fusong County. At the same time, due to the expansion of the city and the improvement of corresponding infrastructure, more and more people gather in the city and surrounding urban areas, while the population in rural and remote areas is becoming more and more sparse, resulting in a relatively concentrated population distribution. On the other hand, the crater of Changbai Mountain is a volcanic landform with a high elevation, and most of the nearby areas are primeval forests with inconvenient transportation, which are not suitable for human habitation and have a lower population distribution.

Methods and Data
This study establishes, analyzes and validates a collapse hazard evaluation model, through the steps of field investigation, data collection, factor correlation analysis and evaluation index selection. The order relation analysis method (G1) and weight of evidence and coefficient of variance (WOE-CV) methods are used to evaluate the collapse geological hazard. Among them, the G1 method is a subjective weight evaluation model improved by the analytic hierarchy process (AHP) [20], which can effectively negate the difference between evaluation indicators. Weight of evidence (WOE) and coefficient of variation (CV) are objective evaluation models. The WOE-CV combination model can optimize the accuracy of the evaluation index weight and improve the accuracy of the evaluation results. The specific process is shown in Figure 3.

Index Selection
Collapse phenomena are caused by the joint action of various factors [21]. These include landform, formation lithology and hydrogeology, etc., while the induced factors include volcanic earthquake, rainfall, human engineering activities, etc. Based on these two aspects, according to the characteristics of the study area, this paper selects 12 indicators including NDVI, elevation, slope, curvature, lithology, slope, distance from fault, distance from river, distance from road, annual average precipitation, land use type, and volcanic earthquake to build the collapse geological hazard assessment system ( Figure 4).
In order to analyze the influence degree of Changbai Mountain volcanic earthquakes on the collapse geological hazard in the study area, this paper uses the transcendental probability of the volcanic earthquake intensity greater than VI degrees to express the impact of Changbai Mountain volcanic earthquakes on the stability of the slope in the study area [22]. The specific calculation process is as follows: Gutenberg Richter recurrence relationship- where M is the magnitude, N (M) is the cumulative frequency of earthquake occurrence, and a and b are constants. The relationship between the frequency of volcanic activity and the magnitude of earthquakes associated with volcanic activity [23] lgN(VEL) = −0.2286H + 1.6286 (2) where H is the magnitude of the earthquake associated with volcanic activity in Changbai Mountain, and N (VEL) is the frequency of volcanic activity. Through the analysis of the statistical results of earthquake-induced collapse phenomena in various regions, it is found that a seismic intensity of VI is the minimum needed to trigger a collapse. Therefore, the seismic intensity a enuation formula in the study area can be expressed as follows [24]: where L is the earthquake intensity, T is the magnitude at the time of an earthquake, and G is the epicentral distance. In order to reflect the impact of volcanic earthquakes on the collapse phenomena in the study area, this paper uses the probability of a Changbai Mountain volcanic earthquake of an intensity exceeding VI, and the research time span is 5000 years. The specific formula is as follows: where K represents the magnitude of Changbai Mountain volcanic earthquake required to cause an intensity of VI in the study area, and is the probability that the intensity exceeds VI.
See Table 1 for the selection principles of evaluation indicators. Table 1. Selection principles of evaluation indicators.

Criteria Evaluation Indicators Principle of Selection
Topographic features

Elevation
The vegetation coverage, sunlight exposure area and water catchment of slopes at different elevations, with varying impacts on the rock and soil mass of slopes [25]. Slope When the slope exceeds a certain value, it will intensify the gravity effect, easily causing the sliding and collapse of rock layers and soil, especially under the influence of natural disasters such as rainy seasons or earthquakes. Therefore, slope is one of the important factors in the occurrence of collapse geological hazard [26].

Curvature
The value of the curvature reflects the shape of the slope, including concave slopes, convex slopes, straight slopes, etc. Convex slopes will form a phenomenon of stress concentration at the slope corner, resulting in a decrease in the stability of the slope. The middle part of the concave slope is prone to strong erosion, which is also not conducive to the stability of the slope [27].
Aspect Different slope directions are exposed to sunlight and precipitations for different time ranges, which further affects the distribution of vegetation on the slope, the weathering degree of rock and soil, and then the occurrence of geological hazard [28].

Rock and soil mass Lithology
Formation lithology is the material condition for forming geological hazards. Different formation lithologies have different shear and weathering resistance of rock mass. According to the analysis and statistics of the formation lithology in the study area, geological hazards of collapse mainly occur in hard rock layers [29].
Geological structure Distance from fault The strength of the rock and soil around the fault zone is low, and the weathering resistance is weak, which can lead to unstable joint surfaces, and then induce the occurrence of collapse geological hazards [30].

Meteorology and hydrology Distance from river
Drainage affects the stability of the slope, mainly manifested in scouring and eroding the slope toe, changing the original slope gradient and reducing the strength of the slope [31].
Human engineering activities Distance from road Through statistical analysis, it is found that more than 80% of the collapse geological hazard in the study area is affected by human engineering activities. Among them, the construction of roads is the primary factor affecting the stability of slopes, with the excavation, backfilling and stacking of materials during the construction process [32].
Vegetation cover

NDVI
The vegetation normalization index can significantly reflect the spatial distribution of vegetation in the study area. The smaller the NDVI value is, the sparser the vegetation is. The higher the NDVI value is, the higher the vegetation coverage is. Areas with high vegetation cover are less susceptible to erosion than bare areas [33].

Land use type
The type of land use can objectively reflect the natural environment and human engineering activities in the study area, both of which will affect the occurrence of geological hazards to a certain extent [34,35].

Inducing factor
Volcanic earthquake The study area is close to the crater of Changbai Mountain, and volcanic activity has occurred frequently in recent years, so volcanic earthquake is selected as the inducing factor for consideration.

Average annual precipitation
Rainfall may induce slope instability and deformation. Part of the rainfall will form surface runoff and scour the slope surface. At the same time, part of the rain will penetrate into the slope, changing the stress balance inside the slope, and further inducing the occurrence of collapse geological hazards [36]. The research data in this paper mainly consist of spatial vector data, remote sensing image data, and meteorological and hydrological data. The data sources of each evaluation index are shown in Table 2. The elevation, curvature, slope and aspect are extracted from DEM grid data using ArcGIS software. After obtaining the above data, the ArcGIS software is used for projection transformation and resampling, so that the data can reach a unified resolution and spatial coordinate system. Other information comes from the Jilin Provincial Geological Environment Monitoring Station.

Building an Evaluation Model
According to the field survey and analysis of historical data in the study area, most of the collapse sites are mainly concentrated in densely populated villages, towns and roads. The reason is that the expansion of towns and the construction of related ancillary facilities have changed the original terrain characteristics. Human factors are an im-portant reason for the occurrence of collapses in the study area. In this paper, the subjective weight method (G1 method) and the objective weight method (WOE-CV) are combined to evaluate the collapse geological hazard in the study area based on remote sensing satellite images, field survey data and historical data of the study area. On one hand, the ranking of index importance and weight calculation are more consistent with the actual situation in the study area through the order relation analysis method; on the other hand, the weight of evidence is used to obtain the original data to compensate for the impact of subjective evaluation, making the final evaluation results more accurate and reasonable. Finally, the two weights are combined using the Euclidean distance function to improve the accuracy of the model.

Order Relation Analysis Method (G1 Method)
Order relation analysis (G1) is a subjective evaluation model that is mainly used to rank the importance of two adjacent indicators by means of expert scoring, so as to further determine the weight of each index. When too many evaluation indicators are selected, the G1 method can, to a certain extent, avoid the disadvantage that the weight is difficult to be determined due to too many evaluation indicators and a large amount of data, and effectively avoid the problem of excessive differences between indicators. The calculation process is clear, which greatly simplifies the calculation steps [37]. The algorithm steps of the G1 method are as follows: (1) Deterministic order relation Assuming that the evaluation object is A, the corresponding evaluation indicators are named X1, X2, X3 … Xm. If the importance of Xi is not inferior to Xj based on a certain evaluation criterion or evaluation object, it is expressed as Xi ≥ Xj (the symbol ≥ means not inferior to).
Based on a certain criterion or evaluation object, if the following relationships exist between indicators, it is said that a relationship of importance between evaluation indicators X1, X2, X3 … Xm is determined. X1 * ≥ X2 * ≥ … ≥ Xm * The specific operation steps are as follows: Assume there are m factors. The decision-maker first selects the most important factor from the m factors and marks it as X1 *.
The decision-maker continues to select the most important factor from the remaining m − 1 factors and mark it as X2 *.
Repeat the above action until all the factor markers are completed.
rk is a mood operator, and the nine-level operator method is commonly used in China. The importance assignment of specific indicators is illustrated in Table 3. The expert gives the assignment of rk, and uses the following formula to calculate the weight value W of the index. The calculation formula is as follows: Suppose there are m indicators in total, and Wm represents the weight value of the m-th indicator. The subjective weight of other m−1 indicators is calculated by the following formula:

Weight of Evidence (WOE)
Weight of evidence is an objective weight calculation method, which is mainly driven by the data of Bayesian conditional probability, and conducts superposition analysis by collecting various existing information. The model covers image analysis, mathematical statistics, artificial intelligence, etc., and is organically combined to effectively predict the study area based on the GIS platform. Weight of evidence was initially used in the prediction of mineral target area [38], and then gradually applied by some scholars to the evaluation of collapse and landslide geological hazards.
(1) Calculate a prior probability The prior probability is expressed as the probability of collapse in the study area, which is calculated by the distribution of existing hazard points. The author uses ArcGIS to divide the study area into cells; the cell division standard is to divide the study area into several grid cells with the same shape and equal areas. In order to maximize the accuracy of the data, according to the resolution of the original data, this paper divides the whole study area into cells with an area of 30 m × 30 m. The regional statistical tool in the regional analysis tool is used to conduct joint statistics on the existing hazard point data and the divided cell data to determine the number of units containing geological hazard points.
Assuming that the study area has a total of A cells, the number of cells containing geological hazard points is determined to be N through the zoning statistical function, then the formula for calculating the prior probability is as follows: A is the total number of cells in the study area, N is the number of cells containing collapse geological hazard points, and P(D) is a prior probability.
(2) Weight calculation We assume that the whole study area is divided into A grid cells with the same size and area within the study area. These grid cells are classified and counted according to the classification of influence factors. Among them, the number of grids with geological hazards is recorded as D, the number of grids without geological hazards is D -, and B is the number of grids with geological hazards under a certain impact factor classification. The number of grids without geological hazards under this impact factor classification is B -. For any index weight classification, the weight is calculated as: where represents the weight value of the area where the impact factor exists, and represents the weight value of the area where the impact factor does not exist. The larger the value of , the closer the relationship with the occurrence of hazards. The smaller the value of , the less close the relationship with the occurrence of hazards. P represents conditional probability, that is, the probability of occurrence and non-occurrence of geological hazards under different conditions.
In the formula, C is the comprehensive weight value, and the larger the C value, the higher the contribution of the evidence factor to the occurrence of hazard. The lower the C value, the lower the contribution of the evidence factor to the occurrence of hazard.
(3) Calculate a posteriori probability For m evidence factors, the probability of geological hazards occurring in any unit K in the study area is: where = ln⌈ /(1 − )⌉, is a prior probability, and is the weight value of the ith factor.
Calculation of posterior probability: The magnitude of the posterior probability value reflects the hazard of geological hazards. The greater the posterior probability is, the higher the hazard is. The smaller the posterior probability is, the lower the hazard is. See Table 4 for the distribution of evidence weight factors and their weight values. Because the final result of the weight of evidence is the posterior probability of the secondary evaluation factor, the weight value of the primary evaluation factor is consid-ered. Therefore, based on the weight of evidence, this paper uses the coefficient of variation method to calculate the weight value of the first-level evaluation factor of collapse hazard.

Coefficient of Variance (CV)
The coefficient of variation, also known as the standard deviation rate, is to calculate the standard deviation and average value of a group of information by using the information contained in each factor, and use the ratio of the two as an objective weighting method for the weight value of the influence factor. The calculated weight value can effectively reflect the gap of evaluation indicators [39].
The specific calculation steps are as follows: (1) Standardized treatment.
Because the dimensions of each factor are different, in order to avoid the impact on the results, the first step is to standardize the first-level evaluation indicators as dimensionless.
Negative indicators: Aij is the i-th, standardized value of the j-type evaluation factor, Aij is the i-th, original value of the j-type evaluation factor, ( ) is the minimum value of the i-th evaluation factor, and ( ) is the maximum value of the i-th evaluation factor.
(2) Calculate the coefficient of variation and weight value. Average: Standard deviation: Coefficient of variation: Weight value: On the basis of the posterior probability of the existing second-level evaluation indicators, combined with the coefficient of variation method, the weight value of each first-level evaluation index is further calculated, as shown in Table 5.

Combination Weight
For the evaluation of collapse geological hazards in the study area, first, according to the actual situation in the study area, fully consider the characteristics of the evaluation index themselves, and calculate the subjective weight. Secondly, consider the existing data information, objectively evaluate the existing data within the study area, and reflect the data characteristics of each evaluation index through objective weight.
In order to ensure that the evaluation results are both subjective and objective, the two weight calculation methods are combined. On the basis of the Lagrangian function and European distance function, a combined weight model is constructed to calculate the combined weight value.
(1) Build a combination weight model.
where is the combined weight value, and are the subjective weight value and objective weight value, respectively, and a and b are the optimization coefficients of subjective weight and objective weight, respectively.
(2) Based on the Euclidean distance function, establish the equations.
(3) Calculate the combined weight value.
According to the weight value calculated by the subjective weight G1 and objective weight WOE-CV, the combined weight value of each evaluation index is calculated by simultaneous equations. See Table 6 for details.

Evaluation Index Correlation Analysis
In the selection of an evaluation index for collapse geological hazards, if the correlation coefficient of the two groups of evaluation indices is too high, it will have a great impact on the final result. Therefore, in order to ensure the independence and objectivity of each evaluation index, improve the accuracy of the prediction model, and reduce the interaction between the factors, the Pearson correlation coefficient method is used to analyze the correlation of the evaluation index.
Correlation coefficient, also known as product correlation coefficient, is a linear correlation coefficient that can be used to reflect the degree of correlation between two continuous variables with normal distribution. The size of the correlation coefficient reflects the strength of the correlation between indexes. When the correlation coefficient between indices is less than 0.6, it is considered that the correlation is weak and has li le impact on the accuracy of the evaluation model [40].
In this paper, the ArcGIS software's spatial Analyst Toolset multivariate analysis band set statistical tool is used to conduct statistical analysis on the grid data of the influencing factors, and calculate the covariance and correlation matrix. The origin drawing software is used to draw a hot spot map of the index correlation, and the results are shown in Figure 5. According to the analysis and statistics of various data in the figure, the correlation coefficient between volcanic earthquake and elevation and distance is slightly higher, about 0.57. The correlation coefficients among other evaluation indexes are less than 0.6, which indicates that the correlation and influences of the selected evaluation indexes are not significant.

Comparative Analysis of Evaluation Models
The results show that the hazard point density of the three models increases with the increase in hazard grade. At the same time, through the horizontal comparison of hazard point density, it can be found that in the high-and very high-hazard areas, the hazard point density was significantly improved after the combination of the G1 method and the WOE-CV method through the Euclidean distance function. It shows that the accuracy of the combined prediction model is higher, and the prediction result is more reasonable, as shown in Table 7. Combined with Figures 1, 3 and 5, the numbers of hazard points in high-and very-high hazard areas predicted based on the combined weight model are relatively dense, mainly distributed in the northeast and west of the study area and along the national highway of Changbai County, showing a zonal distribution, which is consistent with the actual survey results. The distribution of faults and roads in this area is relatively dense, and the geological environment is poor, which is the main reason for the frequent occurrence of environmental geological hazards in the study area.

Hazard Map Analysis of Collapse Geological Hazard
According to the calculation results of the combined weights, the grid data of each evaluation index are superposed and summarized using the grid calculator tool through ArcGIS, and finally a hazard assessment zoning map of collapse geological hazards in the study area is obtained. In order to more intuitively reflect the distribution of hazard areas in the study area, the hazard level in the study area is divided into four levels-low, medium, high and very high-by the natural breakpoint classification method. See Figure 6 for details. The statistical results of each zone are used to conduct statistical analysis on the classified area, area proportion, number of hazard points and proportion of hazard points in each zone, and the following table is obtained. See Table 8 for specific results. The results show that the high-hazard area and very-high hazard area account for 36.48% of the total area of the study area, and the number of hazard points in the high-hazard area and very high-hazard area accounted for 87.72% of the total number of hazard points in the study area; 47.59% of the collapse hazard points are located in high-hazard areas, and 40.13% of the hazard points are located in very-high hazard areas, indicating that the accuracy of the model is high. From the distribution of hazard levels, the medium-and low-hazard areas are mainly distributed in the north of the Changbai Mountain Scenic Area and the northwest of Changbai County, while the high-hazard area and very high-hazard area are mainly distributed near the national and provincial trunk roads and densely populated areas in the study area. This further shows that the collapse geological hazards in the study area are mainly affected by human engineering activities.

Model Accuracy Verification
The receiver operating characteristic curve is an index reflecting the sensitivity and specificity of continuous variables based on a series of different binary data, which can effectively negate the impact of human factors on the accuracy of results, and has good objectivity [41]. In the last several years, it has been widely used to evaluate the accuracy of geological hazard prediction models. The area of the ROC curve is expressed by the AUC value, which ranges from 0 to 1. When the value is closer to 1, the accuracy of the model is higher.
In  The AUC values of the G1 and WOE-CV methods are 0.739 and 0.753, respectively, and the AUC value of the combination is 0.817, which shows that the accuracy of the model is significantly improved by the combination of subjective and objective objects through the European distance function.

Conclusions
Based on an actual filed survey of the study area, at the same time, taking into account the various anomalies that have occurred around the Changbai Mountain Tianchi volcano in the study area in recent years, this article selected 12 indicators, namely, slope, aspect, elevation, curvature, lithology, land use type, NDVI, annual average precipitation, volcanic earthquake, distance from road, river and fault, to evaluate collapse geological hazards in the study area, and drew a hazard zoning map. The transcendental probability of Changbai Mountain volcanic earthquakes greater than VI was used to reflect the relationship between Changbai Mountain volcanic earthquakes and collapse geological hazards in the study area. Then, the G1 and WOE-CV models were combined, and the models before and after the combination were compared to draw the following inferences.
The AUC value of the combined model reached 0.817, which is significantly higher than that of the previous G1 method (0.739) and the WOE-CV method (0.753), indicating that the combined model is more accurate and can objectively and effectively reflect the actual situation of the study area. At the same time, it also shows that the evaluation method based on the combined subjective and objective weights is suitable for the hazard analysis of collapse geological hazards in the study area.
According to the collapse geological hazard map established according to the combined weight model, the low-, medium-, high-and very high-hazard areas account for 19.33%, 44.19%, 33.85% and 2.63% of the total area of the study area, respectively. At the same time, high-hazard areas and very high-hazard areas account for 87.72% of the total number of known collapse geological hazard points in the study area; 47.59% of the collapse hazard points are located in the high-hazard area, and 40.13% of the hazard points are in the very high-hazard area.
Based on the in-field photos and hazard zoning maps of the hazard points, it can be inferred that most of the hazard points are mainly distributed near the highway, with some being important transportation arteries. Once a geological hazard occurs, it will cause road damage, affect traffic, and endanger the safety of pedestrians passing by. There are also some collapse geological hazard points distributed around villages and farmland. Although there are some warning signs and protective measures are enforced on this site, investigations have found that some protective works, such as gabion retaining walls and defense nets, are in disrepair, with weakened protective functions and significant safety hazards.
By comparing the hazard zoning map and the distribution of known geological hazards of collapse in the study area, it is found that low-and medium-hazard areas are far away from densely populated areas, such as cities and towns, while high-hazard areas and very high-hazard areas are mainly distributed in densely populated areas. Once a collapse geological hazard occurs, it will cause huge economic and property losses. Therefore, it is recommended to take corresponding protective measures for key areas, strengthen hazard monitoring and a ach importance to the deployment of hazard prevention and mitigation work.
In the process of calculating the volcanic earthquake evaluation index in this study, the earthquake intensity a enuation formula in eastern China is used, which is employed across a relatively wide range, and there may be certain deviations in the calculation accuracy. At the same time, in the process of acquiring remote sensing images, there are problems such as the large time span and insufficient clarity of the image data. Therefore, in future research work, the applicable range of the seismic intensity a enuation formula can be further narrowed to support more comprehensive data analyses and more accurate results. At the same time, advanced remote sensing monitoring methods, such as synthetic aperture radar (SAR), unmanned aerial vehicle remote sensing technology, system integration technology, etc., can be used to eliminate data biases or errors, improve the efficiency of data analysis, and thus improve the accuracy of geological hazard assessment results.
Author Contributions: Writing-original draft preparation, methodology, formal analysis, and validation: Zengkai Pan and Qiuling Lang; writing-review and editing: Yichen Zhang, Jiquan Zhang, Chenglong Yu and Chenyang Wu. All authors have read and agreed to the published version of the manuscript. Acknowledgments: The authors would like to thank the editor and anonymous reviewers for their comments and suggestions, which helped to significantly improve this paper.

Conflicts of Interest:
The authors declare no conflict of interest.