Geospatial modeling of the tropical cyclone risk in the Guangdong Province, China

Abstract Tropical cyclones (TCs) are prominent natural hazards in the Guangdong Province, China; hence, comprehensive risk assessment is vital to reduce potential losses due to TC hazards. A multi-indicator system was developed by incorporating 19 criteria under 3 risk components and a combination of the analytic hierarchy process and index of entropy was adopted to generate a TC risk map. Spatial statistics were applied to determine the spatial association and significant TC risk hotspots. The TC risk map showed that ∼40% of the study area is classified as a high or extremely high TC risk. Areas with the highest risk are mainly concentrated in the eastern region (ER), western region (WR), and Zhuhai and Shenzhen in the Pearl River Delta (PRD), whereas other areas in the PRD and mountainous region (MR) in the inland are classified as zones with the lowest risk. The results of the spatial autocorrelation analysis showed a significant positive spatial clustering in the spatial pattern of the TC risk, suggesting a varying spatial distribution and clustering pattern of the TC risk level that can be used as preliminary planning framework for tropical cyclone risk management. The results of this study provide theoretical reference for the design of disaster mitigation measures.


Introduction
Tropical cyclones (TCs) are strong cyclonic disturbances within synoptic systems that occur in tropical or subtropical oceans in association with strong winds, torrential rain, and storm surges, making them one of the most lethal and destructive natural disasters worldwide (Hoque et al. 2016;Hoque et al. 2018;Mansour 2019). The coastal region of China is one of the major cyclone-prone areas with a highly concentrated population and economy (Yin et al. 2013;Xu et al. 2015;Zhang et al. 2017;Zhong et al. 2018). Based on historical records from the Japan Meteorological Agency (JMA) from to 1951-2018, $1782 TCs were generated in the western North Pacific (WNP) Basin, 30% of which have made landfall in the southeast coastal region of China over the past 70 years, resulting in huge economic losses, casualties, house collapses, and crop damage (Chou et al. 2020;Wang et al. 2016;Cai et al. 2020). The results of numerous studies showed that the frequency of very intense TCs will likely increase in the future due to climate change (Chan and Liu 2004;Webster et al. 2005;Emanuel et al. 2008;Bacmeister et al. 2018;Sun et al. 2017). Hence, effective risk assessment of TCs is required to mitigate the effects of TCs on the local population, economy, and environment in coastal areas.
The disaster risk is complex and varies in different regions because it is driven by natural and anthropogenic disturbances. Hence, accurate information is crucial for the selection of reliable indicators for TC risk assessment (Azaz 2010;Hoque et al. 2017b). From the perspective of disaster system theory, the risk is generally associated with the interaction between the potential of hazard-induced factors and the vulnerability of the environment, people, and infrastructure exposed to the risk (Rashid 2013;Dewan 2013;Hoque et al. 2017b;Zhang et al. 2017). In reality, the regional disaster reduction capacity is an essential indicator for risk assessment, which may offset the hazard and vulnerability of TCs or overlapping relationships between them Hoque et al. 2017b;Hoque et al. 2018;Mansour 2019;Sajjad and Chan 2019;Nguyen et al. 2019;Sajjad et al. 2020). Previous research (Gao et al. 2014;Xu et al. 2015;Zhang et al. 2017;Alam et al. 2020;Cai et al. 2020) on TC risk assessment primarily focused on the characteristics of TCs occurring in an area and the disaster-bearing body. However, the mitigation capacity has been integrated in few TC risk assessment studies (Hoque et al. 2017a;Sajjad and Chan 2019).
Currently, methods based on index systems are the most widely used in TC risk assessment (Alam et al. 2020;Hoque et al. 2017b;Hoque et al. 2018;Mansour. 2019;Sajjad and Chan 2019;Sajjad et al. 2020) . Multi-criteria evaluation as a basis for the decision-making process is associated with the identification, quantification, and weighting of criteria, which is vital for tropical cyclone risk management (Mansour 2019). In numerous studies, either subjective (e.g., analytic hierarchy process, AHP) and objective methods (e.g., index of entropy, IoE; principal component analysis, PCA) or a combination of both have been used to allocate the relevant weights to each risk modeling index Mazumdar and Paul 2018;Kabenge et al. 2017;Cai et al. 2020;Wang et al. 2020). Among these methods, AHP and IoE are two commonly used subjective and objective methods Du et al. 2020). The AHP is an effective tool that can be used in the spatial multi-criteria decision-making process and is highly recommended for TC risk assessment (Yin et al. 2013;Zhang et al. 2017;Alam et al. 2020). However, it is substantially affected by the subjective imprecision of decision makers (Chen et al. 2015, Roy andBlaschke 2015;Hoque et al. 2019;Mansour 2019). In contrast, the IoE method can be used to determine the index weight based on the objective variation of each indicator but lacks subjective opinions based on prior knowledge (Liu and Chen 2019).
The spatial statistics method using Moran's I and G statistics is an effective tool to analyze the spatial association and significance level of geospatial phenomena (Mazumdar and Paul 2018;Majumder et al. 2019;Hu et al. 2019;Jha and Gundimeda 2019). The occurrence of TCs is uncertain with risks depending on the interactions among physical, socio-economic, and environmental factors. In other words, the TC risk can be regarded as a regional variable with significant spatio-temporal differences. Hence, spatial statistics based on TC risk assessment can be used to reveal spatial agglomeration characteristics to understand the spatial dimensions of the TC risk. To the best of our knowledge, spatial statistics have been applied in few studies to assess the spatial cluster pattern of the TC risk. In two previous reports (Sajjad et al. 2020;Mazumdar and Paul 2018), a spatial statistics approach was applied to evaluate the spatial heterogeneities of the typhoon risk in mainland China and to identify hotspots in Odisha, India, that are vulnerable to TCs.
Guangdong is an important region along the southeastern coast of China, which is frequently affected by tropical cyclones. Recent research results showed that TCs landfall in the Guangdong Province first and exit the area the latest and the economic loss in this region caused by severe TCs from 1994 to 2016 exceeded 200 billion yuan Wang et al. 2020;Xu and Deng 2020). In addition, the frequency of TCs made landfall in Guangdong is the highest and the Pearl River Delta (PRD) has been identified as the hotspot with the highest risk of TCs in coastal mainland China (Liu and Chen 2019;Sajjad et al. 2020). Among recent catastrophic TCs that made landfall in the Guangdong Province, such as Hagupit (2008), Usagi (2013), Rammasun (2014), Mujigae (2015), Hato (2017) and Mangkhut (2018) are responsible for considerable direct economic losses and human causalities (Yu et al. 2013;Fu et al. 2020).
To develop a systematic risk assessment procedure and identify high-risk regions affected by TCs in Guangdong, POI data were used in this study to map the mitigation capacity and an integrated AHP-IoE method was applied to conduct TC risk assessment. Based on the use of remote sensing, meteorological, socio-economic, and various observation data, a multi-indicator system incorporating the mitigation capacity was established, which can be used to evaluate the TC risk level in the Guangdong Province. Furthermore, Moran's I index and Gi statistics were employed to identify the spatial pattern of the TC risk in Guangdong. In this study, the spatial distribution of the TC risk level, spatial pattern of the TC risk, and mitigation measures are analyzed.

Study area
The Guangdong Province is located in the southernmost part of mainland China between 20. 22 -25.52 N and 109.65 -117.32 E. The total terrain includes an area of 179800 km 2 and is mainly covered by arable land, forests, and coastal wetlands (Liu and Chen 2019). Guangdong is the most economically developed and populous region in southern China, with a total gross domestic product (GDP) of RMB 9.7 trillion and a permanent resident population of 113.46 million by the end of 2018 (Wang et al. 2021). Guangdong is characterized by tropical and subtropical monsoon climate. The occurrence of TCs is extremely common from May-August. In addition, Guangdong is located at the northern verge of the South China Sea with a long coastline, which makes it the main area prone to TCs. To better assess the TC risk in different areas, the study area was divided into four zones based on their location and economic development ( Figure 1): Pearl River Delta (PRD), Eastern Region (ER), Western Region (WR), and Mountainous Region (MR).

Data acquisition
Multi-source datasets (TC, remote sensing, meteorological, vector, and socio-economic data) from various domestic and international data platforms and institutions were used in this study. The TC data include past TC information such as the location, minimum central pressure, duration, moving distance, and intensity of each TC every 6 h. Remote sensing data include digital elevation model (DEM), normalized difference vegetation index (NDVI), land use and land cover. Meteorological data include the wind speed and total precipitation. Vector data include road networks, railway networks, water networks, coastlines and point of interest (POI) data such as the medical, public, and educational infrastructure and charitable organizations in the cities of Guangdong. Socio-economic data include the population density, GDP, and historical disaster loss. Disaster loss data include casualties, direct economic losses, wrecked vessels, and the aquaculture loss area from 2014 to 2018. Detailed information on all data sources is provided in Table 1.

Material and methods
The risk assessment process used in this study consisted of four phases. First, a multi-indicator system was established based on the hazard, vulnerability, and mitigation capacity of the risk components. Second, we applied the combined AHP-IoE weighting method to assess the TC risk. Third, we validate the resulting TC risk map by comparing the historical disaster loss statistics. Finally, Moran's I and G statistics were used to analyze the spatial cluster pattern of the TC risk and identify high-risk hotspots in the Guangdong Province. The flowchart of the TC risk assessment of this study is shown in Figure 2.

Index system establishment
Because there is no consensus on the selection criterion for the conditioning factor used for TC risk assessment, 19 relative evaluation indicators under three risk component categories (hazard, vulnerability, and mitigation capacity) were examined based  Zhang et al. 2021). These risk criteria layers were obtained via data processing using projection, interpolation, clipping, distance analysis, and overlay analysis in ArcGIS. All factor layers were uniformly transformed into the Krassovsky ellipsoid and Albers projection at a 1 km grid scale, while indicators under each risk category were classified using the natural break statistical method, as listed in Table 2.

Hazard
A hazard is an event and phenomenon (like TC) that may affect the life, property, socio-economic activities, or environment (UNISDR 2009;Banholzer et al. 2014). The  probability of TC occurrence is determined by the characteristics of TCs and meteorological factors. Hence, six indices of TC-induced factors were identified to capture the TC hazard pattern ( Figure 3): frequency, TC track length, duration, central pressure, wind speed, and total precipitation. The frequency, track length, and duration of TCs are important indicators for the detection of potential hazards of TCs within the spatio-temporal range. The higher frequency of TC, the higher the degree of hazard to the region. The longer length and duration time, the greater the adverse impact of TCs on local areas. We determined the layer of frequency, TC track length, and duration by calculating the number of TC, the total length of TC track, and the duration of the continuous life history of each TC that crossed the Guangdong province between 1951 and 2018. Central pressure, wind speed, and rainfall correspond to the intensity of TCs. The layer of TC minimum central pressure can be established via interpolation based on the value of the minimum central pressure of TC every six hour. Due to limited availability of wind speed and precipitation data for each TCs, the yearly average wind speed data and total 24-hour rainfall data between 2010 and 2018 from the China Meteorological Administration were used to determine the spatial layers of wind speed and total precipitation.

Vulnerability
Vulnerability is the susceptibility of an element to environmental change or the extent to which a community and environment may be affected by a particular hazard (Rashid 2013;UNISDR 2009). In total, nine socio-economic and environmental factors, including the slope angle, NDVI, drainage density, proximity to highways, proximity to railways, proximity to the coastline, population density, GDP, land use and land cover were used to visualize the vulnerability across the Guangdong Province ( Figure 4). The slope angle, NDVI, and river density are important measures of physical vulnerability: A lower slope angle, lower vegetation coverage, and high river density suggest a high vulnerability to TC hazards. Population density, GDP, and land use and cover correspond to socio-economic vulnerability. The layer of slope is calculated from the Digital Elevation Model (DEM) using slope analysis. The layers of population density, GDP, and NDVI were extracted from spatial grid datasets of the Chinese population, GDP, and vegetation index based on a 1 km spatial resolution. Land use and cover were extracted from multi-period land use and land cover remote sensing datasets in China (CNLUCC). Furthermore, proximity to roads, railways, and coastline are important indicators of the susceptibility of exposed elements: The closer to the roads, railways, and coastline, the higher the impact. To produce the layers of proximity to coastline, road, railways, and drainage density, proximity analysis was performed using ArcGIS software.

Mitigation capacity
The mitigation capacity represents the ability to cope with extreme hazard events and implement planned mitigation measures to minimize disaster losses (Nguyen et al. 2019;. As a new spatial data source, POI data are of advantage because of the large sample sizes, detailed information, wide coverage, high recognition accuracy, and easy access, which improve the accuracy of the detection of urban microscale hotspots (Hu et al. 2016;Yao et al. 2017). In this study, four indicators, that is, the proximity to hospitals, proximity to emergency shelters, social assistance level, and education level, were selected to depict the mitigation capacity of TCs ( Figure 5). Healthcare facilities, such as public hospitals, clinics, and health centers, provide medical services in disaster areas. Emergency shelters are also effective mitigation measures, which provide important emergency and life-saving shelters during TCs. The level of social assistance corresponds to the distribution of charities, which provide relief and subsidies for TC disaster sufferers. The education level is quantified by the distribution density of higher education institutions and represents the local awareness of coping with natural disasters, which is a useful indicator of the disaster reduction. The layers of proximity to hospital, emergency shelter, and social assistance and education levels are created by kernel density analysis using POI data of healthcare facilities, emergency shelters, charities, and high education institutions.

Combination weighting method
In light of the problems of single weighting methods and to obtain more reasonable weighting results, combination weighting methods are applied to determine the weight of each TC hazard factor. Because of its simplicity and frequent use in Geographic Information System (GIS), weighted linear combination (WLC) was used in this study by aggregating all weighted layers. Based on the principle of the minimum relative information entropy minF ¼ P n j¼1 W j ðlnW j À lnW 1j Þ þ P n j¼1 W j ðlnW j À lnW 2j Þ s:t P n j¼1 W j ¼ 1; W j > 0, j ¼ 1, 2, :::, n , ( the combination weight equation is as follows ): where W j is the combined weight of the jth criterion and W 1j and W 2j represent the AHP weighting and IoE weighting of the jth index, respectively. The WLC is conducted by multiplying the indicators by the relative combined weights of the TC risk (Ogato et al. 2020): where W i is the combined weight for the ith indicator and X i represents the value of the ith indicator. The weight of each influencing factor was calculated based on the combination weight method, as shown in Table 3.

Risk assessment model
Based on existing literature (Hoque et al. 2018;Hoque et al. 2019;Mansour 2019;Alam et al. 2020) and natural disaster risk theory, a well-accepted risk assessment model for TCs was adopted in this study. The composite risk assessment model considers the hazard-induced components; vulnerability of the hazard-bearing body, including natural and socio-economic components; and capacity of mitigation or resilience. The following equation was used to generate the TC risk map: The risk index was then standardized using a common scale of 0-1 and the following equation: where S is the standardized value, x is the value of a single cell in the spatial raster layer, and min and max represent the minimum and maximum values, respectively. Based on the standardized risk index, the risk was classified into five levels: very low, low, moderate, high, and very high. For more accurate and explicit visual assess the spatial characteristics of the TC risk, the study area was divided into 20 Â 20 km evaluation units with a total of 547 grids and the value of the grid center point was used as the TC risk index of each risk unit (Das et al. 2020). Using the spatial statistics technique of ArcGIS 10.4 software, the characteristics of the spatial distribution of the TC risk level at the grid scale were further analyzed.

K-Means cluster method
The K-means clustering method was applied to determine the damage level. Although several methods are currently available for cluster analysis, the K-means cluster method is widely used clustering algorithm due to its simple implementation and efficiency (Tang et al. 2020). K-means cluster is an unsupervised clustering method that can cluster the input samples with similar characteristics into one class. The main steps is as follows (Xu et al. 2018): 1. Randomly select the k cluster centers c i : 2. For each sample point, calculate the distance between each sample point and cluster center. 3. Classify the point nearest to the same central point as the representative cluster. 4. Update the center point for each cluster class. 5. Repeat the iteration until the best clustering is obtained. The goal of K-means algorithm is to minimize the sum of the squared error (SSE). The objective function J is determined as follows (Jain 2010): kx j À c i k 2 , 1 j n and 1 i k Where x j is the j th sample point, c i is the centroid of the i th cluster, k is the number of cluster centers, and n is the number of sample point. The SSE value is smaller and the clustering is optimal when the sample is divided into K clusters.

Spatial statistics method
Spatial autocorrelation analysis, including global and local spatial autocorrelation analysis, is a spatial statistical method that is mainly used to test whether a certain attribute value in a geospatial region significantly correlates with the attribute value of the same feature in adjacent spatial areas Jing et al. 2020). In this study, Moran's I index and Getis-Ord Gi Ã statistics were used to investigate the spatial autocorrelation of TC risk levels and identify high-risk hotspots in the study area.
The global Moran's I index was used to measure the spatial autocorrelation of an element in the entire area. The global Moran's I index is calculated as follows (Anselin 2010;Zou et al. 2018): where W ij ¼ 1, adjacent D p , otherwise ; N is the number of units; X i and X j are the attribute values of the units i and j, respectively; XX is the average value of the attribute; W ij denotes the spatial weight matrix, and the value of I ranges from À1 to 1.
The local Moran's I index was used in this study to describe spatial clusters and spatial outliers of TC risk levels within units, which can be displayed in the form of the local indicator of spatial association (LISA) graphics . The local Moran's I index is calculated as follows: where I i is the LISA index and S 2 is the sample variance. An I i value greater than zero indicates a positive spatial autocorrelation in the local spatial unit, which is represented by a "high-high" and "low-low" spatial aggregation pattern. In contrast, an I i value below zero indicates a negative spatial autocorrelation of the local spatial unit, presenting a "high-low" and "low-high" spatial aggregation pattern. Hotspot analysis was carried out using Getis-Ord Gi Ã statistics to identify the statistically significant hotspots and cold spots in TC-affected and unaffected regions. The z-score of Getis-Ord Gi Ã statistics is helpful for the determination of locations with spatial clusters of high or low values, with a confidence level of 99% (a < 0.01), 95% (a < 0.05), and 99% (a < 0.1) (Das et al. 2020;Hussain and Islam 2020). The Getis-Ord Gi Ã statistics can be calculated as follows (Ord and Getis 2010): where n is the number of features; X j and X i are the attribute values of features i and j (j 6 ¼ i), respectively; W ij is the spatial weight of features i and j; S is the standard deviation, and G i Ã is the calculated z-score. The higher the statistically significant positive z-score is, the closer is the high-value cluster and vice versa.

Hazard analysis
The hazard map (Figure 6a) shows that the areas with the highest hazard risk are Zhanjiang, Huizhou, and Heyuan, whereas the hazard risk in Qingyuan, Zhongshan, Zhuhai, Dongguang, Shanwei, Shantou, and Chaozhou is relatively low. The area with an extremely high hazard level accounts for 22.08% of the total area, mainly distributed in the southwest of the WR and northeast of the ER. Areas with very low and low hazard levels are distributed across the whole province, except for the WR, covering 19.66% and 16.45% of the total area, respectively. Maoming, Jiangmen, and Guangzhou were classified as high-risk areas and Yunfu, Yangjiang, Zaoqing, and Meizhou were classified as moderate-risk areas, covering 15.70% and 26.11% of the study area, respectively (Figure 7). The results show that areas with high hazard levels strongly correlate with the high frequency of TC landfalls or TCs passing through these areas and are prone to TC with high intensity. Most of the areas in the northeastern and southwestern parts of the Guangdong Province or close to the coast are characterized by a relatively moderate hazard level. Based on the topographic characteristics of Guangdong (Figure 1), the entire area is higher in the north and lower in the south. Furthermore, the mountains present a northeast-southwest alignment parallel to the coastline and these areas are thus more likely affected by heavy rainfall and strong winds.

Vulnerability analysis
The spatial distribution of the vulnerability of TCs in Guangdong decreases from coastal to inland regions, with high values in coastal areas and low values in the inland region (Figure 6b). Approximately 18.90% of the study area belongs to very highly vulnerable zones, which mainly are located along the coast. Areas with very low and low vulnerability levels cover 18.7% and 18.89% of the study area, respectively ( Figure 7). Areas with very high or high vulnerability levels include Zhangjiang, Maoming, Yangjiang, Jiangmen, Zhuhai, Zhongshan, Shenzhen, Southeast Huizhou, Shanwei, Jieyang, Shantou, and Chaozhou, whereas areas with very low or low vulnerability levels include Shaoguan, Qingyuan, Zaoqing, and parts of Heyuan, Meizhou, and Yunfu, accounting for 18.70% and 18.89% of the study area, respectively. The results also show that most areas in Yunfu, Heyuan, Meizhou, Foshan, and Dongguan and the central part of Guangzhou have medium vulnerability levels, accounting for 20.62% of the study area. These results indicate that most areas near the coast is more vulnerable to TC occurrences.

Mitigation capability analysis
The map of the mitigation capability (Figure 6c) shows that the spatial distribution of the mitigation capacity is centered in the PRD urban agglomeration and gradually weakens toward the surrounding cities. Therefore, the mitigation capacity in the PRD is higher, especially in Guangzhou, Shenzhen, Foshan, Dongguang, and Zhongshan, which were classified as zones with very high mitigation capacities. The periphery of these cities was classified as zone with a high mitigation capacity, accounting for 29.52% of the study area. Very low or low mitigation capacities are mainly distributed in Qingyuan, Maoming, and the periphery of Zhangjiang, Zaoqing, Shaoguan,Shanwei, and Jieyang, covering 5.80% and 14.92%, respectively (Figure 7).
The results indicate that most areas in the Guangdong Province have a moderate mitigation capacity and the ability to prevent and reduce disasters. Areas with a higher mitigation capacity are mainly those areas with a more developed economy, advanced medical facilities, and professional rescue agencies.

TC risk assessment
The composite risk map of TCs was produced by combining the hazard, vulnerability, and mitigation capacity maps. As shown in Figure 6d and Figure 8, areas near the coast, including Shantou, Shanwei, Jieyang, and Chaozhou in the ER as well as Zhanjiang, Maoming, and Yangjiang in the WR, exhibit high-or very high-risk levels, especially Maoming and Zhanjiang occupy the largest area proportion of the very high risk level, followed by Shanwei, Yangjiang, Jieyang, Jiangmen and Chaozhou. However, except for Zhuhai and Shenzhen, which exhibit very high-risk levels, most areas in the PRD exhibit the lowest risk level. In addition, the periphery of Yunfu and Yangjiang in the MR, Jiangmen, Zaoqing, and Huizhou in the PRD, Shanguan, Heyuan, and Meizhou in the MR exhibit a medium risk level and the centers of areas in Shaoguan, Meizhou, and Heyuan in the MR, Huizhou, Zaoqing, Zhongshan, and Jiangmen in the PRD, and Yunfu in the WR exhibit a low risk level. As shown in Figure 7, the distribution of different risk levels in the Guangdong Province is homogeneous; approximately 19.68% of the study area exhibits an extremely high-risk level and areas with a high, medium, or low risk levels account for 20.08% of the study area.

Validation of the TCs risk map
Currently, no absolute methods for verifying the generated TC risk map are available, hecne, absolute TC risk assessment is not possible. Therefore, the TC risk map was compared with historical disaster loss statistics derived from the Guangdong Marine Disaster Bulletin. Due to limited data availability and information integrity, historical disaster loss data from 2014 to 2018 were collected. In this study, the direct economic loss, TC-affected population, aquaculture loss area, and wrecked vessels were used as indicators of TC damage loss, while the K-means clustering method was applied to determine the level of TCs damage loss based on the cluster number. The TC damage loss level map (Figure 9) shows that the damage caused by TC hazards is more serious in the areas near the coast. These areas include Zhanjiang, Maoming, Shantou, Shanwei, and Jieyang in the ER and WR, which suffered relatively heavy losses due to TCs, whereas the damages in Guangzhou, Dongguan, and Foshan in the PRD as well as Shaoguan and Qingyuan in the MR were very small, which basically agrees with the results obtained from the TC risk map.

Spatial autocorrelation analysis of the integrated TCs risk
The results of the global spatial autocorrelation report show that the global Moran's I index of the TC risk is 0.71, the z-score is 47.5, and the p-value is below 0.01, which means that there is a 99.9% likelihood that the significant spatial clustering phenomenon observed in this study is not random. Based on the results of the local Moran's statistics (Figure 10), the spatial association of the TC risk presents high-high and low-low aggregation and effects of high-low aggregation and low-low aggregation are not notable, which indicates that there the risk units are spatially homogeneous. The results show that high-high aggregation areas are mainly located in both the ER and WR. These areas have high TC risk values surrounded by similar high-risk values. Low-low aggregation areas are mostly located in the PRD and MR, indicating that the risk values in these areas and adjacent units are low. The LISA significance map shows that the significance level of most of the high-high and low-low aggregation regions is 0.01. A significance level of not significant is mainly located in these spatial outlier areas.
Similar to the LISA cluster results, the results of the hotspots analysis ( Figure 11) show that the WR and ER exhibit a statistically significant spatial cluster of high TC risk hotspots (red color) with 99%, 95%, and 90% confidence levels, whereas regions not affected by TCs present statistically significant spatial clusters of low TC risk cold spots (blue color). In addition, in agreement with the spatial distribution of the TC risk, areas with high-risk spatial clustering are mainly located in Zhanjiang, Maoming in the WR, Shantou, Shanwei, Jieyang, and Chaozhou in the ER and areas with lowrisk spatial clustering are scattered in Guangzhou, Foshan, and Dongguan in the PRD, Shaoguan, Meizhou, and Heyuan in the MR, corresponding to high-risk (hotspots) and low-risk (cold spots) regions. Therefore, the results obtained from the Getis-Ord Gi Ã statistics confirm that the high-high and low-low aggregation areas of the TC risk derived from the local spatial autocorrelation analysis are statistically significant hotspots and cold spots, respectively.

Discussion
Risk assessment is an important basis for disaster prevention as well as post-disaster mitigation and risk reduction. Guangdong is one of the main TC-prone zones in China, but only limited studies have focused on the identification of the TC risk by integrating multi-criterion methods and spatial statistics. In this study, we evaluated the TC risk level and spatial pattern of the TC risk in Guangdong by investigating the systematic interaction between prominent physical and anthropological factors and the integrated mitigation capacity. Additionally, in this study, an approach combining AHP-IoE and spatial statistics were used for the first time to conduct TC risk assessment. The results were compared with historical disaster loss records. The results of the validation indicate that the integrated AHP-IoE method and spatial statistics are feasible for the evaluation of the TC risk, which can be used as references for future work on risk prevention and disaster reduction. For example, the risk map generated in this study shows that the TC risk level in Guangdong decreases from both ends of the coastal areas to the inland areas, which is consistent with the findings of previous studies (Li and Li 2013;Cai et al. 2020). Based on our results, the areas with the highest risk are mainly concentrated in the ER, WR, Zhuhai, and Shenzhen in the PRD near the coastline, whereas other areas of the PRD and MR in the inland were classified as zones with the lowest risk. These results agree with an earlier study in which assessed the risk level of TC disasters in the Guangdong Province. The authors showed that areas with the highest TC risk are located along the western, central and eastern coast in Guangdong (Zhang et al. 2011).
Another important finding is that the risk areas affected by TCs regionally vary, which depends on various factors such as the regional geographical location, terrain characteristics, and socio-economic development. We observed that several areas with a lower risk had a higher TC risk, especially in the PRD. For instance, Shenzhen and Zhuhai exhibited a high efficiency with respect to disaster mitigation and low TC hazard levels; however, these areas remained at high risk due to high population density and proximity to the coastline, which aggravates the vulnerability to TCs. Guangzhou, Foshan, Dongguang, and Zhongshan, were classified as very low risk areas, and are highly efficient in coping with TCs, despite having highly concentrated population. As the most frequent landfall area of TCs, Zhanjiang in the WR has become another very high-risk area due to the very high hazard level of exposure to severe TCs and very high vulnerability with the poor infrastructure to resist TCs. Shantou, Shanwei, Jieyang, and Chouzhou in the ER were classified as high-risk areas. This is mainly due to their geographical locations near the coastline, lack of mitigation measures, which increases their vulnerability to TCs, and lack of efficient coping measures. The centers of Yunfu, Shaoguan, Heyuan, and Meizhou in the MR were classified as areas with the least risk, corresponding to the low hazard level and vulnerability due to the high elevation and being far from the coastline; thus, these areas are less affected by TCs. This suggests that disaster prevention and mitigation measures should be implemented according to local conditions. Based on our results, areas with a high TC risk in the PRD (Zhuhai and Shenzhen) should continue to reduce the vulnerability of the disaster-bearing body to minimize the effects of TCs on local areas. For example, a more flexible evacuation plan should be constructed in areas with very high risk. In addition, such a plan should be followed in case of a TC occurrence to avoid secondary disasters caused by crowds. Alternatively, public schools, Children's palaces, museums, parks, and other public places can be used as open spaces for emergency disaster prevention when it becomes necessary to transfer highly densed populations, which can be integrated with early warning, emergency channel activation, and evacuation command systems (Mansour 2019). In areas near the coastline in the ER (Shantou, Shanwei, Jieyang, Chaozhou) and WR (Zhanjiang, Maoming), effective tropical cyclone precaution plans and mitigation measures, such as the strengthening of dams, delaying of economic activity, closing traffic roads, and construction of emergency shelters, should be implemented to reduce the TC risk (Haque and Blair 1992;Quader et al. 2017). More importantly, the TC mitigation capacity and awareness of disaster threats must be increased by regularly organizing lectures on TCs and conducting evacuation and refuge drills.
Our study has some limitations such as the irrationality of the factor index selection, lack of detailed historical meteorological data and disaster loss data, and different spatial resolution of available data from various data sources. These limitations can be addressed by improving the data acquisition technology. In addition, due to the difficulty with respect to the data acquisition, other indicators of the mitigation capacity criteria, such as the number of medical facilities, medical staff, and financial support, have not been considered. However, the findings of this study have implications for the risk reduction of TCs in the Guangdong Province. In summary, the results of this study provide theoretical basis for the design of disaster reduction plans and TC prevention and mitigation measures and thus are beneficial for TC risk management. The integrated AHP-IoE method and spatial statistics based on GIS techniques are useful for TC risk assessment and can be applied to other disaster research. Furthermore, the hazard level, vulnerability, and mitigation capacity maps constructed for the Guangdong Province can also be used as baseline information for local authorities to understand the risk status of different regions and formulate special risk reduction strategies based on information on disaster-inducing factors, disaster-breeding environments, disaster-bearing bodies, and the disaster reduction capacity in this region. A more reliable TC risk map is warranted and the impact and recovery of TCs should be assessed using additional data and advanced processing techniques.

Conclusion
In this study, an integrated TC risk map was constructed by combining the AHP and IoE methods to identify the degree of the TC risk in the Guangdong Province and analyze the spatial pattern of the TC risk at the 1km grid scale using the spatial statistics of the ArcGIS platform. Nineteen indicators representing hazard level, vulnerability, and mitigation capacity were aggregated using the WLC method to evaluate the hazards of TCs-induced factors, vulnerability of the disaster-bearing body, and mitigation capacity of TC disasters, respectively. Our results indicate that the TC risk level decreased from both ends of coastal areas to the inland areas, while $40% of the areas had an high or very high TC risk. Areas with the highest risk are mainly concentrated in the ER, WR, and Zhuhai and Shenzhen in the PRD near the coastline, whereas areas in the PRD and MR in the inland were classified as zones with the lowest risk. The spatial pattern of the TC risk presents a significant positive spatial clustering. Areas near the coast are subject to the highest TC risk, especially the coastal areas of Zhanjiang, Maoming, Shantou, Shanwei, Jieyang, and Chaozhou, which are significant hotspots with high TC risk levels. In contrast, Guangzhou, Dongguan, Foshan, Shaoguan, Heyuan, and Meizhou were classified as zones with the lowest risk (i.e., cold spots). The present TC risk assessment study suggests that models based on a GIS-integrated multi-index evaluation method and spatial statistics are applicable and reliable and our results can be used as reference for designing disaster mitigation and adaptation policies.

Disclosure statement
No potential conflict of interest was reported by the authors.

Funding
This study was supported by Talent introduction project of the Jinan University, and Innovative Research Team of the Department of Education of the Guangdong Province (Grant No.2020KCXTD005).

Data availability statement
The data presented in this study are openly available in: