Study on Feng Shui (Geomantic) Suitability Evaluation of Mausoleums in Nanjing City Based on GIS

Taking the mausoleums in the main urban area of Nanjing city as the study object, with the support of GIS technology, this paper extracts eight influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) corresponding to geomantic elements (“Long”, “Sha”, “Shui”, and “Xue”) in the geographical environments of mausoleums. The AHP is applied in the correlation weight method to calculate the weights of the eight influencing factors, while the information content method is applied in the grading method to calculate the numerical index of each factor. The feng shui (geomantic) suitability of the mausoleums in the main urban area of Nanjing is evaluated using the AHP-weighted information content method and by combining expert knowledge and experience with mathematical statistical techniques. According to the feng shui (geomantic) suitability evaluation results obtained for the mausoleums in the study area, the relatively low-suitability area accounts for 14.55% of the entire study area, the low-suitability area accounts for 25.40%, the suitable area accounts for 29.13%, the relatively high-suitability area accounts for 22.00%, and the high-suitability area accounts for 8.92%. Finally, through a sensitivity analysis of the model, a verification analysis applied using random mausoleums and a feng shui (geomantic) analysis of the Ming Xiaoling Mausoleum, the feng shui (geomantic) suitability evaluation results obtained for the mausoleums analyzed herein are found to be consistent with the actual distribution of the mausoleum sites. The research shows that the AHP-weighted information method based on GIS is suitable for evaluating the feng shui (geomantic) suitability of mausoleums.


Introduction
Tombs are places where the remains of the dead are stored and include architectural shrines, burial vaults, church monuments, mausoleums and so on [1,2]. In China, mausoleums are traditional buildings that refer to the burial places of emperors, princes, and prominent figures. The design, site selection, and layout of mausoleums are all deeply influenced by feng shui [3][4][5]. Feng shui is the Chinese art or practice of positioning objects, especially graves, buildings, and furniture. Furthermore, in Chinese philosophy, feng shui is regarded as the rule governing spatial arrangements and orientations [6,7]. Therefore, studying the feng shui of mausoleums can deepen the understanding of the spatial arrangements of mausoleums and the scientific understanding of geomancy.
Analyses of the geomancy of mausoleums mainly involve the four geomantic elements of "Long", "Sha", "Shui", and "Xue" to judge spatial arrangements [8,9]. The judgments involved in geomancy must be combined with an analysis of the surrounding geographical environment [10][11][12][13]. "Long" is a mountain range and refers to the back of a mausoleum. The analysis of the mountain potential and the overall elevation of the mountains are important conditions to consider when placing a mausoleum. "Sha" refers to the landscape surrounding the mausoleum and is used to "gather Qi". Analyzing the mountains before and after the mausoleum is placed and visualizing the landscape of the surrounding image are other conditions to consider when placing a mausoleum. "Shui" is the outfall and river, and the search for an outfall and the observation of the flow form are conditions for the mausoleum placement. "Xue" is the specific location of the mausoleum within the surrounding environment. In a subsidence environment, the point at which the most "Qi" can be gathered is chosen [14][15][16]. At present, the methods used to evaluate the feng shui (geomantic) suitability of mausoleums are mainly qualitative [17][18][19], including the methods used by geomantic experts to obtain geomantic assessment results in the region of study through empirical judgments of the geomantic elements ("Long", "Sha", "Shui", and "Xue") in combination with the geographical environment in the field. These assessment methods rely on expert geomantic experience to a high degree [20,21].
The geomantic elements ("Long", "Sha", "Shui", and "Xue") correspond to the influencing factors (surface cutting depth, elevation, surface roughness, surface peaks, outfall, river, slope, aspect) of the geographical environment. Mausoleums, as a kind of geographical and cultural heritage site, have been studied deeply by many scholars in recent years [22][23][24]. The methods used to research mausoleums include not only traditional archaeological techniques but also spatial information technologies with GIS technologies as the core [25,26]. In terms of research content, previous studies have included not only spatial analyses of heritage sites and evolution pattern analyses of heritage sites and their surrounding environments but also some analyses of heritage risk identification and risk assessments [27,28]. Therefore, GIS technology can be effectively applied to evaluate the feng shui (geomantic) suitability of the geographical environment of a mausoleum based on the powerful associated spatial data acquisition, analysis and processing functions, and the complexity of the geographical environment of the studied mausoleum.
Therefore, geomantic mausoleum judgments can be regarded as suitability assessments of the geographical environments surrounding mausoleums [29,30]. The methods used to conduct suitability assessments of these geographical environments mainly include the correlation weight method of each influencing factor and the grading method of a single influencing factor. The correlation weight method of each influencing factor refers to the method in which the correlation of each influencing factor and the relationship between each influencing factor and the geographic environment are first determined and then the geographic environment is assessed by calculating the grading differences in the obtained spatial characteristics. These methods include the weight of evidence analysis method [31,32], logistic regression method [33,34], fuzzy comprehensive evaluation method [35,36], and analytic hierarchy process [37,38]. The weight of evidence analysis method examines the conditional independence among different influencing factors, calculates the weight coefficient of each influencing factor, and then conducts a statistical synthesis on the influencing factor layers [39,40]; in the weight evidence analysis method, it is difficult to resolve the influence of non-independence conditions on the evaluation results, and the correlations among influencing factors are ignored [41][42][43]. The logistic regression method has limited application conditions, as it requires a linear relationship between influencing factors and logistic probability, and the degree of fit and accuracy are not high when performing calculations; further, this method ignores the complex relationships among the geographical environment factors of the analyzed mausoleum [44,45]. The fuzzy comprehensive evaluation method uses the degree of membership to describe the fuzzy boundary and comprehensively evaluates the membership level of the evaluated object from multiple influencing factors [46,47]. However, this method involves complicated calculations, and the determination of the influencing factor weights is highly subjective and ignores the correlations among influencing factors [48,49]. In contrast, the analytic hierarchy process takes the research object as a system and analyzes the system in a comprehensive thinking process. This method requires less quantitative data and information and has become a simple and practical method used to assess the geographic environments of mausoleums [50,51].
The grading method of single influencing factor refers to the quantitative evaluation method in which the classification of each influencing factor is based on the existing geographic object; this method is performed under the assumption that no correlations exist among the influencing factors. This method reveals the degree of the geographical environment by calculating the spatial characteristic value of each influencing factor and reveals the relationship between the classification of each influencing factor and the geographical environment. These methods include the frequency ratio method [52,53] and information content method [54,55]. The frequency ratio method first extracts the point density and areal density of the analyzed mausoleums, uses the frequency ratio to calculate the probability using the influencing factors, and finally superimposes the obtained data to obtain the evaluation results [56,57]. However, for mausoleum data, it is difficult to obtain accurate area information and the reliability of the calculated results is therefore low; additionally, this method ignores the complexity of the geographical environment [58,59]. The information content method calculates the spatial characteristics of each influencing factor based on the information content representing the given assessment level of the target mausoleum's geographic environment and uses these characteristics to measure the intimacy of each influencing factor with the mausoleum's geographic environment [60][61][62]. Therefore, under the assumption that no correlations exist among the various influencing factors, the information content method is more appropriate than the frequency ratio method for determining the geographic positions of mausoleums.
In the methods in which the correlation weight of each influencing factor is analyzed, the knowledge and experience of geomantic experts can be applied to the studied mausoleums to determine the correlation between each influencing factor and the relationship between each influencing factor and the geomancy of the mausoleums [63,64]. However, nonquantitative and empirical expert judgments can easily affect the evaluation results, and the objectivity of the assessment results cannot be tested. Methods in which individual influencing factors are graded can quantitatively calculate the relationship between a single influencing factor and the geomancy of mausoleums based on geographical environment data characterizing the analyzed mausoleums, thus avoiding the subjectivity of the assessment index weight value [65,66]. However, this method tends to ignore the correlations between influencing factors. Therefore, this paper proposes a method in which the correlation weight method of each influencing factor and the grading method of a single influencing factor are combined to assess the geomancy of mausoleums. This method not only considers the knowledge and experience of feng shui experts, but also considers the correlation among the influencing factors, so that the feng shui evaluation results of mausoleum are more scientific.
Based on GIS technology, this paper takes the main urban area of Nanjing city, Jiangsu Province as the study area and extracts eight types (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) of geomantic mausoleum assessment factors by combining: geomancy theory; the analytical hierarchy process (AHP) is applied in the correlation weight method of each influencing factor to calculate the weights of the eight assessment factors; the information content method (I) is applied in the grading method of a single influencing factor to calculate the numerical index of each factor; a comprehensive weighting method is used to construct a geomancy assessment model of mausoleums (AHP-I); and the weight of each factor is optimally combined with a single-factor grading numerical index to assess the geomancy of mausoleums in the study area. tain, and Qixia Mountain in the western Jiangning County. The study area is dominated by low mountains and gentle hills, with large undulations and complex and changeable topography. Dense rivers and abundant water resources are found in the territory. Nanjing is the ancient capital of the Six Dynasties and has a long history; it has long been the political, economic, and cultural center of Southern China. Emperors, princes, and prominent figures are buried in this city. The main urban area of Nanjing mainly includes eight areas: Pukou District, Qixia District, Gulou District, Jianye District, Xuanwu District, Qinhuai District, Yuhuatai District, and Jiangning District.

Research Area and
According to the cultural relic protection unit information published by China's State Administration of Cultural Heritage, this paper selected the mausoleums in national, provincial, and municipal cultural relic protection units within the study area as research objects (42 mausoleums in total). These mausoleums were divided into categories such as emperor mausoleums, dukes under an emperor mausoleum, and general ministers' mausoleums, as shown in Table 1 and Figure 1. The spatial information and attributes of these mausoleums were collected and entered into a database. In the mausoleum database, each mausoleum is represented by a data point, and its corresponding attribute information includes its name, longitude, latitude, dynasty, and the corresponding spatial attribute values of 8 influencing factors, as shown in Figure 2. The spatial distribution of mausoleums in the study area is shown in Figure 3.

Selection of Assessment Factors and Data Sources
Geomantic elements ("Long", "Sha", "Shui", and "Xue") correspond to influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) in the geographical environment, and the judgment of mausoleum geomancy can be regarded as an assessment of the suitability of the mausoleum environment.
"Long", in feng shui theory (hypothesis), refers to mountains and the trends of mountains, and mainly corresponds to feng shui suitability judgments through the observation of the overall situation of a large mountain range in the geographical environment and the specific situation of the surrounding mountains. In geography, the surface roughness is an index reflecting the impact of the land surface on sand activities and its weakening impact on the wind speed, as well as the ability of the land surface to resist wind erosion. According to the surface roughness, the wind resistance of the geographical environment can be judged. Strong wind resistance in a given area can ensure that the "Qi" of the local environment is not dispersed and that the feng shui suitability is high. On the other hand, the elevation of a given area can express the relative height of the local mountains. Therefore, the surface roughness and elevation metrics can be used as the basic criterion for analyzing "Long" in Feng Shui.
"Sha", according to feng shui theory (hypothesis), refers to the low and gentle hills around a mausoleum. Its function is to let the "Long" surrounding the overall environment cover the strong wind and protect the "Qi" of the small environment to achieve wind accumulation. In geography, the surface cutting depth refers to the erosion degree of the land surface by natural factors (such as water, wind, and sand). The degree of surface cutting can represent the hill or shallow basin density. The relief of hills can be judged using the hill density. On the other hand, the peak point of a local area can be extracted to more intuitively judge the locations of hills. Through observations of the locations of hills, it can be determined whether the hills around and in front of a mausoleum are surrounded by similar topography in the geographical environment. Therefore, the analyses of the surface cutting depth and peak point can be used as the basic criterion of "Sha" in feng shui applications.

Selection of Assessment Factors and Data Sources
Geomantic elements ("Long", "Sha", "Shui", and "Xue") correspond to influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) in the geographical environment, and the judgment of mausoleum geomancy can be regarded as an assessment of the suitability of the mausoleum environment.
"Long", in feng shui theory (hypothesis), refers to mountains and the trends of mountains, and mainly corresponds to feng shui suitability judgments through the observation of the overall situation of a large mountain range in the geographical environment and the specific situation of the surrounding mountains. In geography, the surface roughness is an index reflecting the impact of the land surface on sand activities and its weakening impact on the wind speed, as well as the ability of the land surface to resist wind erosion. According to the surface roughness, the wind resistance of the geographical environment can be judged. Strong wind resistance in a given area can ensure that the "Qi" of the local environment is not dispersed and that the feng shui suitability is high. On the other hand, the elevation of a given area can express the relative height of the local mountains. Therefore, the surface roughness and elevation metrics can be used as the basic criterion for analyzing "Long" in Feng Shui.
"Sha", according to feng shui theory (hypothesis), refers to the low and gentle hills around a mausoleum. Its function is to let the "Long" surrounding the overall environment cover the strong wind and protect the "Qi" of the small environment to achieve wind accumulation. In geography, the surface cutting depth refers to the erosion degree of the land surface by natural factors (such as water, wind, and sand). The degree of surface cutting can represent the hill or shallow basin density. The relief of hills can be judged using the hill density. On the other hand, the peak point of a local area can be extracted to more intuitively judge the locations of hills. Through observations of the locations of hills, it can be determined whether the hills around and in front of a mausoleum are surrounded by similar topography in the geographical environment. Therefore, the analyses of the surface cutting depth and peak point can be used as the basic criterion of "Sha" in feng shui applications.
In feng shui theory (hypothesis), "Shui" mainly refers to the outfall and shape of water in a local environment. The outfall is the outlet of flowing water in a basin and the lowest point of the whole basin. The feng shui suitability can be judged according to the distance from the outfall. The basic requirement for the water shape in areas with high feng shui suitability is "buckling embracing". The feng shui suitability can thus be judged according to the distance from the dendritic river system. Therefore, the use of the two geographical environmental factors of the outfall and river system can serve as the basic criterion of "Shui" in feng shui applications.
"Xue", in feng shui theory (hypothesis), refers to the determination of the specific location of a mausoleum. "Xue" refers to finding the most "Qi-gathering" position in the subsidence environment to build a mausoleum. The premise of "gathering Qi" is that "Qi" should be transported through the rolling "Long". Therefore, "Xue" and "Long" are closely related. The influence of mountain slopes and aspect on feng shui suitability should be considered. Aspect is the direction facing the slope, and its numerical value reflects the projection direction of the sun on the surface. The aspect can be divided into shaded slopes and sunny slopes, and different solar directions have different influences on the location of mausoleums. Therefore, the consideration of known mausoleum points, sunny slopes, and shaded slopes can be used as the basic criterion of "Xue" in feng shui applications.
In summary, these eight influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) can be used as geomantic suitability evaluation factors for mausoleums. The sources of the influencing factor data are shown in Table 2. The corresponding table of geomantic elements and influencing factors is shown in Table 3.     "Long" Indicates the trends of mountains and mountain ranges.

Surface roughness
Indicates the weakening effect of the land surface on the wind speed and the windsand interactive activity. Can be used to describe the "Long" component in the feng shui of mausoleums Elevation Indicates the relative height of a mountain.
Can be used to describe the "Long" component in the feng shui of mausoleums Surface cutting depth Indicates the density of hills. Can be used to describe the "Sha" component in the feng shui of mausoleums Table 3. Cont.

Geographical Definition References [67-70] Graphs
Graphs References [8,9,71] "Sha" Indicates low hills and their forms around and in front of a mausoleum

Surface cutting depth
Indicates the density of hills.
Can be used to describe the "Sha" component in the feng shui of mausoleums "Long" Indicates the trends of mountains and mountain ranges.
tivity. Can be used to describe the "Long" component in the feng shui of mausoleums Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval; "Long" Indicates the trends of mountains and mountain ranges. Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval;

Surface peaks
Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Long" of mountains and mountain ranges. Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval; "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

Elevation
Indicates the relative height of a mountain. Can be used to describe the "Long" component in the feng shui of mausoleums "Sha" Indicates low hills and their forms around and in front of a mausoleum Surface cutting depth Indicates the density of hills. Can be used to describe the "Sha" component in the feng shui of mausoleums Surface peaks Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval; mountain ranges. ums Elevation Indicates the relative height of a mountain. Can be used to describe the "Long" component in the feng shui of mausoleums "Sha" Indicates low hills and their forms around and in front of a mausoleum Surface cutting depth Indicates the density of hills. Can be used to describe the "Sha" component in the feng shui of mausoleums Surface peaks Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval;

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin

Elevation
Indicates the relative height of a mountain.
Can be used to describe the "Long" component in the feng shui of mausoleums "Sha" Indicates low hills and their forms around and in front of a mausoleum Surface cutting depth Indicates the density of hills. Can be used to describe the "Sha" component in the feng shui of mausoleums Surface peaks Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum Indicates an area higher than the surrounding points within the analysis scope of a specific field. Can be used to describe the "Sha" component in the feng shui of mausoleums "Shui" Indicates outfall and water shape.

Outfall
Indicates the outlet of flowing water in the basin, corresponding to the lowest point of the whole basin

River system
Indicates the water network system composed of all rivers, lakes, and other water bodies in the basin "Xue" Indicates the specific location of the mausoleum

Slope
Indicates the degree of steepness and gentleness of surface units. For mausoleum points, a buffer zone is established by taking the slope value as the distance interval; "Long" Indicates the trends of mountains and mountain ranges.
describe the "Long" component in the feng shui of mausoleums Indicates the direction of mountains from high to low. Extracting the slope direction of mausoleums can distinguish between sunny and shaded slopes; this information can be used to describe the "Xue" component in the feng shui of mausoleums

AHP-Weighted Information Content Method
In this paper, a weighted AHP information method is proposed to evaluate the feng shui suitability of mausoleums. This weighted method comprehensively evaluates the influencing factors of various geographical environments as well as the contribution degree within the quantified classification range from the overall and partial dimensions. Its essence is to quantify the cognitive evaluation of the geographical environment in terms of feng shui through the use of GIS mathematical statistics.
(1) First, the weight value of the target geographical environment influencing factor is calculated by the AHP method to characterize the importance of the influencing factor to the feng shui suitability of the analyzed mausoleum. The eight influencing factors are arranged, their relative importance is compared, and a judgment matrix is constructed. The matrix is shown in (Equation 1).

AHP-Weighted Information Content Method
In this paper, a weighted AHP information method is proposed to evaluate the feng shui suitability of mausoleums. This weighted method comprehensively evaluates the influencing factors of various geographical environments as well as the contribution degree within the quantified classification range from the overall and partial dimensions. Its essence is to quantify the cognitive evaluation of the geographical environment in terms of feng shui through the use of GIS mathematical statistics.
(1) First, the weight value of the target geographical environment influencing factor is calculated by the AHP method to characterize the importance of the influencing factor to the feng shui suitability of the analyzed mausoleum.
The eight influencing factors are arranged, their relative importance is compared, and a judgment matrix is constructed. The matrix is shown in (Equation (1)).
The weight calculation formula used to determine the weights of various factors influencing the mausoleum geographical environment is shown as follows (Equation (2)): (2) Second, the information value of each geographical environment influencing factor in its different classification range is calculated by the information content method to represent the contribution degree of an influencing factor in various classifications.
The calculation formula of the information method is shown as follows (Equation (3)): In Equation (3), T i is the influencing factor (x i ) corresponding to the number of mausoleums distributed in each level unit; T is the total number of tombs in the study area; S i is the influencing factor (x i ) corresponding to the total number of grid cells; and S is the total number of grid cells in the study area.
(3) Finally, the comprehensive weighted method is used to calculate the compound operation to obtain the feng shui evaluation results of the analyzed mausoleum in the study area, as shown in the following equation (Equation (4)). R = ∑ W i I ij = W 1 I 1j + W 2 I 2j + W 3 I 3j + W 4 I 4j + W 5 I 5j + W 6 I 6j + W 7 I 7j + W 8 I 8j (4) In Equation (4), W 1 -W 8 correspond to the importance degrees of various influencing factors, as calculated by the AHP method, considering the location of mausoleums, and I 1 -I 8 are the contribution values of various geographical environmental influencing factors, as calculated by the information method for different classification ranges.
According to the distribution characteristics of the mausoleums in the study area, the geographical environmental influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) of mausoleum geomancy are graded as follows: (f) River system: a multiring buffer zone is established for the river systems, which are divided into: >3 km, 2-3 km, 1-2 km, 0-1 km. (g) Slope: for the mausoleum points, multiple buffer zones are established with the slope as the distance: <10 km, 10-15 km, 15-20 km, and >20 km. (h) Aspect: the slope direction of each mausoleum is extracted; then, the slope direction is divided into one of the following categories according to the slope: flat slope, shaded slopes, semi-shaded slopes, semi-sunny slopes, and sunny slopes. Figure 4 shows the classification results of the assessment factors of the feng shui of mausoleums. divided into four levels: 1-1.000064, 1.000064-1.0016, 1.0016-1.04, >1.04.
(b) Elevation: according to the natural breakpoint method, the elevation is divided into four levels: 0-50 m, 50-100 m, and >100 m.
(c) Surface cutting depth: according to the natural breakpoint method, the surface cutting depth is divided into four levels: 0-7.3 m, 7.3-19.8 m, 19.8-39.5 m, and >39.5 m.
(e) Outfall: a multiring buffer zone is established for the outfall, which is divided into the following categories: >6 km, 3-6 km, and 0-3 km.
(h) Aspect: the slope direction of each mausoleum is extracted; then, the slope direction is divided into one of the following categories according to the slope: flat slope, shaded slopes, semi-shaded slopes, semi-sunny slopes, and sunny slopes. (a)

Calculation of the AHP Method
After establishing the assessment index system using the extracted eight factors influencing mausoleum geomancy, a judgment matrix of the mausoleum-related geomantic assessment factors in the study area was constructed according to the actual situation of study area combined with the experience of experts. This matrix is shown in Table 4. After normalizing the judgment matrix, the weight of each assessment factor was obtained, as shown in Table 5.

Calculation of the AHP Method
After establishing the assessment index system using the extracted eight factors influencing mausoleum geomancy, a judgment matrix of the mausoleum-related geomantic assessment factors in the study area was constructed according to the actual situation of study area combined with the experience of experts. This matrix is shown in Table 4. After normalizing the judgment matrix, the weight of each assessment factor was obtained, as shown in Table 5.  To prevent logical contradictions in the scoring matrix, it was necessary to test the consistency of the matrix and the coordination between the assessment factor weights. The test calculation is shown in Table 6, and CR = 0.096 < 1, indicating that the consistency test passed.

Calculation of the Information Content Method
According to the classification of the eight influencing factors of mausoleum geomancy combined with the current discovery of mausoleum distribution points, the information content method was used to calculate the amount of information contained within each influencing factor classification result. The results are shown in Table 7.  Table 7, T i represents the number of mausoleums within a given hierarchical influencing factor in the study area, T i /T represents the number of mausoleums within a given hierarchical influence factor in the study area as a proportion of the total number of mausoleums, and S i /S represents the proportion of the number of grading influencing factor units in the study area to the total number of units. Rank represents the order of the amount of information corresponding to each influencing factor from small to large.

Calculation of the Comprehensive Weighting Method (AHP-I)
The comprehensive weighting method (AHP-I) refers to a combination of the weights of the eight analyzed influencing factors, as calculated using the AHP method, with the information contents that were graded by each influencing factor, as calculated using the information content method (I). In this paper, the weighting calculations were combined with the relevant calculations shown in Tables 4 and 6. The calculation equation is shown as follows (Equation (5)). The analysis results are shown in Figure 5.  Table 7, Ti represents the number of mausoleums within a given hierarchical influencing factor in the study area, Ti/T represents the number of mausoleums within a given hierarchical influence factor in the study area as a proportion of the total number of mausoleums, and Si/S represents the proportion of the number of grading influencing factor units in the study area to the total number of units. Rank represents the order of the amount of information corresponding to each influencing factor from small to large.

Calculation of the Comprehensive Weighting Method (AHP-I)
The comprehensive weighting method (AHP-I) refers to a combination of the weights of the eight analyzed influencing factors, as calculated using the AHP method, with the information contents that were graded by each influencing factor, as calculated using the information content method (I). In this paper, the weighting calculations were combined with the relevant calculations shown in Tables 4 and 6. The calculation equation is shown as follows (Equation 5). The analysis results are shown in Figure 5. According to the geomantic mausoleums assessment results obtained with the comprehensive weighting method (AHP-I), the level of mausoleum geomancy in the main urban area of Nanjing was divided into five levels (relatively low suitability, low suitability, suitability, relatively high suitability, and high suitability) using the natural breakpoint method, which is commonly used in mathematical statistics. By combining the results with GIS-derived spatial statistics, the area of each level was counted. The results According to the geomantic mausoleums assessment results obtained with the comprehensive weighting method (AHP-I), the level of mausoleum geomancy in the main urban area of Nanjing was divided into five levels (relatively low suitability, low suitability, suitability, relatively high suitability, and high suitability) using the natural breakpoint method, which is commonly used in mathematical statistics. By combining the results with GIS-derived spatial statistics, the area of each level was counted. The results showed that the relatively low-suitability area accounted for 14.55% of the overall study area, the low-suitability area accounted for 25.40%, the suitable area accounted for 29.13%, the relatively high-suitability area accounted for 22.00%, and the high-suitability area accounted for 8.92%. These results are shown in Figure 6.
showed that the relatively low-suitability area accounted for 14.55% of the overall study area, the low-suitability area accounted for 25.40%, the suitable area accounted for 29.13%, the relatively high-suitability area accounted for 22.00%, and the high-suitability area accounted for 8.92%. These results are shown in Figure 6. Following Table 4 and in conjunction with the list of weighted information contents derived for mausoleums in the main urban area of Nanjing (Table 7), the characteristics of the relatively high-suitability and high-suitability regions within Nanjing city were determined using the following steps: (a) The surface roughness reflects the weakening effect of the surface on the wind speed and the influence of the land surface on wind-sand activity. Surface roughness is a corresponding influencing factor of "Long". The surface roughness ranges from 1.0016 to 1.22 in the feng-shui-suitable mausoleum area.
(b) The elevation reflects the terrain in a small region. Elevation is a corresponding influencing factor of "Long". Clearly, the elevation ranges from 1 to 100 m in the fengshui-suitable mausoleum area.
(c) The surface cutting depth reflects the hill density and is a corresponding influencing factor of "Sha". The surface cutting depth ranges from 19.8 to 39.5 in the feng-shuisuitable mausoleum area.
(d) Surface peaks reflect the area higher than the surrounding points within the analysis scope of a specific field. Surface peaks are a corresponding influencing factor of "Sha". The surface peaks range from 0 to 2 km in the feng-shui-suitable mausoleum area.
(e) The outfall reflects the outlet of flowing water in a basin and corresponds to the lowest point in the entire basin. The outfall is a corresponding influencing factor of "Shui". The outfall ranges from 3 to 6 km in the feng-shui-suitable mausoleum area.
(f) The river system reflects the water network system composed of all rivers, lakes and other water bodies in the analyzed basin. The river system is a corresponding influencing factor of "Shui". The river system ranges from 0 to 2 km in the feng-shui-suitable mausoleum area.
(g) The slope reflects the steepness and gentleness degrees of the analyzed surface units and is a corresponding influencing factor of "Xue". The slope ranges from 10 to 15 km in the feng-shui-suitable mausoleum area.
(h) The aspect reflects the direction faced by mountains from high to low. The extracted slope directions of mausoleums can be used to distinguish between sunny and shaded slopes. Aspect is a corresponding influencing factor of "Xue". Feng-shui-suitable mausoleum areas mainly comprise sunny slopes and semi-sunny slopes. Figure 6. The area proportion of each suitability level to the total study area.

Analysis and discussion
Following Table 4 and in conjunction with the list of weighted information contents derived for mausoleums in the main urban area of Nanjing (Table 7), the characteristics of the relatively high-suitability and high-suitability regions within Nanjing city were determined using the following steps: (a) The surface roughness reflects the weakening effect of the surface on the wind speed and the influence of the land surface on wind-sand activity. Surface roughness is a corresponding influencing factor of "Long". The surface roughness ranges from 1.0016 to 1.22 in the feng-shui-suitable mausoleum area. (b) The elevation reflects the terrain in a small region. Elevation is a corresponding influencing factor of "Long". Clearly, the elevation ranges from 1 to 100 m in the feng-shui-suitable mausoleum area. (c) The surface cutting depth reflects the hill density and is a corresponding influencing factor of "Sha". The surface cutting depth ranges from 19.8 to 39.5 in the feng-shuisuitable mausoleum area. (d) Surface peaks reflect the area higher than the surrounding points within the analysis scope of a specific field. Surface peaks are a corresponding influencing factor of "Sha". The surface peaks range from 0 to 2 km in the feng-shui-suitable mausoleum area. (e) The outfall reflects the outlet of flowing water in a basin and corresponds to the lowest point in the entire basin. The outfall is a corresponding influencing factor of "Shui". The outfall ranges from 3 to 6 km in the feng-shui-suitable mausoleum area. (f) The river system reflects the water network system composed of all rivers, lakes and other water bodies in the analyzed basin. The river system is a corresponding influencing factor of "Shui". The river system ranges from 0 to 2 km in the feng-shuisuitable mausoleum area. (g) The slope reflects the steepness and gentleness degrees of the analyzed surface units and is a corresponding influencing factor of "Xue". The slope ranges from 10 to 15 km in the feng-shui-suitable mausoleum area. (h) The aspect reflects the direction faced by mountains from high to low. The extracted slope directions of mausoleums can be used to distinguish between sunny and shaded slopes. Aspect is a corresponding influencing factor of "Xue". Feng-shui-suitable mausoleum areas mainly comprise sunny slopes and semi-sunny slopes.

Sensitivity Analysis for the AHP-Weighted Information Content Method
This paper analyzed the sensitivity of the AHP-weighted information method. In a sensitivity analysis, the values (336 in total) of parameters (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) are appropriately increased or decreased in the model simulation calculations [72,73]. By comparing the simulation results, the parameters with high sensitivities among the eight influencing factors can be identified to determine the influence of these parameters on the suitability evaluation simulation results.
In this paper, the sensitivity of the parameters was calculated individually and consecutively using the differential method; that is, one influencing factor was changed each time, and the influence of this change in one influencing factor on the suitability evaluation results (S i sensitivity) was tested. The calculation equations are shown in Equations (6) and (7).
where W i is the correlation weight of influencing factors, I ij is the value of each influencing factor, and C represents the other unchanged influencing factors when one influencing factor is changed; here, C is regarded as a constant. R * is the output value after the parameter change, R represents the output value before the parameter change, and ∆ i represents the variation amplitude of parameter I. Forty-two different assignments were made for each influencing factor, and 336 operations were carried out. Finally, the sensitivity of each influencing factor (S i ) was calculated. The larger the average values were, the more sensitive the influencing factors were to the suitability evaluation results. The calculation results are shown in Figure 7. According to the calculation results, the sensitivity value of elevation was the highest, at 0.044. The order of the other influencing factors, from large to small, was the surface cutting depth, aspect, slope, outfall, surface peaks, surface roughness and river system.

Random Mausoleums Used to Validate the Model Assessment and Analysis Results
According to the mausoleum geomantic assessment results obtained by the AHP-I method, 26 random mausoleum sites were introduced as verification data in this paper, and the proportion of random mausoleum sites in each suitability level area was calculated. The spatial distribution of the verification data is shown in Figure 8. The 26 random mausoleum sites included 19 general minister tombs (such as the Li Ruiqing tombs, Chen Zuoqi tombs, and Tao Xingzhi tombs), 3 family cemeteries (such as the Zhao family cemetery and the Zhang Jinghou family cemetery) and 5 martyr cemeteries. The sensitivity analysis of the AHP-I method shows that the sensitivity values of the eight influencing factors (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect) were all within the range of 0-0.5, indicating non-sensitive parameters. This shows that in the feng shui suitability evaluation of the mausoleum model, the eight selected influencing factors were relatively appropriate, and the influence of the calculated spatial eigenvalues of the influencing factors on the feng shui suitability evaluation of the mausoleum results was within an acceptable range. Therefore, the use of the AHP-I method to evaluate the feng shui suitability of mausoleums not only avoids the subjectivity of traditional feng shui judgment methods but also scientifically explains the four feng shui elements from the perspective of mathematical statistics, more closely linking feng shui to other disciplines such as geography. In the future, we may combine feng shui theory with other disciplines and technologies, such as machine learning methods [74,75].

Random Mausoleums Used to Validate the Model Assessment and Analysis Results
According to the mausoleum geomantic assessment results obtained by the AHP-I method, 26 random mausoleum sites were introduced as verification data in this paper, and the proportion of random mausoleum sites in each suitability level area was calculated. The spatial distribution of the verification data is shown in Figure 8. The 26 random mausoleum sites included 19 general minister tombs (such as the Li Ruiqing tombs, Chen Zuoqi tombs, and Tao Xingzhi tombs), 3 family cemeteries (such as the Zhao family cemetery and the Zhang Jinghou family cemetery) and 5 martyr cemeteries.

Random Mausoleums Used to Validate the Model Assessment and Analysis Results
According to the mausoleum geomantic assessment results obtained by the AHP-I method, 26 random mausoleum sites were introduced as verification data in this paper, and the proportion of random mausoleum sites in each suitability level area was calculated. The spatial distribution of the verification data is shown in Figure 8. The 26 random mausoleum sites included 19 general minister tombs (such as the Li Ruiqing tombs, Chen Zuoqi tombs, and Tao Xingzhi tombs), 3 family cemeteries (such as the Zhao family cemetery and the Zhang Jinghou family cemetery) and 5 martyr cemeteries. By comparing the mausoleum geomantic assessment results with the distribution of random mausoleum sites, the mausoleum geomantic assessment results were found to correspond well with the mausoleum distribution. As shown in Figure 9, of the 26 random tomb sites, 3 general minister tombs were located in relatively low-suitability and low- By comparing the mausoleum geomantic assessment results with the distribution of random mausoleum sites, the mausoleum geomantic assessment results were found to correspond well with the mausoleum distribution. As shown in Figure 9, of the 26 random tomb sites, 3 general minister tombs were located in relatively low-suitability and low-suitability areas, accounting for 11.54% of the overall number of sites. In contrast, five tombs were located in suitable areas, eight tombs were located in relatively high-suitability areas, and ten tombs were located in high-suitability areas, accounting for 19.23%, 30.77%, and 38.46%, respectively. Both the family cemeteries and the martyr cemeteries were located in high-suitability mausoleum geomancy areas.
suitability areas, accounting for 11.54% of the overall number of sites. In contrast, five tombs were located in suitable areas, eight tombs were located in relatively high-suitability areas, and ten tombs were located in high-suitability areas, accounting for 19.23%, 30.77%, and 38.46%, respectively. Both the family cemeteries and the martyr cemeteries were located in high-suitability mausoleum geomancy areas.

Analysis of Ming Xiaoling Mausoleum Geomancy
The Ming Xiaoling Mausoleum is the mausoleum of the emperor (Zhu Yuanzhang) and queen (Ma) of the Ming Dynasty. As the head of the Ming and Qing emperors' mausoleums in China, the Ming Xiaoling Mausoleum unified humanities and nature harmoniously and was listed as a world cultural heritage site in 2003 [76,77].
The Ming Xiaoling Mausoleum is located at Dulongfu (a large circular mound with a diameter of 325 to 400 meters and a height of approximately 70 meters) in the southern foothills of Zijin Mountain, an eastern suburb of Nanjing. Zijin Mountain has three peaks: the east, middle and west peaks. In geomancy, it is called "Huagai Three Peaks". Among the three peaks, the middle peak is the leading peak. The Ming Xiaoling Mausoleum is located within the foothills to the south of the middle peak. Looking north from Dulongfu, at the site where the Ming Xiaoling Mausoleum is located, the central mountain range of Zijin Mountain is curved; the arc mouth faces south, and the height of the mountain slowly drops towards Dulongfu. Dulongfu has an elevation above 118 m (in Table 7, this site belongs to rank 7), a surface roughness above 1.019 (in Table 7, it belongs to rank 1), a high terrain, and a winding mountain range. The surface cutting depth in the Ming Xiaoling Mausoleum is 17.78 (in Table 7, it belongs to rank 10), and the surrounding surface peak buffer area is 0-2 km (in Table 7, it belongs to rank 27). The superior "Long" and "Sha" of this area are important conditions that allowed it to become a geomantic treasure.
Water flows out of the mountains in the northeast direction of Dulongfu, divided into three royal rivers, respectively, from the front of the Minglou in Fangcheng, the Wenwu side door, and under the Royal River Bridge on Shixiang road; all of these rivers flow into a former lake in the southwest direction of Dulongfu. The surrounding outfall buffer area of the Ming Xiaoling Mausoleum is 0-3 km, belonging to the first classification level (in Table 7, this site belongs to rank 11), and its surrounding river system buffer area is 1-2 km (in Table 7, it belongs to rank 14). In eastern Dulongfu, Zijin Mountain extends southwards towards several hills. The slope type of Ming Xiaoling Mausoleum is between 180°-

Analysis of Ming Xiaoling Mausoleum Geomancy
The Ming Xiaoling Mausoleum is the mausoleum of the emperor (Zhu Yuanzhang) and queen (Ma) of the Ming Dynasty. As the head of the Ming and Qing emperors' mausoleums in China, the Ming Xiaoling Mausoleum unified humanities and nature harmoniously and was listed as a world cultural heritage site in 2003 [76,77].
The Ming Xiaoling Mausoleum is located at Dulongfu (a large circular mound with a diameter of 325 to 400 meters and a height of approximately 70 meters) in the southern foothills of Zijin Mountain, an eastern suburb of Nanjing. Zijin Mountain has three peaks: the east, middle and west peaks. In geomancy, it is called "Huagai Three Peaks". Among the three peaks, the middle peak is the leading peak. The Ming Xiaoling Mausoleum is located within the foothills to the south of the middle peak. Looking north from Dulongfu, at the site where the Ming Xiaoling Mausoleum is located, the central mountain range of Zijin Mountain is curved; the arc mouth faces south, and the height of the mountain slowly drops towards Dulongfu. Dulongfu has an elevation above 118 m (in Table 7, this site belongs to rank 7), a surface roughness above 1.019 (in Table 7, it belongs to rank 1), a high terrain, and a winding mountain range. The surface cutting depth in the Ming Xiaoling Mausoleum is 17.78 (in Table 7, it belongs to rank 10), and the surrounding surface peak buffer area is 0-2 km (in Table 7, it belongs to rank 27). The superior "Long" and "Sha" of this area are important conditions that allowed it to become a geomantic treasure.
Water flows out of the mountains in the northeast direction of Dulongfu, divided into three royal rivers, respectively, from the front of the Minglou in Fangcheng, the Wenwu side door, and under the Royal River Bridge on Shixiang road; all of these rivers flow into a former lake in the southwest direction of Dulongfu. The surrounding outfall buffer area of the Ming Xiaoling Mausoleum is 0-3 km, belonging to the first classification level (in Table 7, this site belongs to rank 11), and its surrounding river system buffer area is 1-2 km (in Table 7, it belongs to rank 14). In eastern Dulongfu, Zijin Mountain extends southwards towards several hills. The slope type of Ming Xiaoling Mausoleum is between 180 • -270 • (in Table 7, it belongs to rank 6), and the aspect is a sunny slope (in Table 7, it belongs to rank 20). The Ming Xiaoling Mausoleum chose Dulongfu as the "Xue" and is located in an area hidden from the wind and enjoys the presence of water. The Ming Xiaoling Mausoleum has four geomantic elements ("Long", "Sha", "Shui", and "Xue") and can be called a veritable geomantic treasure, as shown in Figure 10. Table 7, it belongs to rank 6), and the aspect is a sunny slope (in Table 7, it belongs to rank 20). The Ming Xiaoling Mausoleum chose Dulongfu as the "Xue" and is located in an area hidden from the wind and enjoys the presence of water. The Ming Xiaoling Mausoleum has four geomantic elements ("Long", "Sha", "Shui", and "Xue") and can be called a veritable geomantic treasure, as shown in Figure 10.

Conclusions
Based on GIS applications, this paper proposed a method that combines the correlation weighting method of each influencing factor with the classification method of individual influencing factors to assess mausoleum geomancy, which is different from the methodology of Chinese feng shui. First, the four geomantic elements of "Long, Sha, Shui, and Xue" were transformed into corresponding environmental influencing factors of mausoleums (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect). Then, the analytical hierarchy process (AHP) was used within the correlation weighting method of each influencing factor to calculate the weights of the eight influencing factors, and the information content method (I) was applied in the classification method of individual influencing factors to calculate the grading numerical index of each influencing factor. Finally, a mausoleum geomantic assessment model (AHP-I) was constructed by the comprehensive weighting method, and the mausoleum geomancy in the main urban area of Nanjing was assessed. Then, combined with the sensitivity analysis of the model, a random mausoleum site assessment analysis and a Ming Xiaoling Mausoleum geomantic assessment analysis, two methods were used to verify the assessment results of the model. The following conclusions were drawn: (1) According to the results of the mausoleum geomantic assessment conducted in the main urban area of Nanjing, the relatively low-suitability area accounted for 14.55% of the overall studied mausoleums, the low-suitability area accounted for 25.40%, the suitable area accounted for 29.13%, the relatively high-suitability area accounted for 22.00%, and the high-suitability area accounted for 8.92%.
(2) The results of the mausoleum geomantic assessment were consistent with the actual distribution of mausoleums, and based on GIS and AHP applications, the weighted information method is applicable for evaluating the feng shui suitability of mausoleums.

Conclusions
Based on GIS applications, this paper proposed a method that combines the correlation weighting method of each influencing factor with the classification method of individual influencing factors to assess mausoleum geomancy, which is different from the methodology of Chinese feng shui. First, the four geomantic elements of "Long, Sha, Shui, and Xue" were transformed into corresponding environmental influencing factors of mausoleums (surface roughness, elevation, surface cutting depth, surface peaks, outfall, river system, slope, and aspect). Then, the analytical hierarchy process (AHP) was used within the correlation weighting method of each influencing factor to calculate the weights of the eight influencing factors, and the information content method (I) was applied in the classification method of individual influencing factors to calculate the grading numerical index of each influencing factor. Finally, a mausoleum geomantic assessment model (AHP-I) was constructed by the comprehensive weighting method, and the mausoleum geomancy in the main urban area of Nanjing was assessed. Then, combined with the sensitivity analysis of the model, a random mausoleum site assessment analysis and a Ming Xiaoling Mausoleum geomantic assessment analysis, two methods were used to verify the assessment results of the model. The following conclusions were drawn: (1) According to the results of the mausoleum geomantic assessment conducted in the main urban area of Nanjing, the relatively low-suitability area accounted for 14.55% of the overall studied mausoleums, the low-suitability area accounted for 25.40%, the suitable area accounted for 29.13%, the relatively high-suitability area accounted for 22.00%, and the high-suitability area accounted for 8.92%. (2) The results of the mausoleum geomantic assessment were consistent with the actual distribution of mausoleums, and based on GIS and AHP applications, the weighted information method is applicable for evaluating the feng shui suitability of mausoleums.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy issues.