Next Article in Journal
Improved Hydrological Loading Models in South America: Analysis of GPS Displacements Using M-SSA
Previous Article in Journal
Assessing Marginal Shallow-Water Bathymetric Information Content of Lidar Sounding Attribute Data and Derived Seafloor Geomorphometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weighted Information Models for the Quantitative Prediction and Evaluation of the Geothermal Anomaly Area in the Plateau: A Case Study of the Sichuan–Tibet Railway

1
Key Laboratory of Digital Earth Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2
University of Chinese Academy of Sciences, No.19(A) Yuquan Road, Shijingshan District, Beijing 100049, China
3
School of Earth Science and Resources, China University of Geosciences (Beijing), Beijing 100083, China
4
Beijing Key Laboratory of Land and Resources Information Research and Development, Beijing 100083, China
5
China Railway Eryuan Engineering Group Co. Ltd., Chengdu 610031, China
6
China Railway First Survey and Design Institute Group Co. Ltd., Xi’an 710043, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(9), 1606; https://doi.org/10.3390/rs13091606
Submission received: 6 March 2021 / Revised: 30 March 2021 / Accepted: 14 April 2021 / Published: 21 April 2021

Abstract

:
The prediction of geothermal high-temperature anomalies along the plateau railway will be helpful in the construction of the project and its later management. Taking the Sichuan–Tibet railway as the study area and based on Landsat8 thermal infrared images, map data, and measured data regarding the cause and distribution of geothermal high-temperature anomalies, through correlation analysis, we selected six impact factors including the LST, combined entropy of geological formation, fault density, buffer distance to rivers, magnetic anomaly, and earthquake peak acceleration as the input maps of the model. The index-overlay information model, the weights of the entropy information model, and the weights of the evidence information model were established to quantitatively predict the geothermal anomaly in the study area, and the prediction maps were divided into four classes. The results show that the weights of the evidence information model achieved a high prediction accuracy; the success index and the ratio of the high anomaly area reached 0.0053% and 0.872, respectively, and the spatial distribution of the geothermal points is basically consistent with the prediction results. This research can act as a reference for the design and construction of the Sichuan–Tibet railway.

1. Introduction

High-temperature geothermal anomalies will seriously affect the health and safety of operators and the efficiency of mechanical equipment in the process of railway construction, reduce the durability of building materials, and threaten the stability and safety of deep excavation tunnels [1]. The Sichuan–Tibet railway is a major strategic project as part of China’s implementation of a comprehensive transportation system for the development of the western region. The whole line starts from Chengdu in the east and passes Ya’an, Kangding, Changdu, Linzhi, Shannan, and Lhasa in the west, with a total length of about 1543 km. The railway must pass through an area with the most complex geological and geomorphic environments in the world. It has the characteristics of significant terrain height differences, strong plate activity, frequent mountain disasters, and a sensitive ecological environment [2,3,4,5,6,7]. Moreover, due to the large length and buried depth of the railway tunnel, its construction and maintenance process is greatly affected by high-temperature disasters [8]. Therefore, it is very important to predict and evaluate the geothermal high-temperature anomaly area using modern space exploration technology and a spatial analysis method.
In the study of geothermal anomaly prediction, common methods include the thermal infrared remote sensing inversion method, the mathematical statistical model, the spatial information superposition analysis method, and the geophysical method [9,10,11]. These methods have some problems, such as their single-factor consideration, the complex field and indoor work required, their large proportion of human subjective factors, and their low recognition accuracy [12]. Therefore, it is very important to build a multifactor prediction model to improve the automation and delineation accuracy in identifying geothermal high-temperature anomaly areas. Data-driven or knowledge-driven models based on a Geographic Information System (GIS) and Remote Sensing (RS), such as multiple regression analysis, logistic regression analysis, the information method, and artificial neural networks, have been widely used in metallogenic prediction and evaluation [13], groundwater resource evaluation and prediction [14], landslide risk assessment [15], pollution risk assessment [16], and other fields. For example, Proll-Ledesma used faults, hot springs, resistivity, and other factors to build a Boolean, index-overlay, and fuzzy combination model to evaluate the Los Azufres geothermal area in Mexico; the results show that the highly geothermal area prediction is consistent with the existing geothermal point location [17]. Coolbaugh combined multiple factors to build a logistic regression model to predict and evaluate geothermal areas in Nevada and the western basin; 170 potential geothermal zones were found in Nevada [18]. Using seismic, active fault, Bouguer gravity anomaly, and magnetic anomaly data; a tectonic index-overlay model; and a weight of evidence model, Tüfekçi et al. evaluated the geothermal area of the Anatolia continent where Turkey is located [19]. Noorollahi et al. developed a model for geothermal resource exploration by using fault, hot spring, geothermal gradient, and geothermal flow data to predict the geothermal resources in Akita city and Iwate county in northern Japan [20]. Yousefi selected volcanic rocks, faults, hot springs, hydrothermal alteration zones, magnetic anomalies, and other factors for analysis and calculation [21]. Moghaddam used volcano, fault strike, hot spring, hydrothermal alteration zone, geothermal gradient, and geothermal flow data to build a weight of evidence and credibility function model and used Boolean index-overlay, multilevel index-overlay, and a fuzzy logic model to re-evaluate his findings; he concluded that the fuzzy logic model had the highest accuracy [22,23]. Therefore, it is feasible to quantitatively predict the geothermal anomaly area by evaluating the influence factors closely related to the geothermal area. Among all kinds of geological status and disaster assessment methods, the information model is a commonly used bivariate statistical method which can effectively deal with many factors and difficult-to-quantify natural conditions [24]. However, it is rarely used in the prediction and evaluation of geothermal high-temperature anomaly areas, and the research on the weighted information model is not sufficient.
The purpose of this study is to predict and evaluate the geothermal high-temperature anomaly areas along the Sichuan–Tibet railway. The distribution of geothermal sites is closely related to the Earth’s heat flow [25,26], faults, earthquakes [27,28], water, Bouguer gravity, and magnetic anomalies [29]. Therefore, taking the land surface temperature (LST), combined entropy of geological formation, fault density, buffer distance to river, magnetic force, and earthquake peak acceleration as the input impact factors, based on the information model combined with information entropy and weight of evidence theory, an index-overlay, weights of entropy, and evidence information model is established to predict the distribution of geothermal high-temperature anomalies in the study area. This paper will provide suggestions and decision support for the design and construction of the Sichuan–Tibet railway. The technology and methodology flow chart of this paper is shown in Figure 1.

2. Study Area

The study area covers the whole Sichuan–Tibet railway, which is located between 93.5° E–103.5° E and 28.5° N–32.5° N in Southwest China. It mainly includes the eastern part of the Tibet Autonomous Region and the western part of Sichuan Province, with a total area of 399,760 km2 (Figure 2a). The area is dominated by the plateau mountain monsoon climate due to the high mountains acting as a barrier to water vapor. The mean annual air temperature is 2.4–12.6 °C and the annual total precipitation is 417–935 mm, with 80% of it occurring between May and September [30].
This area has steep terrain, including the Langla, Boshula, Nyainqntanglha, Sejila, Himalayas, and other high mountains, with an average altitude of more than 4000 m. The surface water is intensively distributed, and the Yarlung-Zangbo River, Nujiang River, Lancang River, and Jinsha River run through the area. High and deep canyons are highly developed, the relative height difference of the canyons is mostly 2000–3000 m, and the maximum is more than 5000 m. The basic topographic conditions are shown in Figure 2b.
The geological conditions in the area are complex, with many types of formation lithology; sedimentary rocks, metamorphic rocks, and intrusive rocks are widely distributed. The Quaternary (Q) strata covered by the surface are mainly sandy soil and gravel soil; the underlying bedrock developed from the Proterozoic to the Cenozoic period. The main active faults include the Lancang River, Nujiang River, Bianba-Luolong, Jiali, and Miling fault zones, which are deep, large, and concentrated. The earthquake activity is frequent and characteristic of a high-intensity earthquake-prone region. Additionally, the study area is located in the Mediterranean-South Asia geothermal anomaly zone, one of the most intense areas of geothermal activity in the Chinese mainland, and the temperature of the outlying hot springs is about 40 °C [31]. To sum up, the geological and geomorphic conditions of the surrounding area of the Sichuan–Tibet railway are intricate and the impact of geothermal high-temperature disasters on the railway design and construction is obvious. Thus, in this area it will be highly useful to carry out geothermal high-temperature anomaly prediction research, so this space range is listed as the study area.

3. Methodology

3.1. Factor Selection

To examine the genesis and spatial distribution of the geothermal anomalies in this area, our original data included Landsat8 images, a lithology map, a fault map, a river map, a magnetic map, Bouguer gravity data, and earthquake acceleration data. In order to highlight the spatial relationship between these data and geothermal anomalies, they were transformed into a land surface temperature map, a combined entropy of geological formation map, a buffer distance to fault map, a fault density map, a buffer distance to river map, a magnetic map, a Bouguer gravity map, and an earthquake peak acceleration map.

3.1.1. Geothermal Points

The hot and warm springs exposed on the land surface were effective for the characterization of geothermal high-temperature anomalies. These springs and their surrounding areas can be regarded as known geothermal anomaly areas and were thus used in the calculation of subsequent prediction models. The geothermal points in the study area were obtained from field measurements, mainly comprising measurements of exposed hot spring points and drilling points (Figure 3), with a maximum drilling depth of 1575 m. Hot springs measure the temperature of water at a certain depth with portable noncontact thermometers, and drilling points measure the temperature of bedrock at different drilling depths. When the average temperature of a point reaches 30 °C, it is considered as a geothermal anomaly point. There are 249 points in total, with the highest temperature reaching 94 °C and an average temperature of 41 °C (Figure 2).

3.1.2. Land Surface Temperature

The land surface temperature (LST) is mainly influenced by solar radiation and geological activities, which can extract heat information from the Earth’s surface, thus providing a basis for the determination of geothermal anomalies [32]. LST information generated from remote sensing images has been successfully applied to the study of geothermal resources and disasters. Therefore, it is feasible to use an LST map as a representative factor influencing terrestrial heat flow. Thermal infrared remote sensing is an important means to obtain LST. Satellite systems that can realize temperature inversion by using thermal infrared bands include MODIS, Landsat8, Landsat7 ETM+, ASTER, NOAA/AVHRR, etc. Among them, Landsat8 data have many advantages: firstly, these data can be obtained directly (http://glovis.usgs.gov/, accessed on 15 August 2020); secondly, thermal infrared band 10 has a high spatial resolution of 30 m and can effectively obtain the fine thermal field landscape and surface temperature anomalies. LST inversion methods using remote sensing images can be divided into single-channel algorithms [33], multichannel algorithms [34], multiangle algorithms [35], multiphase algorithms [36], and hyperspectral inversion algorithms [37].
Based on the radiative transfer equation, Qin et al. established a single-window algorithm for LST inversion using only one thermal infrared band [38]. The algorithm only needs three parameters—namely, surface emissivity, atmospheric transmittance, and average atmospheric temperature—to retrieve LST. For the surface emissivity, the normalized difference vegetation index (NDVI) threshold method was used to calculate it. In the scale of the Landsat8 satellite, the Earth’s surface can be roughly regarded as being covered by four types of ground features: water, vegetation, buildings, and soil. According to the different values of NDVI, we can determine the proportions of various ground features in the mixed pixels to obtain the surface emissivity. The atmospheric transmittance, according to previous studies, mainly depends on the atmospheric water vapor content, and w can be obtained through the MOD05_L2 product. For mid-latitude winter, the relationship between atmospheric transmittance and atmospheric water vapor content is:
Τ = 0.9228 − 0.0735w,
where T is the atmospheric transmittance and w is the atmospheric water vapor content.
For the average atmospheric temperature, Qin et al. proposed a linear relation for the approximation of the average atmospheric temperature from the near surface air temperature. For mid-latitude winter, the average atmospheric temperature could be approximated as:
Ta = 19.2704 + 0.9112T_0,
where Ta is the average atmospheric temperature and T_0 is the near-surface air temperature.
The single window method is widely used and has a high inversion accuracy, and most of the satellite sensor data with a thermal infrared band can be used for LST inversion by this algorithm [38]. The formula can be expressed as:
Ts = {a(1 − Ci − Di) + [b(1 − Ci − Di) + Ci + Di] Tb − DiTa}/Ci,
C i = ε i τ i ,
D i = ( 1 τ i ) [ 1 + τ i ( 1 ε i ) ] ,
where Ts is the retrieved LST; Tb is the brightness temperature (k) obtained by the sensor; Ta is the average atmospheric temperature (k); a and b are the linear regression coefficients, which are related to the temperature range of the study area; C and D are intermediate variables; ε is the surface emissivity; and τ is the atmospheric transmittance [39].
The inversion results of LST in plateau and mountainous areas are significantly affected by solar radiation heating, which covers the effect of geological activities on LST changes and needs to be eliminated by topographic correction. Empirical statistics is a commonly used terrain correction method. There is a certain linear relationship between the solar radiation and the radiance value of image pixels. Through regression analysis, the linear relationship between them was established, and then the radiation energy received by the slope pixels was corrected to the horizontal position through the regression relationship so as to eliminate the effect of terrain. The calculation formulae are as follows:
cos ( i ) = cos ( z ) cos ( S ) + sin ( z ) sin ( S ) cos ( Φ x - Φ n ) ,
LT = mcos(i) + b,
L H = L T - [ mcos ( i ) + b ] + L T ¯ ,
where i is the effective incidence angle of the sun; z is the solar zenith angle; Φ x is the sun azimuth angle; S is the pixel slope angle; Φ n is the pixel aspect angle; L T is the radiance value of the ground object before correction; m and B are the parameters obtained by regression analysis; L H is the radiance value after correction; L T ¯ is the radiance value of ground features in a flat area.
The study area needs 15–20 images for full coverage; a total of 104 winter Landsat8 images were obtained in the periods November 2013–February 2014, December 2014–January 2015, December 2015–February 2016, November 2016–January 2017, December 2017–February 2018, and November 2018–January 2019. LST values were obtained by using the single-window algorithm, and terrain correction was carried out. The multiyear average LST (Figure 4a) was calculated as the input data of the geothermal anomaly prediction model. Since the study only focuses on high-temperature anomalies, the lower threshold of LST was set to −10 °C. It can be seen that the low-temperature anomalies were mainly concentrated in the mountains with high altitudes, while the high temperature anomalies were distributed in the valleys and plains.

3.1.3. Combined Entropy of Geological Formation

The comprehensive variable constructed by relative entropy can reflect the lithologic combination characteristics in the unit area, which can quantitatively characterize the degree of homogeneity or variation in a geological structure in a certain range [40,41]. Entropy anomaly can be used as an important index to reveal geological anomaly. The combined entropy of geological formation is the entropy anomaly of various geological bodies or different attributes of the same geological body in terms of unit area or volume, which can be calculated using the lithology data (http://geocloud.cgs.gov.cn/, accessed on 2 October 2020). In a certain range, the higher the entropy value is, the higher the variation degree of a geological structure is, and the higher the possibility of geothermal anomaly is. The lithology data in this study can be divided into 323 categories. The calculation of combined entropy can be divided into the following steps (Figure 5).
First, the lithologic map needs to be divided into several grid units. The long axis direction, shape, and size are the elements of unit division. The direction of the long axis of the grid unit should be consistent with the direction of the regional tectonic line. The main fault direction of the study area is southeast, so the long axis of the unit is also southeast. The shape of the grid unit should correspond to the shape of the stratum. The stratum in the area has an obvious strike, and the ratio of width to length is about 1/2–1/3, so the width/length ratio of the grid unit is 1/2. After determining the grid units, each independent lithology area in the unit is calculated, then the sum of their area in the unit is calculated and the ratio x i (i = 1,2,3…, n) of each lithology area in the unit to the unit area is calculated. On this basis, the combined entropy of geological formation can be evaluated by the following formula:
E jk = i = 1 n x i lnx i / lnn ,
where n is the lithology type existing in the grid unit and j and K are the row number and column number of the unit. The calculated combined entropy of geological formation is shown in Figure 4b. The high values for the combination entropy are concentrated in the linear area with sharp changes in lithology, which has a high spatial correlation with the pattern of geological faults.

3.1.4. Buffer Distance to Fault and Fault Density

It is well known that deeper faults form seepage channels where groundwater can flow towards the deep crust and reach high temperatures. When the heated water returns to the surface again through shallow faults, geothermal occurrences are possible [42]. Based on the fault data (Figure 2b, http://geocloud.cgs.gov.cn/, accessed on 5 October 2020), the buffer distance to fault and fault density can be used to determine the relationship between geothermal points and active thermal geology. In the buffer distance map, each pixel contains the vertical distance information from the nearest fault. Figure 4c clearly shows that there are 208 geothermal points in total, accounting for 83.53% of the total number of points, which are distributed within 4 km from the fault. It can be seen from Figure 4d that there are 154 geothermal points, accounting for 61.85% of all points, which are distributed in the area with a fault density greater than 0.05 km/km2.

3.1.5. Buffer Distance to River

The different migration modes of water channels have different effects on geothermal formation. The downward migration of groundwater leads to a decrease in rock temperature deep underground, where it is not easy to generate geothermal high-temperature anomalies; when the groundwater heated by the surrounding rock circulates upward under certain geological conditions, it will cause the local temperature of surrounding rock to rise and form a high-temperature anomaly in the shallow part of the surface [43]. Therefore, it is feasible to use the distance to river system as the influencing factor for geothermal areas. The data of rivers in the study area with a scale of 1:3 million were divided into 3 classes, including 3-level rivers (basin area of 1–10 thousand km2), 4-level rivers (basin area of 0.1–1 thousand km2), and 5-level rivers (basin area <0.1 thousand km2) [44]. In order to study the spatial relationship between the river system and geothermal activity, the river was converted into buffer distance (Figure 4e). Statistics show that 85.54% of geothermal points (213) are distributed within a range of less than 4 km from the river, indicating a close relationship between river and geothermal anomalies. The abundant water resources in the study area provide a material and energy basis for the occurrence of geothermal high-temperature anomalies.

3.1.6. Magnetic Anomaly

Magnetic anomalies are often used to indicate the active area of underground hot water and the areas with significant changes in tectonic stress, so the magnetic layer can be used as the influence factor [45]. In theory, lower aeromagnetic values are closely related to geothermal area; in the range of geothermal water activities, the magnetic properties of rocks are reduced due to thermal alteration; on the other hand, under the action of tectonic stress, the magnetic properties of rocks along the stress direction will be weakened, so the magnetic properties of rocks in the tectonic fracture zone will also be reduced. According to the magnetic anomaly map (Figure 4f) [46], most of the geothermal points are distributed in the medium magnetic region (−70–−80 nT), and the points with high magnetism (>100 nT) and low magnetism (<−50 nT) are relatively few.

3.1.7. Bouguer Gravity

Bouguer gravity anomaly is the area where the density of crustal material changes sharply along the horizontal direction, which is a sign of the existence of a graben system [47]. The reason for the formation of this is the uneven distribution of underground rock mass and mineral density, or a density difference between the geological body and surrounding rock. Bouguer gravity anomaly can be used to understand regional structures and delineate large fault structures and possible hydrothermal activities. The Bouguer gravity anomaly can reflect the variation of crustal thickness in different areas, and the rock density decreases with the increase in temperature. The distribution gradient of Bouguer gravity anomaly is consistent with the distribution of hot springs, so it is regarded as one of the impact factors. A Bouguer gravity map is shown in Figure 4g.

3.1.8. Earthquake Peak Acceleration

In the seismicity area and active fault zone, due to the high permeability and circulation of water the hydrothermal activity is closely related to the seismic area. The seismicity can be studied by means of the intensity of the earthquake, and then the spatial relationship between the hydrothermal activities in the region can be found [48,49]. The index of earthquake peak acceleration is used to represent the intensity of seismic activity. From the map of earthquake peak acceleration (Figure 4h, http://geocloud.cgs.gov.cn/, accessed on 12 October 2020), we can see that there are more geothermal points distributed in the areas with medium and high earthquake peak acceleration.

3.2. Factor Reclassification and Independence

3.2.1. Factor Reclassification

The reclassification of impact factors is an important part of the subsequent model prediction, and the appropriate threshold is needed to obtain the reasonable impact factor classification map. Generally speaking, the more classes the factor map divides, the more objective and reliable the analysis results are. However, too many classes will increase the computational burden. Therefore, in the process of reclassification, using the natural breakpoint method, each factor map determines seven classes, and the area between the intervals is neither too large nor too small. In this way, all the factors obtain the threshold for the reclassification of impact factors (Table 1). Additionally, when using a square grid as the evaluation unit, the size of the grid should be considered. In theory, the smaller the grid, the more similar the geological environment of each point in the grid, and the higher the accuracy of evaluation and prediction. However, if the size is too small, the calculation efficiency will be seriously reduced. Therefore, when selecting the grid size, the identity of geological environment in the unit, the suitability of the size proportion of the study area, and the computer processing ability are taken into account. According to the experiment, Tang et al. put forward the corresponding grid size selection standard (Table 2) [50]. Combined with the scale of the study area, this study selects 100 × 100 m as the grid size and resamples the factor map after reclassification (Figure 6).

3.2.2. Independence Test

The establishment of an information model needs a strict assumption—that is, that the selected influence factors are independent of each other [51]. In this paper, based on the results of the factor reclassification and factor analysis, the independence was judged—that is, the correlation coefficient was used to measure the correlation of factors, the correlation between each factor was analyzed by calculating the covariance and correlation matrix, and the factors with a large correlation were excluded to ensure the effectiveness of the model application. When the correlation coefficient of each factor satisfies |R|≤ 0.3, it can be considered as a weak correlation or no correlation [52]. From Table 3, it can be seen that the correlation between he Bouguer gravity and magnetic map is high, so the Bouguer gravity factor is excluded; the buffer distance to the fault and fault density are generated by fault data, so the correlation is the highest. When |R| reaches 0.771, the buffer distance to the fault factor is excluded. Finally, the LST, combined entropy of geological formation, fault density, buffer distance to rivers, and magnetic and earthquake peak acceleration are retained as the influencing factors in the model.

3.3. Model Establishment

Based on the information model combined with a variety of objective weighting methods, the index-overlay information model, weights of entropy information model, and weights of evidence information model are proposed. These models are used to predict and evaluate the geothermal anomaly in the study area.

3.3.1. Information Model

The main idea of the information model applied to the evaluation of geothermal anomaly area is to combine the geological conditions closely related to the formation of high-temperature geothermal anomalies, convert the values of these influencing factors into information values, and calculate the information value provided by each influence factor on the generation of high-temperature geothermal anomalies to analyze the relationship between the two so as to combine the geological geomorphology, remote sensing information, and other elements to build a prediction model for geothermal high-temperature anomalies [24]. The information model is a bivariate statistical analysis method, geothermal anomaly is affected by many factors that will form a kind of “optimal combination of factors” to cause the occurrence of a geothermal anomaly [53]. Its formation is related to the quantity and quality of influencing factors; the greater the amount of information there is, the greater the possibility of forming a geothermal anomaly.
When the information model is used to evaluate the geothermal anomaly area, n factors affecting geothermal energy are selected and expressed as X i (i = 1,2,3…, n). The geothermal information provided by single factor X i is I i and the total geothermal information provided by n factors is I. I i and I are expressed as follows:
I i = log 2 P [ Y / X i ] P ( Y ) ,
I = log 2 P [ Y / X 1 X 2 X n ] P ( Y ) .
Based on 3.99 × 107 grids with equal area of 100 × 100 m in the study area, the amount of information provided by each class of influencing factor X i (i = 1,2,3…, n) to geothermal anomaly is I i , which can be expressed as follows:
I i = log 2 N ij / N S ij / S ,
where S is the total number of grids in the study area, S ij is the number of grid units with X ij of j class, n is the total number of geothermal points in the study area, and N ij is the number of geothermal points of X ij with class J.
The total information value I of the information model is as follows:
  I = i = 1 n I i = i = 1 n log 2 N ij / N S ij / S ,
where n is the number of influence factors.
After calculation, the information values of each factor are shown in Table 4. The + and −information values indicate the positive and negative correlation with geothermal high-temperature anomalies. For example, the information value of the LST 1st class is −3.145, which indicates that this class will greatly lessen the occurrence of geothermal high-temperature anomalies, while the information of the seventh class is 0.926, and the geothermal temperature of this class will promote the occurrence of geothermal anomaly. With the increase in the series of LST, the combined entropy of geological formation, the fault density, and the earthquake peak acceleration, the amount of information provided increases, which indicates that these factors have a positive correlation with geothermal high-temperature anomalies as a whole. The larger their values are, the more likely it is that a geothermal high-temperature anomaly will occur. Among them, the positive value of the information provided by LST, fault density, and earthquake peak acceleration is very prominent, up to 1.203. There is a negative correlation between the buffer distance to rivers and magnetic anomalies and geothermal high-temperature anomalies. With the increase in their values, the possibility of geothermal anomaly decreases.

3.3.2. Index-Overlay Information Model

Due to the difference in the quality and quantity of geological and geophysical data, the contribution of influencing factors to the formation of geothermal anomaly is different. However, the traditional information model simply superimposes the information values of each factor, and it is difficult to reflect the distinction of the factors’ influence on geothermal anomalies. Therefore, several objective weighting methods are used to establish the weighted information model and discuss the different contribution degrees of each factor. The index-overlay is a method that directly overlays the information value of the factors [54], and its expression is as follows:
Index = i = 1 n W i × I i = i = 1 n W i × log 2 N ij / N S ij / S ,
where Index is the superposition value of weighted information and W i is the weight of each influence factor. In the index-overlay model, all the factors have equal effects, W i = 0.167.

3.3.3. Weights of Entropy Information Model

Using the theory of information entropy, objective weight can be obtained to reflect the importance of factors. According to the theory of information entropy, the smaller the entropy value of a factor is, the greater the change in the density of geothermal points in each classification is, and the higher the corresponding weight is, the greater the effect of this factor will be on the prediction and evaluation objectives [55]. The calculation formula for geothermal activity is:
d ij = N ij S ij ,
where d ij is the geothermal occurrence rate of the j level in the ith factor layer and the meanings of N ij and S ij are consistent with Formula 10.
The formula used to calculate the normalization value for points can be expressed as follows:
K ij = d ij j = 1 n d ij ,
where K ij is the normalization value for the jth class of the ith map.
According to the informational entropy theory, the entropy is:
  H i = j = 1 m K ij lnK ij lnm ,
where H i is the theoretical value for the ith map. If K ij = 0, then ln Kij = 0 and H i = 0.
The weight of the ith map is:
W ij = ( 1 H i ) × n n i = 1 n H i .
The weights of entropy information mode can be expressed as follows:
  Entropy = i = 1 n W i × I i = i = 1 n W i × log 2 N ij / N S ij / S ,
where Entropy is the calculated result for the pixels and W i is the weight of the ith map (Table 4).

3.3.4. Weights of Evidence Information Model

The weight of evidence is a quantitative evaluation method based on Bayesian statistics. Through the spatial correlation analysis of some geothermal anomalies that have occurred and the related factors affecting the formation of geothermal anomalies, the contribution rate (weight) of each factor to the geothermal anomaly can be determined and used to calculate the possibility of geothermal anomaly. Finally, each factor is given the corresponding weight, and the geothermal anomaly is obtained by superposition. Its advantage is that the interpretation of comprehensive weights is relatively intuitive and easy to understand. By introducing the weight of evidence into the information model, more objective and accurate factor weights can be obtained, which makes the prediction results more accurate [56].
The probability of the random occurrence of training points in the whole study area is defined as the prior probability and is applicable to the whole research area, and its calculation formula is as follows:
  P ( d ) = M T ,
where P(d) is an a priori probability, M is the number of evaluation units with geothermal points, and T is the total number of evaluation units in the study area.
The weight of the evidence layer can be defined as:
W ij k = { In P ( e ij | d ) P ( e ij | d ¯ ) , k = + In P ( e ij ¯ | d ) P ( e ij ¯ | d ¯ ) , k = 0 , k = 0
where P(x|y) is the conditional probability of X phenomenon when y phenomenon occurs; e ij is the number of grids in J class of the ith evidence factor; d is the number of grids with geothermal anomaly; e ij ¯ is that j class of the ith evidence factor does not occur; d ¯ is no geothermal anomaly; k is the state of the ith evidence layer in the unit, k = + represents positive correlation weight, and k = is negative correlation weight. When the positive weight is greater than 0 or the negative weight is less than 0, it means that the factor is positively correlated with the geothermal anomaly. When the positive weight is less than 0 or the negative weight is greater than 0, this means that the factor has little influence on the geothermal anomaly. When the positive weight or negative weight is 0, this means that the factor is not related to geothermal anomaly.
Contrast (C) represents the correlation between the evidence and training points. The greater the C is, the closer the correlation is. Positive is defined as favorable to the target, and negative is unfavorable to the target. The contrast can be used as the final weight and combined with the information value. The calculation formula of C is as follows:
C ij = W ij + - W ij .
The weights of evidence information mode can be expressed as follows:
Evidence = i = 1 n C i × I i = i = 1 n C i × log 2 N ij / N S ij / S
where evidence is the predicted value of the weight of evidence method and C i is the weight of the ith factor (Table 5).

3.3.5. Classification Maps of Prediction

In order to compare and analyze the models more intuitively, the information series models need to be further classified. If the model value is lower than 0, this indicates that there is nothing to do with the geothermal high-temperature phenomenon. Therefore, the 1st threshold value is set to 0 and values greater than 0 are divided into three classes according to the natural breakpoint method. The specific segmentation threshold is shown in Table 6, and the classification result is shown in Figure 7. Four in the maps shows the high geothermal anomaly area, where high-temperature anomalies occur frequently, and this is the most serious area affected by geothermal high-temperature disaster; three and two, respectively, indicate the middle and low geothermal anomaly areas, and the possibility and severity of geothermal anomalies in the region are reduced in turn; one indicates that there is no abnormal area, the probability of occurrence of geothermal anomaly is very low, and it is not affected by geothermal high temperature in this area.
According to the classification map (Figure 7), the overall prediction results are relatively close, and the geothermal anomaly area is generally in the northwest–southeast direction, which is highly similar to the distribution of faults and water systems. The highly abnormal areas are mainly distributed in the middle of Ganzi and Changdu and the north of Linzhi. The fault structures and rivers are highly developed in these areas, with a high LST and frequent earthquakes. Among them, the geothermal anomaly areas in Ganzi and Linzhi coincide with the Sichuan–Tibet railway line greatly, which has a severe impact on the design and construction of the railway and is the key area to prevent the occurrence of geothermal anomaly disasters. In some of the details, the difference in the prediction results is more significant. The high abnormal area obtained by the weights of entropy information model and the weights of evidence information volume model is larger, and the classification of the prediction results of the weights of evidence model is more significant, and we can clearly see the spatial correlation between the medium and low abnormal areas and high abnormal areas. The index-overlay model can obtain the smallest abnormal area, and the subtle abnormal changes cannot be effectively displayed in the classification map.
Additionally, there is a certain area of false anomaly in the eastern part of the study area, which is due to the impact of the urban heat source of Ya’an and Chengdu, causing the LST of the area rise. In the model with a large weight of LST, there will be a certain degree of interference in the judgment of geothermal anomalies.

4. Model Assessment

4.1. Analysis of Success Index

The occurrence rate of geothermal training points can be used to assess the effectiveness of the information model—that is, to analyze the success index of models. Compared with the prior probability of geothermal occurrence (0.00063%), it can indicate whether the prediction results are accurate [19]. The prior probability is the ratio of the total number of geothermal points to the total number of grids in the study area. According to the classification of the prediction results, a success index analysis table (Table 7) was obtained. It can be seen from the table that the occurrence rate of geothermal points will also increase with the increase in the anomaly degree of the geothermal anomaly area, which shows that the evaluation results of the information series model are reliable. In the prediction of highly abnormal areas, the weights of the entropy information model do not reflect the advantages, and the success index is lower than that of the index-overlay model; the results of the weights of evidence information model in high and medium abnormal areas are better, indicating that the prediction results are more accurate.

4.2. Analysis of Area Ratio

Area ratio analysis can be used to quantify the prediction results of the model. The effectiveness of the model was evaluated by calculating the area covered by the prediction probability function [57]. According to the attribute values of the prediction results of the information model, the results were arranged from large to small and divided into 25 classes by the equal quantile method. The number of geothermal points in each class and the grid area in the classification interval were assessed, the cumulative percentage curve of the grid units and points was calculated, then the prediction function curve was obtained. The expression of area ratio is as follows:
λ = (2A − P)/(2 − 2P),
where P is the prior probability (0.00063%); A is the area surrounded by the prediction probability function and X-axis; λ is the area ratio of the prediction accuracy of the model. The closer the value is to 1, the more accurate the prediction result is. The calculation results are shown in Figure 8, and the area ratio of the models is 0.859, 0.846, and 0.872, respectively, which indicates that the overall accuracy of the model is higher, but the difference between these models is small. Additionally, the accuracy of the weights of evidence information model is slightly higher than that of the other two models.

4.3. Ground Verification

In order to further verify the effectiveness of the proposed model, we used twelve ground measurement points (Figure 9). Verification points were located in Guoqing, Changdu City, Litang, Ganzi Prefecture and Bayi District, Linzhi City about 2100–8700 m away from the railway line and close to the river, with flat terrain, bare soil, and low vegetation; Points 1–6 belonged to the high geothermal anomaly area predicted by the weights of evidence model and points 7–12 were in the low geothermal anomaly area. The temperature of water and bedrock at a certain depth within 10 m2 around the verification point is measured, and the specific situations of some verification points are shown in Table 8. It can be seen that the temperature range of the verification points in the high anomaly area is significantly higher than that in the low anomaly area. The ground test can verify that the prediction results of the weights of evidence information model are reliable and accurate to a certain extent.

5. Discussion and Conclusions

5.1. Discussion

High-temperature geothermal phenomena will seriously affect the construction of a plateau railway. In order to predict the area of geothermal high-temperature anomalies, a prediction method for the geothermal anomaly area was constructed based on the weighted information model for the Sichuan–Tibet railway area. Compared with the current prediction research on the geothermal anomaly area, we fully considered the geological factors, natural environment factors, and inducing factors of geothermal high-temperature anomalies and selected a number of impact factors closely related to the geothermal are, thus reducing the one-sidedness of single-factor or multifactor analysis.
Information models are widely used in geological hazard risk assessment [58], metallogenic prediction [59], hydrological feature extraction [60], and resource potential assessment [24]. At present, the superposition of factors is often used to accumulate the information of each factor to obtain the final evaluation results. However, the contribution of each factor to the prediction target cannot be the same, so the weighted superposition can comprehensively analyze the influence and relationship of multiple factors on the research target. The weight of evidence method has been extended to many fields, such as resource prediction [56], disaster analysis [61], and ecological environment research [62]. This weight determination method was introduced into the information model to enrich the research on the weight of information model, and it can quantitatively analyze the role of multiple factors. It can achieve better prediction results and greatly reduce the impact of human subjective factors on the prediction results. The model can provide geological information for the design and construction of the Sichuan–Tibet railway.
This paper optimizes the existing prediction methods used for geothermal anomaly areas, but there are still some problems that need to be fixed: on the one hand, the accuracy of the Bouguer gravity and magnetic anomaly data is relatively low, which will affect the prediction results to a certain extent; on the other hand, the prediction results of the geothermal anomaly areas are only classified and discussed in the macro scale. In the next stage of our research, we will make use of emerging geological remote sensing technologies such as reconstruction of Light Detection and Ranging (LiDAR), and the high and medium abnormal areas will be analyzed with a smaller spatial scale and in more detail in order to form specific guidance for the design and construction of the Sichuan–Tibet railway.

5.2. Conclusions

(1) In terms of impact factor selection, Landsat8 image data, various map data, and measured data were converted into LST, combined entropy of geological formation, buffer distance to faults, fault density, buffer distance to rivers, magnetic anomaly, Bouguer gravity, and earthquake peak acceleration data and combined with the measured geothermal points as the impact factors to predict geothermal high-temperature anomalies. After transformation, the spatial distribution of the factors had a certain relationship with the geothermal area: the more the geothermal points that were distributed in the areas with a higher LST, the higher the combination entropy was, the higher the fault density was, the closer it was to rivers and faults, the lower the magnetic values were, the higher the Bouguer gravity was, and the more severe the earthquakes were. In the factor correlation analysis, it was found that the correlation |R| of some factors was >0.3. Therefore, the buffer distance to faults and Bouguer gravity were removed, and the remaining six impact factors were retained as the input maps of the subsequent prediction model.
(2) Based on the information model, the objective weighted method was introduced into the prediction of the geothermal anomaly area and the index-overlay information model, the weights of entropy information model, and the weights of evidence information model are established, respectively. From the calculation results of the information values, different classes of reclassification factors were found to have different positive and negative effects on geothermal anomaly. Overall, with the increase in LST, combined entropy of geological formation, fault density, and earthquake peak acceleration, the values of the information provided were also increased. On the contrary, there was a negative correlation between the buffer distance to rivers and magnetic anomalies and geothermal high-temperature anomaly—with the increase in their values, the information values decreased. From the weight calculation results, the LST, fault density, and buffer distance to rivers accounted for a large proportion of the weights of entropy model, and the prediction results were mainly based on the information provided by these factors; the weight of LST, combined entropy, and fault density in the weights of evidence is relatively large. The results of the two models show that the core factors in judging geothermal high-temperature anomalies are LST and fault, which represent the geothermal flow and geological structure.
(3) The overall predictions of the three models were relatively similar. The geothermal anomaly area ran northwest–southeast, and the high anomaly area was mainly distributed in the middle of Ganzi and Changdu and the north of Linzhi. The fault structures and river systems in these areas are highly developed, the LST is high, and earthquakes occur frequently. Among them, the geothermal anomaly areas in Ganzi and Linzhi have a great impact on the design and construction of the railway. In the prediction details, the highly abnormal area obtained by the weights of entropy information model and the weights of evidence information model was larger, and the classification result of the latter is more prominent; the abnormal area obtained by the index-overlay model was the smallest, and the classification situation was relatively fuzzy in the map. From the analysis of the success index, the three models achieved a good prediction effect. The weights of evidence information model was the best in the prediction of the high abnormal area and the index reached 0.0053%, far exceeding the prior probability (0.00063%); the index-overlay model was better in the medium abnormal area and the index was 0.0021%. From the analysis of the area ratio, we can see that the prediction results of these models were similar, and the weights of evidence model has a slight advantage.

Author Contributions

W.Z. and Z.C.: methodology, writing-original draft preparation and edit-ing, and software. T.F.: writing-review and editing. D.W.: validation and formal analysis. L.J.: in-vestigation and data curation. S.D.: resources and visualization. X.Z.: supervision and project administration. Q.D. and J.C.: funding acquisition. D.M. and M.B.: data arrangement and processing. All authors have read and agreed to the published ver-sion of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No. 41876210), National Key Research and Development Program of China (No. 2017YFA0603003 and No. 2017YFC0601502) and Institute International Cooperation Projects (E03405020N).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to offer special thanks to NASA and USGS for providing us with the Landsat8 images. We also thank the China Railway Eryuan Engineering Group Co. Ltd. for the support of field investigation and geothermal point data they provided. A special acknowledgement should be expressed to China-Pakistan Joint Research Center on Earth Sciences that supported the implementation of this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peng, J.B.; Cui, P.; Zhuang, J.Q. Challenges to engineering geology of Sichuan-Tibet railway. Chin. J. Rock Mech. Eng. 2020, 39, 1–13. [Google Scholar] [CrossRef]
  2. Guo, C.; Zhang, Y.; Jiang, L.; Shi, J. Discussion on the environmental and engineering geological problems along the Sichuan-Tibet railway and its adjacent area. Geoscience 2017, 31, 4–16. [Google Scholar] [CrossRef]
  3. Lu, C.; Cai, C. Challenges and countermeasures for construction safety during the Sichuan–Tibet railway project. Engineering 2019, 5, 833–838. [Google Scholar] [CrossRef]
  4. Jiang, L.W.; Li, Y.S.; Yi, S.J.; Huang, C.; Liu, K. Engineering Geological Research in Formations-structural features of the Sichuan-Tibet Railway Plate Collision Junction. In Proceedings of the 2016 Academic Exchange Conference on Challenges and Countermeasures of Sichuan-Tibet Railway Construction; People’s Communications Press: Beijing, China, 2016. [Google Scholar]
  5. Xu, Y.D.; Yao, L.K. Some cognitions and thinking about the specific geo-environmental problems along the Sichuan-Tibet railway. J. Railw. Eng. Soc. 2017, 34, 1–5. [Google Scholar] [CrossRef]
  6. Zhu, Y. The Sichuan-Tibet Railway Construction Challenges and Countermeasures; China Communications Press: Beijing, China, 2017. [Google Scholar]
  7. Zheng, Z.X.; Sun, Q.Q. The tunneling engineering of Sichuan-Tibet railway. Tunn. Constr. 2017, 37, 138–143. [Google Scholar]
  8. Yan, J.; He, C.; Wang, B.; Meng, W.; Wu, F.Y. Inoculation and characters of rock bursts in extra-long and deep tunnels located on Yarlung Zangbo suture. Chin. J. Rock Mech. Rock Eng. 2019, 38, 769–781. [Google Scholar] [CrossRef]
  9. Watson, T.; Kruse, F.A.; Hummer, M.S. Thermal Infrared Exploration in the Carlin Trend. Geophysics 1990, 55, 70–79. [Google Scholar] [CrossRef]
  10. Rodgers, J.L. The Epistemology of Mathematical and Statistical Modeling. Am. Psychol. 2010, 65, 1–12. [Google Scholar] [CrossRef] [Green Version]
  11. Criss, R.E.; Champion, D.E. Magnetic Properties of Granitic Rocks from the Southern Half of the Idaho Batholith: Influences of Hydrothermal Alteration and Implications for Aeromagnetic Interpretation. J. Geophys. Res. Solid Earth 1984, 89, 7061–7076. [Google Scholar] [CrossRef]
  12. Yan, B.Z.; Qiu, S.W.; Xiao, C.L.; Liang, X.J. Potential Geothermal Fields Remote Sensing Identification in Changbai Mountain Basalt Area. J. Jilin Univ. 2017, 47, 1819–1828. [Google Scholar] [CrossRef]
  13. Chen, J.P.; Lu, P.; Wu, W.; Zhao, H.Q. A 3D method for predicting blind orebodies, based on a 3D visualization model and its application. Earth Sci. Front. 2017, 47, 1819–1828. [Google Scholar] [CrossRef]
  14. Hou, J. Spatial assessment for groundwater quality based on GIS and improved fuzzy comprehensive assessment with entropy weights. Chin. J. Popul. Resour. Environ. 2013, 11, 135–141. [Google Scholar] [CrossRef]
  15. Fan, L.F.; Hu, R.L.; Zeng, F.C.; Wang, S.S.; Zhang, X.Y. Application of Weighted Information Value Model to Landslide Susceptibility Assessment-A Case Study of Enshi City, Hubei Province. J. Eng. Geol. 2012, 20, 508–513. [Google Scholar] [CrossRef]
  16. Ding, E.J.; Xie, D.T.; Wei, C.F.; Ni, J.P. MUSLE-Based Agricultural Nonpoint Source Pollution Risk Assessment in the Three Gorges Reservoir Area of Chongqing. J. Southwest Univ. 2010, 32, 96–101. [Google Scholar]
  17. Prol-Ledesma, R.M. Evaluation of the reconnaissance results in geothermal exploration using GIS. Geothermics 2000, 29, 83–103. [Google Scholar] [CrossRef]
  18. Coolbaugh, M.F.; Raines, G.L.; Zehner, R.E. Assessment of exploration bias in data-driven predictive models and the estimation of undiscovered resources. Nat. Resour. Res. 2007, 16, 199–207. [Google Scholar] [CrossRef]
  19. Tüfekçi, N.; Süzen, M.L.; Güleç, N. GIS based geothermal potential assessment: A case study from Western Anatolia, Turkey. Energy 2010, 35, 246–261. [Google Scholar] [CrossRef]
  20. Noorollahi, Y.; Itoi, R.; Fujii, H.; Tanaka, T. GIS model for geothermal resource exploration in Akita and Iwate prefectures, northern Japan. Comput. Geosci. 2007, 33, 1008–1021. [Google Scholar] [CrossRef]
  21. Yousefi, R.; Muhamad, M.R. Effects of gold catalysts and thermal evaporation method modifications on the growth process of Zn1-x MgxO nanowires. J. Solid State Chem. 2010, 183, 1733–1739. [Google Scholar] [CrossRef]
  22. Moghaddam, M.; Pourbaige, M.; Tabar, H.K.; Farhadi, N.; Hosseini, S.M.A. Composition and Antifungal Activity of Peppermint (Mentha piperita) Essential Oil from Iran. J. Essent. Oil Bear. Plants 2013, 506–512. [Google Scholar] [CrossRef]
  23. Moghaddam, T.B.; Soltani, M.; Karim, M.R. Evaluation of permanent deformation characteristics of unmodified and Polyethylene Terephthalate modified asphalt mixtures using dynamic creep test. Mater. Des. 2014, 53, 317–324. [Google Scholar] [CrossRef] [Green Version]
  24. Li, J.M.; Zhang, Y.J. GIS based information model for assessment of geothermal potential: Case study of Tengchong county, southwest China. J. Eng. Geol. 2018, 26, 540–550. [Google Scholar] [CrossRef]
  25. Lachenbruch, A.H. Heat flow in the Basin and Range province and thermal effects of tectonic extension. Pure Appl. Geophys. 1978, 117, 34–50. [Google Scholar] [CrossRef]
  26. Vlassova, L.; Perezcabello, F.; Nieto, H.; Martín, P.; Riano, D. Assessment of methods for land surface temperature retrieval from Landsat-5 TM images applicable to multiscale tree-grass ecosystem modeling. Remote Sens. 2014, 6, 4345–4368. [Google Scholar] [CrossRef] [Green Version]
  27. Johnson, L.R. Source mechanisms of induced earthquakes at the Geysers geothermal reservoir. Pure Appl. Geophys. 2014, 171, 1641–1668. [Google Scholar] [CrossRef] [Green Version]
  28. Mcguire, J.J.; Lohman, R.B.; Catchings, R.D.; Rymer, M.J.; Goldman, M.R. Relationships among seismic velocity, metamorphism, and seismic and aseismic fault slip in the Salton Sea Geothermal Field region. Geophys. Res. Solid Earth 2015, 120, 2600–2615. [Google Scholar] [CrossRef] [Green Version]
  29. Crespo, E.; Lillo, J.; Oyarzun, R.; Cubas, P.; Leal, M. The Mazarron basin, SE Spain: A study of mineralization processes, evolving magmatic series, and geothermal activity. Int. Geol. Rev. 2013, 55, 1978–1990. [Google Scholar] [CrossRef]
  30. Kang, S.C.; Xu, Y.W.; You, Q.L.; Flugel, W.A.; Pepin, N.; Yao, T.D. Review of climate and cryospheric change in the Tibetan Plateau. Env. Res. Lett. 2010, 5, 8. [Google Scholar] [CrossRef]
  31. Yang, D.H. Analysis of Main Engineering Geological Problems in Changdu to Linzhi Section of Sichuan-Tibet Railway. Railw. Stand. Des. 2019, 63, 1622–1627. [Google Scholar] [CrossRef]
  32. Vlassova, L.; Perezcabello, F.; Mimbrero, M.; Llovería, R.M.; García-Martín, A. Analysis of the relationship between land surface temperature and wildfire severity in a series of landsat images. Remote Sens. 2014, 6, 6136–6162. [Google Scholar] [CrossRef] [Green Version]
  33. Jiménez-Muñoz, J.C.; Sobrino, J.A. A generalized single-channel method for retrieving land surface temperature from remote sensing data. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  34. Becker, F.; Li, Z.L. Temperature-independent spectral indices in thermal infrared bands. Remote Sens. Environ. 1990, 32, 17–33. [Google Scholar] [CrossRef]
  35. Sobrino, J.; Jiménez-Muñoz, J.C.; Paolini, L. Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440. [Google Scholar] [CrossRef]
  36. Wan, Z.M.; Li, Z.L. A physics-based algorithm for retrieving land-surface emissivity and temperature from EOS/MODIS data. IEEE Trans. Geosci. Remote Sens. 1997, 35, 980–996. [Google Scholar] [CrossRef]
  37. Wang, N.; Wu, H.; Nerry, F.; Li, C.R.; Li, Z.L. Temperature and emissivity retrievals from hyperspectral thermal infrared data using linear spectral emissivity constraint. IEEE Trans. Geosci. Remote Sens. 2001, 49, 1291–1303. [Google Scholar] [CrossRef]
  38. Qin, Z.H.; Zhang, M.H.; Arnon, K. Mono-window Algorithm for Retrieving Land Surface Temperature from Landsat TM6 data. Acta Geogr. Sin. 2001, 456–466. [Google Scholar] [CrossRef]
  39. Jin, D.D.; Gong, Z.N. Algorithms Comparison of Land Surface Temperature Retrieval from Landsat Series Data: A Case Study in Qiqihar, China. Remote Sens. Technol. Appl. 2018, 33, 830–841. [Google Scholar] [CrossRef]
  40. Zhao, P.D.; Chi, S.D. Preliminary view on geological anomaly. Earth Sci. J. China Univ. Geosci. 1991, 16, 241–248. [Google Scholar]
  41. Chi, S.D.; Zhao, P.D. Application of Combined-entropy Anomany of Geological Formations to Delineation of Preferable Ore-finding Area. Geoscience 2000, 423–428. [Google Scholar] [CrossRef]
  42. Revil, A.; Pezard, P.A. Streaming electrical potential anomaly along faults in geothermal areas. Geophys. Res. Lett. 1998, 25, 3197–3200. [Google Scholar] [CrossRef]
  43. Zhang, F.W.; Wang, G.L.; Hou, X.W.; Li, J.H.; Li, Y.J. An Analysis of the Formation of Geothermal Resources and the Effects of Groundwater Circulation on the Wall Rock Temperature Field-Taking the Pingdingshan Mining Field as an Example. Acta Geosci. Sin. 2000, 21, 142–146. [Google Scholar]
  44. Hochstein, M.P.; Hunt, T.M. Seismic, gravity and magnetic studies, Broadlands geothermal field. Geothermics 1970, 2, 333–346. [Google Scholar] [CrossRef]
  45. Institute of Geography; CAS. Atlas of Qinghai-Tibet Plateau; Science Press: Beijing, China, 2015. [Google Scholar]
  46. Aydogan, D. Extraction of lineaments from gravity anomaly maps using the gradient calculation: Application to Central Anatolia. Earth Planets Space 2011, 63, 903–913. [Google Scholar] [CrossRef] [Green Version]
  47. Wang, D.F. Aeromagnetic Series Maps and Specifications of Qinghai-Tibet Plateau and Adjacent Areas; Geological Publishing House: Beijing, China, 2013. [Google Scholar]
  48. Gutenberg, B.; Richter, C.F. Frequency of earthquakes in California. Bull. Seismol. Soc. Am. 1945, 34, 185–188. [Google Scholar] [CrossRef] [Green Version]
  49. Bayrak, Y.; Bayrak, E. Regional variations and correlations of Gutenberge-Richter parameters and fractal dimension for the different seismogenic zones in Western Anatolia. J. Asian Earth Sci. 2012, 58, 98–107. [Google Scholar] [CrossRef]
  50. Tang, G.A.; Yang, X. Experimental Course of Spatial Analysis of ArcGIS; Science Press: Beijing, China, 2004. [Google Scholar]
  51. Zhao, H.; Song, E.X. Improved Information Value Model and Its Application in the Spatial Prediction of Landslides. J. Civ. Archit. Environ. Eng. 2011, 33, 38–44. [Google Scholar] [CrossRef]
  52. Zhang, X.M. Hazard Assessment and Zoning Research of Landslide in Shanxi Province Based on GIS; Chang’an University: Xi’an, China, 2019. [Google Scholar]
  53. Ye, C.F.; Zhang, Y.C.; Xiong, J.N.; Qin, J.X. Hazard Assessment of Mountain Torrent Disaster in Small Watersheds of the Hilly Areas of Hunan Province. J. Geo-Inf. Sci. 2017, 19, 1593–1603. [Google Scholar] [CrossRef]
  54. Tziritis, E.; Lombardo, L. Estimation of intrinsic aquifer vulnerability with index-overlay and statistical methods: The case of eastern Kopaida, central Greece. Appl. Water Sci. 2016. [Google Scholar] [CrossRef] [Green Version]
  55. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; The University of Illinois Press: Urbana, IL, USA, 1947. [Google Scholar] [CrossRef]
  56. Tangestani, M.H.; Moore, F. Porphyry copper potential mapping using the weights-of- evidence model in a GIS, northern Shahr-e-Babak, Iran. Aust. J. Earth Sci. 2002, 48, 695–701. [Google Scholar] [CrossRef]
  57. Lee, S.; Talib, J.A. Probabilistic landslide susceptibility and factor effect analysis. Environ. Geol. 2005, 47, 982–990. [Google Scholar] [CrossRef]
  58. Cong, W.Q.; Li, T.F. Study on the key problems of landslide and debris flow hazard zoning based on GIS. Earth Sci. Front. 2006, 13, 185–190. [Google Scholar] [CrossRef]
  59. Xiao, K.Y.; Zhang, X.H.; Chen, Z.H. Comparison of evidence weight method and information quantity method in metallogenic prediction. Calc. Technol. Geophys. Geochem. Explor. 1999, 21, 223–226. [Google Scholar] [CrossRef]
  60. He, Z.C.; Chi, C.X. Effects of spatial resolution and sampling methods on watershed feature extraction from DEM. J. Glaciol. Geocryol. 2004, 26, 610–616. [Google Scholar] [CrossRef]
  61. Dahal, R.K.; Hasegawa, S.; Bhandary, N.P.; Poudel, P.P.; Nonomura, A.; Yatabe, R. A replication of landslide hazard mapping at catchment scale. Geomat. Nat. Hazards Risk 2012, 3, 161–192. [Google Scholar] [CrossRef]
  62. Romero-Calcerrada, R.; Novillo, C.J.; Millington, J.D.A.; Gomez-Jimenez, I. GIS analysis of spatial patterns of human-caused wildfire ignition risk in the SW of Madrid (Central Spain). Landsc. Ecol. 2008, 23, 341–354. [Google Scholar] [CrossRef]
Figure 1. The technology and methodology flow chart.
Figure 1. The technology and methodology flow chart.
Remotesensing 13 01606 g001
Figure 2. Geographical location and topography of the study area. (a) Location of the study area; (b) topography of the study area (Landsat8 true color synthetic overlaid DEM; spatial resolution: 30 m; acquisition time: 2018.11–2019.1).
Figure 2. Geographical location and topography of the study area. (a) Location of the study area; (b) topography of the study area (Landsat8 true color synthetic overlaid DEM; spatial resolution: 30 m; acquisition time: 2018.11–2019.1).
Remotesensing 13 01606 g002
Figure 3. Field temperature measurement of geothermal points.
Figure 3. Field temperature measurement of geothermal points.
Remotesensing 13 01606 g003
Figure 4. Factor maps.
Figure 4. Factor maps.
Remotesensing 13 01606 g004
Figure 5. The calculation of the combined entropy of geological units. (a) Surface pattern in two dimensions; (b) surface lithological identification—lithological units are identified by color and pattern, S1 and S3 are surfaces belonging to the same lithological unit, S2 and S4 correspond to another lithological unit; (c) the calculation of combined entropy.
Figure 5. The calculation of the combined entropy of geological units. (a) Surface pattern in two dimensions; (b) surface lithological identification—lithological units are identified by color and pattern, S1 and S3 are surfaces belonging to the same lithological unit, S2 and S4 correspond to another lithological unit; (c) the calculation of combined entropy.
Remotesensing 13 01606 g005
Figure 6. Factor reclassification maps.
Figure 6. Factor reclassification maps.
Remotesensing 13 01606 g006aRemotesensing 13 01606 g006b
Figure 7. Classification maps of prediction.
Figure 7. Classification maps of prediction.
Remotesensing 13 01606 g007
Figure 8. Curves of prediction ratio function and its area ratio.
Figure 8. Curves of prediction ratio function and its area ratio.
Remotesensing 13 01606 g008
Figure 9. Verification point location and surrounding environment. (a) Location of verification points, (b) Surrounding environment of verification point 6.
Figure 9. Verification point location and surrounding environment. (a) Location of verification points, (b) Surrounding environment of verification point 6.
Remotesensing 13 01606 g009
Table 1. Reclassification threshold for impact factor maps.
Table 1. Reclassification threshold for impact factor maps.
Impact FactorsClass
1234567
LST (°C)−10.00–6.21−6.21–0.690.69–6.836.83–11.7511.75–16.0216.02–23.5423.54–50.42
Combined entropy of geological formation0–1010–2020–4040–5050–6060–8080–100
Buffer distance to fault (km)0–22–44–66–88–1010–12>12
Fault density (km/km2)0.000–0.0190.019–0.0470.047–0.0690.069–0.0920.092–0.1170.117–0.1480.148–0.266
Buffer distance to river (km)0–22–44–66–88–1010–12>12
Magnetic anomaly (nT)−241–−163−163–−115−115–−75−75–−50−50–−26−26–88–100
Bouguer gravity (mgal)−540–−520−520–−480−480–−440−440–−400−400–−340−340–−300−300–−220
Earthquake peak acceleration (g).0.02–0.080.08–0.140.14–0.200.20–0.240.24–0.300.30–0.340.34–0.40
Table 2. Grid size selection standard.
Table 2. Grid size selection standard.
Study Area (km2)ScaleGrid Size (m)
≥100,000≤1:250,000≥100 × 100
10,000–100,0001:250,000–1:100,000≥50 × 50
1000–10,0001:100,000–1:50,000≥25 × 25
<10,000>1:50,000≥5 × 5
Table 3. Correlation of impact factor maps.
Table 3. Correlation of impact factor maps.
Impact FactorsLSTCombined Entropy Buffer Distance to FaultFault DensityBuffer Distance to RiverMagnetic AnomalyBouguer GravityEarthquake Peak Acceleration
LST1.000−0.0620.127−0.1520.0380.0530.089−0.006
Combined entropy −0.0621.000−0.1580.175−0.1570.035−0.079−0.037
Buffer distance to fault0.127−0.1581.000−0.7710.0640.171−0.0210.263
Fault density −0.1520.175−0.7711.000−0.085−0.160−0.031−0.246
Buffer distance to river0.038−0.1570.064−0.0851.000−0.0340.1500.051
Magnetic anomaly0.0530.0350.171−0.160−0.0341.000−0.406−0.016
Bouguer gravity0.089−0.079−0.021−0.0310.150−0.4061.0000.292
Earthquake peak
acceleration
−0.006−0.0370.263−0.2460.051−0.0160.2921.000
Table 4. Information values for impact factors.
Table 4. Information values for impact factors.
Impact factorsClassPointsGridsIijHiWi
LST127,459,893−3.1450.8032.462
226,526,266−3.012
3266,864,561−0.497
4557,272,7700.194
5856,157,5220.796
6513,924,6830.736
7281,781,7930.926
Combined entropy of geological formation15321,710,503−0.9360.9550.568
2182,378,6330.195
3424,703,7930.360
4262,538,3500.498
5372,590,8930.830
6504,408,9540.599
7231,656,3620.802
Fault density1405,387,5180.1760.9320.847
2554,568,4250.659
3343,905,2750.335
4443,325,4380.754
5332,855,4770.618
6302,452,8990.675
71317,492,456−2.126
Buffer distance to river11125,920,2661.1110.8691.647
2365,608,9090.030
3225,273,756−0.401
4204,772,597−0.396
5184,150,553−0.362
6133,479,987−0.511
72810,781,420−0.875
Magnetic anomaly151,393,066−0.5510.9640.457
292,079,263−0.364
3205,555,534−0.548
4709,336,8000.186
58811,748,2980.185
6527,902,4640.055
751,972,063−0.899
Earthquake peak acceleration174,149,900−1.3060.9800.254
2136,921,418−1.199
36913,669,294−0.210
4839,371,5860.352
5523,792,1040.789
6231,108,6621.203
72974,524−1.110
Table 5. Values of the weights of evidence for impact factors.
Table 5. Values of the weights of evidence for impact factors.
Impact FactorsClass W ij + W ij CijCi
LST1−0.3460.198−0.5443.547
2−0.3120.170−0.482
3−0.4980.078−0.576
40.194−0.0490.243
50.896−0.3501.246
61.535−0.2261.761
71.725−0.1741.899
Combined entropy of geological formation1−0.9240.529−1.4532.203
20.195−0.0140.208
30.360−0.0600.420
40.497−0.0450.542
50.830−0.0940.924
60.599−0.1070.706
70.802−0.0550.856
Fault density1−1.0210.416−1.4373.383
2−0.1600.050−0.211
30.283−0.0380.321
40.586−0.0910.677
50.978−0.0911.069
61.310−0.0951.406
71.515−0.0421.557
Buffer distance to river11.111−0.7371.8481.324
20.730−0.6501.380
3−0.1010.049−0.150
4−0.0960.043−0.140
5−0.2620.035−0.297
6−0.4110.137−0.549
7−0.5740.195−0.768
Magnetic anomaly1−0.2510.015−0.2660.862
2−0.1640.017−0.181
3−0.4480.036−0.484
40.685−0.4641.149
50.707−0.3971.104
60.455−0.3140.769
7−0.6990.530−1.229
Earthquake peak acceleration1−1.4430.089−1.5321.183
2−0.5560.128−0.684
30.063−0.0340.097
40.290−0.1090.398
50.576−0.1850.761
60.623−0.2780.901
70.905−0.3371.242
Table 6. Threshold of model classification values.
Table 6. Threshold of model classification values.
Model Class
1234
Index-overlay information model−1.99–0.000.00–0.240.24–0.570.57–1.32
Weights of entropy information model−16.21–0.000.00–6.506.50–7.217.21–8.81
Weights of evidence information model−31.72–0.000.00–3.843.84–8.418.41–18.94
Table 7. Success indices from the classification maps of prediction.
Table 7. Success indices from the classification maps of prediction.
ClassIndex-OverlayWeights of EntropyWeights of Evidence
GridsPointsIndexGridsPointsIndexGridsPointsIndex
41,689,379740.44%2,696,6871060.39%1,814,461960.53%
33,764,914770.21%5,649,918840.15%4,127,875750.18%
24,397,731520.12%5,009,384310.06%5,189,797410.08%
129,611,165460.02%26,107,200280.01%28,331,056370.01%
Table 8. Verification point details.
Table 8. Verification point details.
NumberMeasurement DateCoordinateAltitude (m)Temperature (°C)
112 June 201994°44′19″ E,317530.8–39.7
29°48′56″ N
212 June 201995°02′22″ E,203350.3–64.6
30°03′54″ N
59 June 2019100°08′29″ E,397131.4–52.3
30°02′16″ N
1012 June 201996°50′44″ E,437719.3–33.2
30°42′49″ N
119 June 201999°20′24″ E,329717.9–23.9
30°16′55″ N
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, W.; Dong, Q.; Chen, Z.; Feng, T.; Wang, D.; Jiang, L.; Du, S.; Zhang, X.; Meng, D.; Bian, M.; et al. Weighted Information Models for the Quantitative Prediction and Evaluation of the Geothermal Anomaly Area in the Plateau: A Case Study of the Sichuan–Tibet Railway. Remote Sens. 2021, 13, 1606. https://doi.org/10.3390/rs13091606

AMA Style

Zhao W, Dong Q, Chen Z, Feng T, Wang D, Jiang L, Du S, Zhang X, Meng D, Bian M, et al. Weighted Information Models for the Quantitative Prediction and Evaluation of the Geothermal Anomaly Area in the Plateau: A Case Study of the Sichuan–Tibet Railway. Remote Sensing. 2021; 13(9):1606. https://doi.org/10.3390/rs13091606

Chicago/Turabian Style

Zhao, Wenbo, Qing Dong, Zhe Chen, Tao Feng, Dong Wang, Liangwen Jiang, Shihui Du, Xiaoyu Zhang, Deli Meng, Min Bian, and et al. 2021. "Weighted Information Models for the Quantitative Prediction and Evaluation of the Geothermal Anomaly Area in the Plateau: A Case Study of the Sichuan–Tibet Railway" Remote Sensing 13, no. 9: 1606. https://doi.org/10.3390/rs13091606

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop