Identification of Typhoon-Vulnerable Areas and Countermeasures in High-Density Coastal Cities: The Case of Macau

: Typhoons are extremely severe weather events which seriously threaten the safety of people’s lives and properties. Therefore, identifying and controlling typhoon disaster hazards have become important research topics. The spatial–temporal characteristics of typhoons are analysed using the typhoon disaster data in Macau from 2000 to 2020. Computational fluid dynamics (CFD) numerical simulation is adopted to understand the 3D urban wind environment. Moreover, the ‘exposure, sensitivity and adaptation’ evaluation model is applied to construct the study framework. To calculate urban disaster vulnerability, the Create Fishnet tool is used to divide the city of Macau into 470 grids. The principal component analysis method is used to reveal the factors that significantly affect the typhoon’s vulnerable areas. Result shows that 31.27% of grids are severely vulnerable. In addition, six principal components are identified, including indicators such as population density, building area ratio, mean elevation and wind speed. This study verifies the feasibility of wind speed data obtained by CFD in the typhoon evaluation model. Moreover, it provides a reliable reference guide for future urban microlevel studies.


Introduction
Extreme weather events threaten the sustainability of coastal cities. Typhoons are one of the most important extreme weather scenarios in coastal communities [1,2].Similar to hurricanes, typhoons are cyclones that primarily affect the western Pacific coast.They also considerably threaten personal safety and urban construction [2].Along with the global trend of rising urban populations, low-elevation coastal cities are at substantial risk when facing typhoon disasters currently [1].The cumulative economic damage caused by Cyclone Amphan in 2020 is as high as USD 13.6 billion [3].Furthermore, the major contributing areas of human casualties and economic losses are the vulnerable areas of coastal cities [4].
Therefore, urban disaster vulnerability (UDV) assessment is necessary in the context of typhoons.It is not only beneficial for urban resilience building but also positive for sustainability [2].However, the current UDV of coastal cities in cyclone contexts is not yet fully satisfactory.Three issues must be addressed.Firstly, many coastal cities have attracted a large number of immigrants due to economic and other factors.This phenomenon is not slowing down at present [5], but the research on coastal high-density cities is limited.Secondly, the spatial layout and structure of the city in the urban spatial pattern are necessary for urban development [6].In particular, the height and distribution of buildings have a significant impact on the distribution and direction of typhoon wind speed [7].However, few studies have explored the urban spatial morphology associated with resilience in the context of typhoons.Thirdly, a notable research gap exists in the academic community in terms of assessing the vulnerability of coastal cities to typhoons.For instance, an assessment framework and indicator system specifically for typhoon vulnerability remains unavailable.Current research focuses more on urban vulnerability.Meanwhile, assessment frameworks for typhoon vulnerability of coastal cities are neglected [6,[8][9][10][11].
In response to the first research question, we believe that a UDV assessment of coastal cities characterised by high population density must be conducted at this stage given the global trend of rising population.Another reason is the continuous population flow to coastal low-elevation cities over the previous decades [1].This condition can potentially maintain increasing population densities in coastal low-elevation cities in the future.For the second research question, computational fluid dynamics (CFD) has a huge advantage in modelling typhoon disasters in three dimensions [12,13].However, to the best of our knowledge, only a few studies have attempted to combine this approach with UDV.We believe that a typhoon simulation using CFD is required in conjunction with urban spatial morphology to assist in UDV assessment.For the third research question, the exposuresensitivity-adaptation model has been proposed as early as 2007 and has been accepted by many studies [14].However, studies that utilised sufficient data across the three dimensions to fully validate the model are limited.Therefore, an assessment framework and indicator system using this model must be developed specifically for typhoon vulnerability.
Macau, a former Portuguese colony, is located in southern China, close to Hong Kong and Guangzhou.It is one of the most densely populated cities in the world [15].In 2022, the population per square kilometre of Macau had exceeded 20,000, reaching 20,620 [16].Macau suffers from frequent typhoon disasters due to its location in the coastal area and its subtropical climate [17].Two of the worst typhoons in Macau's history, Hato and Mangkhut, caused severe economic losses and casualties [18].In addition, Macau has a well-developed economy.At its peak, it ranked first in the world in terms of gross domestic product (GDP) per capita at purchasing power parity [19].Therefore, Macau has a complex urban spatial pattern.The disaster simulation must be completed during a typhoon by using CFD.Finally, the Macau government has conducted sound statistical data disclosure, such as the age structure of the residents as fine as that of single residential buildings [16].This condition contributes to the exposure-sensitivity-adaptation modelling.The content of the appeal reveals that the case study of Macau can provide a key reference for the construction of indicator systems for coastal cities worldwide.

Literature Review
In 2007, the Intergovernmental Panel on Climate Change (IPCC) report indicated that vulnerability assessment is a vital tool in addressing the challenges of climate change.The report defined vulnerability as an integrated concept composed of three dimensions: exposure, sensitivity and adaptive capacity [20], which have been widely applied [21,22].At the urban level, especially when facing the threat of extreme weather events, conducting a comprehensive vulnerability assessment is particularly crucial.Although academic research on vulnerability assessment has progressed, many deficiencies remain and must be improved.
Firstly, a significant gap exists in current academic research regarding the vulnerability assessment and study of high-density urban areas under the impact of typhoons.Deng et al. (2020) suggested that these regions face unique challenges due to high population density, concentrated infrastructure and limited resources [23].However, these characteristics have yet to receive adequate attention in the existing literature.Although many high-density cities are located in coastal areas, not all coastal cities exhibit the characteristics of highdensity urban areas.Therefore, conducting specialised vulnerability assessments for these cities would provide new perspectives and a deeper understanding of the field.Ku et al. (2021) assessed the vulnerability of coastal areas, but their assessment lacked explicit focus on the unique requirements and challenges of high-density urban regions [24].In the field of urban planning, methods for typhoon vulnerability assessment have been widely applied but have yet to comprehensively cover all types of cities, mainly high-density urban areas.The risks faced by cities during typhoons stem not only from the destructive power of the typhoons themselves but also from multiple factors, such as changes in surface cover caused by urbanisation processes, inadequate drainage systems and the lack of flood defence facilities [25].In-depth research on typhoon vulnerability in high-density urban areas becomes particularly critical to more comprehensively and accurately assess their vulnerability and susceptibility during typhoons.
Secondly, in the research on typhoon vulnerability, scholars generally focus on the impact of factors, such as climate change, population density and economic conditions [26], whereas the analysis on urban spatial forms must still be improved.Although He et al. (2017) investigated the coupling mechanism between coastal urban spatial form and mesoscale wind environment, emphasising the critical role of urban form in regulating wind environment, subsequent studies, such as that of Gao et al. (2020), did not delve into the specific impact of urban form on disaster risk when assessing the typhoon disaster risk in Zhuhai City [27,28].Su et al. (2022) highlighted the response and recovery capabilities of urban ecosystems in typhoon disasters and the key role of urban form in reducing vulnerability.However, they failed to investigate how urban spatial layout can reduce typhoon damage, especially in terms of 3D scale analysis [29].Current typhoon vulnerability research is often limited to considerations of 2D spatial factors, with a relative lack of in-depth analysis of urban 3D spatial form [30].The spatial layout and structure of cities, especially the height and distribution of buildings, significantly impact the distribution and direction of wind speed.In addition, the distribution of urban green spaces regulates the urban climate and effectively slows down wind speed during typhoons, providing additional protection for the city [31].Therefore, future research should extend to the analysis of urban 3D spatial form to more comprehensively understand its effects on typhoon vulnerability and provide more precise guidance for urban planning and disaster risk management.
Finally, in assessing the vulnerability of coastal cities to typhoons, the academic community currently faces a notable gap in research, namely, the absence of a dedicated assessment framework and indicator system for typhoon vulnerability.Jeong and Cheong (2012) and Gao et al. (2020) laid some groundwork in this area; they primarily subsume typhoon vulnerability within a broader analysis of urban-coupled vulnerability, failing to fully account for the urban spatial vulnerability factors specific to typhoon disasters [28,32].Kim et al. (2020) attempted to develop a typhoon vulnerability function using loss records from Typhoon Maemi.However, their research should have comprehensively considered the impact of urban spatial layout and environmental characteristics [33].Similarly, Yan et al. (2023), in exploring the spatiotemporal variations in typhoon risk in Guangdong Province, needed to provide a specialised assessment indicator and framework for typhoon vulnerability [34].Existing research tends to incorporate typhoon vulnerability into a more general framework of urban coupled vulnerability [35][36][37], which often overlooks the urban spatial vulnerability factors unique to typhoon disasters.However, the impact of typhoons on cities is deeply influenced by their spatial layout and environmental characteristics [27], and thus, these specific spatial factors must be considered when assessing typhoon vulnerability.This urgently calls for the academic community to develop an assessment framework and indicator system that starts from the perspective of urban space and is aimed explicitly at typhoon vulnerability to more accurately identify and quantify the vulnerability of cities in the face of typhoon threats.To fill this research gap, future studies should focus on constructing an integrated assessment framework that includes specialised indicators for typhoon vulnerability whilst considering the 3D spatial characteristics of cities, such as the height, density and building layout, as well as the distribution of green spaces and water bodies.
In summary, despite substantial research on urban vulnerability, some academic gaps remain.The main gaps are as follows: (1) the absence of typhoon vulnerability research in high-density urban areas and (2) the insufficiency of typhoon vulnerability research, in terms of analysis and discussion, within the urban 3D space.
The present study aims to fill these gaps.Firstly, given that existing research ignored the impact of vulnerability assessment in high-density urban areas, we focus on identifying typhoon disasters and exploring adaptation strategies in high-density regions.This study conducts a spatiotemporal quantitative analysis of all typhoons that affected Macau from 2000 to 2020, investigating the spatial distribution and seasonal characteristics of typhoon disasters and revealing the main tracks of these typhoons.Moreover, line density analysis and seasonal variations are performed through GIS spatial statistics.Secondly, since 3D data can offer more details and precision than 2D data, the typhoon wind environment in Macau is numerically simulated at a 3D scale, considering the main characteristics and features of Macau's topographical surface.Thirdly, we developed an urban disaster vulnerability (UDV) assessment model using the 'exposure-sensitivity-adaptation' framework.This model aims to construct an urban disaster vulnerability assessment framework specifically for typhoon vulnerability, and a corresponding index system is established to systematically assess and quantify the vulnerability to typhoon disasters in cities.The model assesses the entire city of Macau on a grid scale to identify areas vulnerable to typhoon disasters.In addition, principal component analysis (PCA) reveals the significant factors influencing the areas prone to typhoons in Macau.

Study Area
The Macau Special Administrative Region (SAR) of the People's Republic of China, commonly referred to as Macau, is a distinctive geographical entity situated on China's southeast coast.It is not just a pinpoint on a map, at a longitude of 113.5 • east and a latitude of 22.2 • north, but a region with a rich spatial diversity.Macau holds a strategic position, bordered by Zhuhai City in Guangdong Province to the north and facing Hong Kong across the Pearl River Estuary to the east (Figure 1).ronment in Macau is numerically simulated at a 3D scale, considering the m istics and features of Macau's topographical surface.Thirdly, we developed aster vulnerability (UDV) assessment model using the 'exposure-sensitivity framework.This model aims to construct an urban disaster vulnerabilit framework specifically for typhoon vulnerability, and a corresponding ind established to systematically assess and quantify the vulnerability to typhoo cities.The model assesses the entire city of Macau on a grid scale to identify able to typhoon disasters.In addition, principal component analysis (PCA significant factors influencing the areas prone to typhoons in Macau.

Study Area
The Macau Special Administrative Region (SAR) of the People's Repu commonly referred to as Macau, is a distinctive geographical entity situate southeast coast.It is not just a pinpoint on a map, at a longitude of 113.5° e tude of 22.2° north, but a region with a rich spatial diversity.Macau holds a sition, bordered by Zhuhai City in Guangdong Province to the north and Kong across the Pearl River Estuary to the east (Figure 1).
The territory of Macau is a vibrant tapestry of interconnected spaces, e the Macau Peninsula, Taipa Island and Coloane Island.These areas are link and causeways, creating a cohesive urban expanse that transcends individ form a unified urban landscape.This spatial configuration is visually repres ure 2, which illustrates the urban continuum and maritime boundaries, prov prehensive view of Macau's geographical layout and spatial characteristics.The territory of Macau is a vibrant tapestry of interconnected spaces, encompassing the Macau Peninsula, Taipa Island and Coloane Island.These areas are linked by bridges and causeways, creating a cohesive urban expanse that transcends individual points to form a unified urban landscape.This spatial configuration is visually represented in Figure 2, which illustrates the urban continuum and maritime boundaries, providing a comprehensive view of Macau's geographical layout and spatial characteristics.Macau is located in the central part of the Pearl River Estuary, frequently make landfall.The meteorological disasters in Macau a rainstorms.May-November represents the season of frequent ty their activities are characterised by long duration and high freque At the same time, to deeply analyse the specific situation of ty cau, the method of using Macau as the base map is adopted, and visions are conducted in ArcGIS10.8.We used grids as the basic uni invalid grids through precise data processing and finally determin tive grids to be 470, which provided a solid foundation for subse analysis (Figure 2).

CFD wind environment simulation
The wind environment simulation utilises historical typhoo Tropical Cyclone Record List, spanning the years 2000 to 2020.T from the typhoon track network http://typhoon.zjwater.gov.cn/de 17 December 2021), where we manually recorded the transit data o six-hour intervals.These data include the latitude and longitude phoon centre, as well as the wind speed and pressure measurem at the time of the typhoon's occurrence, thereby reflecting near-re Currently, typhoon data is typically collected using anemom tipping-bucket rain gauges, automatic weather stations, weathe ments and CINRAD/CC-type weather radars (Table 1).In the co At the same time, to deeply analyse the specific situation of typhoon disasters in Macau, the method of using Macau as the base map is adopted, and detailed fishing net divisions are conducted in ArcGIS10.8.We used grids as the basic unit of analysis, eliminated invalid grids through precise data processing and finally determined the number of effective grids to be 470, which provided a solid foundation for subsequent typhoon disaster analysis (Figure 2).

CFD wind environment simulation
The wind environment simulation utilises historical typhoon data from the Macau Tropical Cyclone Record List, spanning the years 2000 to 2020.The data were extracted from the typhoon track network http://typhoon.zjwater.gov.cn/default.aspx(accessed on 17 December 2021), where we manually recorded the transit data of the listed typhoons at six-hour intervals.These data include the latitude and longitude coordinates of the typhoon centre, as well as the wind speed and pressure measurements that were obtained at the time of the typhoon's occurrence, thereby reflecting near-real-time conditions.
Currently, typhoon data is typically collected using anemometers with wind cups, tipping-bucket rain gauges, automatic weather stations, weather phenomenon instruments and CINRAD/CC-type weather radars (Table 1).In the context of these historical records, the wind pressure range for typhoons during this period was observed to be between a minimum of 895 Pa and a maximum of 1010 Pa.Wind speed, on the other hand, varied from a low of 10 km/h to a high of 72 km/h.It is important to note that wind pressure, measured in Pascals (Pa), is a physical quantity that indicates the force exerted by the wind on a unit area.A higher wind pressure suggests a stronger typhoon with greater potential for damage.Conversely, wind speed, measured in kilometres per hour (km/h), quantifies the rate at which the wind is moving.The higher the wind speed, the more intense the typhoon and the more significant the impact it can have on structures and the environment.

Equipment Name Explanation
Anemometers with wind cups These measure wind speed by capturing the wind's kinetic energy as it impacts the cups, which rotate at a speed proportional to the wind velocity.

Tipping-bucket rain gauges
These gauges measure rainfall by collecting rain in a bucket that tips when it fills to a certain level, providing a precise measure of precipitation.
Automatic weather stations (AWSs) These are automated systems that measure various weather parameters, including temperature, humidity, atmospheric pressure and wind speed and direction.

Weather phenomenon instruments
These are used to observe and record various weather phenomena such as fog, thunderstorms and visibility conditions.

CINRAD/CC-type weather radars
These radar systems provide detailed information on precipitation distribution, storm movement and intensity, which are crucial for tracking and predicting typhoon paths and impacts.

Identification of typhoon-disaster-vulnerable areas
In the disaster assessment framework system, an increase in exposure and sensitivity increases the region's vulnerability.Thus, a positive correlation is observed amongst exposure, sensitivity and disaster vulnerability.Conversely, increasing adaptability reduces the region's vulnerability.Therefore, this article selects indicators for analysis using the framework of 'exposure-sensitivity-adaptation' combined with the literature (Figure 3).
rotate at a speed proportiona Tipping-bucket rain gauges These gauges measure rainfa a bucket that tips when it fil providing a precise measure Automatic weather stations (AWSs) These are automated system ous weather parameters, inc humidity, atmospheric press and direction.

Weather phenomenon instruments
These are used to observe an weather phenomena such as and visibility conditions.

CINRAD/CC-type weather radars
These radar systems provide on precipitation distribution and intensity, which are cru predicting typhoon paths an

Identification of typhoon-disaster-vulnerable areas
In the disaster assessment framework system, an increase in ex increases the region's vulnerability.Thus, a positive correlation is posure, sensitivity and disaster vulnerability.Conversely, increasin the region's vulnerability.Therefore, this article selects indicators framework of 'exposure-sensitivity-adaptation' combined with th

Data collection and sources
The identification of typhoon-disaster-vulnerable areas is mai data, spatial data and social data collection methods as follows (Ta

Data collection and sources
The identification of typhoon-disaster-vulnerable areas is mainly based on statistical data, spatial data and social data collection methods as follows (Table 2).

Research Method
This section details the research methodology.To visually display the research design, we provide a flow chart (Figure 4) to summarise the entire research process, as follows:

Research Method
This section details the research methodology.To visually display the research design, we provide a flow chart (Figure 4) to summarise the entire research process, as follows:

CFD Wind Environment Simulation
This article mainly studies the wind environment under severe weather caused by typhoon disasters in Macau.Thus, using field measurement and wind tunnel testing research methods is unsuitable.The CFD wind environment simulation method is used for

CFD Wind Environment Simulation
This article mainly studies the wind environment under severe weather caused by typhoon disasters in Macau.Thus, using field measurement and wind tunnel testing research methods is unsuitable.The CFD wind environment simulation method is used for 3D numerical simulation of the entire city of Macau.This approach aims to summarise the characteristics of the wind environment of the Macau typhoon disaster.
ANSYS FLUENT 2019 R1 is used as the simulation software for the typhoon wind ring simulation work in the Macau area.As general-purpose CFD software, FLUENT is powerful and flexible; it can be applied to various complex engineering simulations [47].The framework equations for FLUENT simulation are as follows: where ui represents the average speed in metres per second (m/s); ui' refers to the turbulent pulsation speed in metres per second (m/s); ρ indicates the air density in kilograms per cubic metre (kg/m 3 ); i denotes the buoyancy acceleration in metres per second squared (m/s 2 ); T represents the temperature, whose unit is Kelvin (K); k indicates the turbulent kinetic energy in square metres per square seconds (m 2 /s 2 ); and ε is the turbulent energy dissipation rate in square metres per square seconds (m 2 /s 2 ).Therefore, the k-ε model is implemented in the turbulence closure model of FLUENT simulation.

Principal Component Analysis
Component analysis is a statistical method that converts a group of variables that may have linear correlations into a few linearly independent variables through dimensionality reduction techniques.The set of linearly independent variables after transformation is called the principal components.The basic idea is that multiple original variables (X1, X2, . .., Xp, m ≤ p) with a specific correlation are replaced by a set of linearly independent new comprehensive indicators (Y1, Y2, . .., Ym).At the same time, these new comprehensive indicators reflect the information represented by the original variables to the greatest extent.Moreover, the information between the new comprehensive indicators does not overlap.

Vulnerable Area Identification Model Selection
Goharian et al. (2016) introduced a framework to assess the vulnerability of reservoir systems according to three factors, namely, quantified severity, potential severity and exposure.These factors were estimated using a top-down approach.The vulnerability assessment framework includes sensitivity and adaptive capacity [48].Sensitivity, SoVI and WSACI are estimated using a bottom-up approach and formulated as follows: The research data and process indicate that the vulnerability research applied 'exposuresensitivity-adaptation'.Thus, the second approach is selected as the vulnerability fitting formula.As discussed in the previous literature, vulnerability is generally considered a function of exposure, sensitivity and adaptive capacity.This study adopted a commonly used model to express UDV, as follows: where UDV h is the vulnerability degree of sample h, E h is the exposure degree of sample h, S h is the sensitivity degree of sample h and AC h is the adaptability degree of sample h.UDV h consists of the arithmetic mean of the exposure, sensitivity and adaptability of sample h.It is calculated by adding the three measures and dividing by 3.Each metric value is dimensionless, reflecting the comprehensive vulnerability of sample h to typhoon disasters.

Spatio-Temporal Distribution of Typhoons in Macau
According to the investigation and research, typhoons that make landfall in or affect Macau are mainly prevalent from June to October each year and are affected by typhoon disasters, which is up to five months.Therefore, research on typhoon disasters in Macau is imminent.July to September every year represents the period when typhoons are prevalent in Macau.From 2000 to 2020, the typhoons affecting Macau in these three months accounted for 68% of the total typhoons.The data suggest that typhoon disasters in Macau have prominent seasonal characteristics.At the same time, according to the survey, August exhibits the highest frequency of typhoons.Therefore, according to the seasonal characteristics, disaster identification of typhoon disasters in Macau is considered to deal with frequent typhoon disasters in Macau.

FLUENT Wind Environment Numerical Simulation of Macau Typhoon Disasters
This study delivers a meticulous simulation encompassing the entirety of Macau's urban expanse, meticulously crafted to depict the intricacies of the wind field amidst typhoon conditions.With a strategic approach to simplifying building models, we have adeptly harmonised computational expedience with the fidelity of our simulation, thus capturing a nuanced and comprehensive depiction of Macau's urban topography.
Our simulation model underwent a stringent validation against empirical wind speed data procured from diverse locales throughout Macau (Table 3).A case in point is the Maritime Museum, where our simulation precisely predicted a wind speed of 10.45 m/s, which closely mirrors the actual recorded speed of 11.9 m/s.Further validations were conducted at key sites such as the Memorial Sun Yat-sen Park, where the simulation forecasted a wind speed of 7.496 km/h, nearly indistinguishable from the actual 7.6 km/h, and Coloane, where the simulation's 8.74 km/h was in perfect accord with the actual measurement of 8.6 km/h.At Fortress Hill, the simulation projected a wind speed of 17.49 km/h, which, while slightly above the actual 16.2 km/h, still demonstrates commendable precision (Figure 6).casted a wind speed of 7.496 km/h, nearly indistinguishable from the actua Coloane, where the simulation's 8.74 km/h was in perfect accord with the a ment of 8.6 km/h.At Fortress Hill, the simulation projected a wind speed which, while slightly above the actual 16.2 km/h, still demonstrates comm sion (Figure 6).The study has meticulously examined the interaction between typhoons and Macau's unique regional features by employing simplified building data and leveraging FLUENT software.This approach has facilitated a comprehensive analysis of Macau's 3D urban wind environment and the meticulous simulation of a typhoon disaster tailored to the region's specific characteristics.The integration of these simulated wind speeds into the urban typhoon disaster vulnerability (UDV) assessment framework has notably enhanced our capacity to assess the vulnerability of Macau's districts to the devastating impacts of typhoons.The robust alignment between our simulated and actual wind speeds affirms the model's reliability in identifying areas of heightened risk from severe windstorms and their associated dangers.Furthermore, the detailed examination of Macau's 3D urban wind environment, as depicted in Figures 7-11, yields insightful and actionable intelligence.This analysis is invaluable for formulating strategic disaster preparedness measures and for guiding urban planning to effectively confront the challenges posed by typhoon threats.The study has meticulously examined the interaction between typhoon unique regional features by employing simplified building data and levera software.This approach has facilitated a comprehensive analysis of Mac wind environment and the meticulous simulation of a typhoon disaster region's specific characteristics.The integration of these simulated wind s urban typhoon disaster vulnerability (UDV) assessment framework has not our capacity to assess the vulnerability of Macau's districts to the devastat typhoons.The robust alignment between our simulated and actual wind the model's reliability in identifying areas of heightened risk from severe w their associated dangers.Furthermore, the detailed examination of Mac wind environment, as depicted in Figures 7-11, yields insightful and act gence.This analysis is invaluable for formulating strategic disaster measures and for guiding urban planning to effectively confront the challe typhoon threats.The study has meticulously examined the interaction between typhoon unique regional features by employing simplified building data and levera software.This approach has facilitated a comprehensive analysis of Mac wind environment and the meticulous simulation of a typhoon disaster region's specific characteristics.The integration of these simulated wind s urban typhoon disaster vulnerability (UDV) assessment framework has not our capacity to assess the vulnerability of Macau's districts to the devastat typhoons.The robust alignment between our simulated and actual wind the model's reliability in identifying areas of heightened risk from severe w their associated dangers.Furthermore, the detailed examination of Mac wind environment, as depicted in Figures 7-11, yields insightful and act gence.This analysis is invaluable for formulating strategic disaster measures and for guiding urban planning to effectively confront the challe typhoon threats.

Numerical Simulation of a Typhoon Disaster in Macau Using FLUENT for W ment
According to the evaluation model, the specific parameters of exposu and adaptability of 470 grids are calculated (Figure 2), and the correspondi UDV values are obtained.This study further adopts the standard deviat objectively and reasonably discuss the statistical distribution of each catego ple.The standard deviation category shows the difference between a feat value and the mean and the degree to which the behaviour differs from UDV value is divided into five vulnerability levels: very low, low, med heavy.
In accordance with the evaluation of UDV in Macau, the vulnerability v difference are analysed.The different effects of grids on the three dimen (exposure, sensitivity and adaptability) are highlighted by analysing the ponents.When interpreting the dimensions and principal components PCA, the criterion of 'eigenvalues greater than 1' is used.These principal co

Numerical Simulation of a Typhoon Disaster in Macau Using FLUENT ment
According to the evaluation model, the specific parameters of e and adaptability of 470 grids are calculated (Figure 2), and the corresp UDV values are obtained.This study further adopts the standard d objectively and reasonably discuss the statistical distribution of each c ple.The standard deviation category shows the difference between value and the mean and the degree to which the behaviour differs f UDV value is divided into five vulnerability levels: very low, low, heavy.
In accordance with the evaluation of UDV in Macau, the vulnerab difference are analysed.The different effects of grids on the three d (exposure, sensitivity and adaptability) are highlighted by analysing ponents.When interpreting the dimensions and principal compon

Numerical Simulation of a Typhoon Disaster in Macau Using FLUENT for Wind Environment
According to the evaluation model, the specific parameters of exposure, sensitivity and adaptability of 470 grids are calculated (Figure 2), and the corresponding normalised UDV values are obtained.This study further adopts the standard deviation method to objectively and reasonably discuss the statistical distribution of each category in the sample.The standard deviation category shows the difference between a feature's attribute value and the mean and the degree to which the behaviour differs from the mean.The UDV value is divided into five vulnerability levels: very low, low, medium, high and heavy.
In accordance with the evaluation of UDV in Macau, the vulnerability value and level difference are analysed.The different effects of grids on the three dimensions of UDV (exposure, sensitivity and adaptability) are highlighted by analysing the principal components.When interpreting the dimensions and principal components obtained from PCA, the criterion of 'eigenvalues greater than 1' is used.These principal components and important indicators reveal the critical indicators of the sample UDV clearly and accurately in each dimension.

Exposure Analysis
PCA is used for the data of four exposure indicators: the proportion of land use and buildings, wind speed, average elevation and population density (Table 4).Moreover, the two principal component solutions explain 59.43% of the variance in the dataset.The first principal component includes the building footprint, average elevation and population density, integrating the current urban development and population expansion levels in Macau and therefore referred to as the urban development factors.
The second principal component comprises the average wind speed, which reflects the influence of natural disaster factors under typhoon disasters and is therefore named the natural environment factor.The exposure of Macau to typhoon disasters can be analysed more intuitively using the principal component table.Meanwhile, based on the data for each principal component, the principal component scores of 470 grids were divided into five levels using the natural-break method in ArcGIS.
These high exposures are confined to ultra-high population densities and high building footprints, resulting in a large number of residents and fewer open areas.This condition exacerbates the exposure to typhoon hazards in Macau.In addition, the average typhoon wind speed causes severe exposure of these grids, as global climate change increases likelihood of natural disasters.Macau, a high-density coastal city, is more vulnerable to typhoon hazards and more likely to cause economic losses and human casualties.Figure 12 shows that the extremely high exposure of the Chopsticks base area is caused by the high population density.
ISPRS Int.J. Geo-Inf.2024, 13, x FOR PEER REVIEW 14 of 22 The first principal component includes the building footprint, average elevation and population density, integrating the current urban development and population expansion levels in Macau and therefore referred to as the urban development factors.
The second principal component comprises the average wind speed, which reflects the influence of natural disaster factors under typhoon disasters and is therefore named the natural environment factor.The exposure of Macau to typhoon disasters can be analysed more intuitively using the principal component table.Meanwhile, based on the data for each principal component, the principal component scores of 470 grids were divided into five levels using the natural-break method in ArcGIS.
These high exposures are confined to ultra-high population densities and high building footprints, resulting in a large number of residents and fewer open areas.This condition exacerbates the exposure to typhoon hazards in Macau.In addition, the average typhoon wind speed causes severe exposure of these grids, as global climate change increases likelihood of natural disasters.Macau, a high-density coastal city, is more vulnerable to typhoon hazards and more likely to cause economic losses and human casualties.Figure 12 shows that the extremely high exposure of the Chopsticks base area is caused by the high population density.

Sensitivity Analysis
In this study, PCA analysis was used for five sensitivity indicators: the proportion of the elderly population over 65 years old, the proportion of the female population, the proportion of the population under 14 years old, the total population and the age of buildings.In addition, the three principal components cumulatively explained 77.29% of the vari-

Sensitivity Analysis
In this study, PCA analysis was used for five sensitivity indicators: the proportion of the elderly population over 65 years old, the proportion of the female population, the proportion of the population under 14 years old, the total population and the age of buildings.In addition, the three principal components cumulatively explained 77.29% of the variance in the dataset (Table 5).The first principal component is the proportion of the elderly population over 65 years old and the proportion of the female population, collectively reflecting the current situation of vulnerable groups in Macau.The second principal component is the proportion of the population under 14 years old and the total population, reflecting the impact of the vulnerable population factor under a typhoon disaster.It is, therefore, referred to as the vulnerable population factor.The third principal component is the age of the buildings, integrally reflecting the current situation of Macau; it is termed the building factor.In accordance with the principal component table, the sensitivity of Macau to typhoon disasters is analysed more intuitively.The principal component scores of 470 grids are divided into five levels using the natural-break method in ArcGIS.
These high sensitivities are exacerbated by the high proportion of the elderly population, the proportion of the female population and the proportion of the population under 14 years old.As a result, many inhabitants are considered vulnerable, exacerbating the sensitivity of Macau to typhoon hazards.In addition, the total population is another reason for the severe exposure of these grids.The higher population size of Macau leads to a higher number of people at risk during typhoon hazards and more vulnerable to receiving typhoon hazards.Figure 13 shows that St. Anthony's Parish is extremely sensitive given the high population density within the area.The first principal component is the proportion of the elderly population over 65 years old and the proportion of the female population, collectively reflecting the current situation of vulnerable groups in Macau.The second principal component is the proportion of the population under 14 years old and the total population, reflecting the impact of the vulnerable population factor under a typhoon disaster.It is, therefore, referred to as the vulnerable population factor.The third principal component is the age of the buildings, integrally reflecting the current situation of Macau; it is termed the building factor.In accordance with the principal component table, the sensitivity of Macau to typhoon disasters is analysed more intuitively.The principal component scores of 470 grids are divided into five levels using the natural-break method in ArcGIS.
These high sensitivities are exacerbated by the high proportion of the elderly population, the proportion of the female population and the proportion of the population under 14 years old.As a result, many inhabitants are considered vulnerable, exacerbating the sensitivity of Macau to typhoon hazards.In addition, the total population is another reason for the severe exposure of these grids.The higher population size of Macau leads to a higher number of people at risk during typhoon hazards and more vulnerable to receiving typhoon hazards.Figure 13 shows that St. Anthony's Parish is extremely sensitive given the high population density within the area.6).This component contained two factors for vegetation

Analysis of UDV Vulnerability in Macau
The UDV value of vulnerability in Macau was calculated using the Moran's index.Its value was 0.155 with a z-score of 15.76, and the probability of randomly generating this clustering pattern was less than 1%.Therefore, this value exhibits spatial autocorrelation.
Meanwhile, ArcGIS is used to analyse the local autocorrelation of vulnerability UDV values in Macau and LISA regional aggregation analysis to visualise the characteristics of its spatial autocorrelation (Figure 15).The figure shows the four main types of spatial autocorrelation: high-high, low-low, high-low and low-high, as well as a blank nonsignificant grid.By contrast, high-low and low-high show high or low values, with spatial outliers reflecting negative local spatial autocorrelation.The figure illustrates that St. Anthony's Parish exhibits a high value of aggregation centre, indicating that the UDV value of the area requires attention and improvement.This area is undoubtedly the most vulnerable when a typhoon disaster strikes.
ISPRS Int.J. Geo-Inf.2024, 13, x FOR PEER REVIEW 17 ArcGis 10.7 was used to calculate the mean and standard deviation of UDV va Visualisation was conducted by highlighting the UDV values above and below the m using the colour bands (Figure 16).In this study, the similarity statistics of exposure, sitivity and adaptive capacity were calculated for each grid and merged into the categ closest to the squared Euclidean distance.They were defined as the Macau vulnerab UDV types: habitable, safe, generally safe, average, highly vulnerable, very vulne and hazardous.In addition, the two subtypes with the highest frequency in each U type were selected as samples.On this basis, the grid was divided into 19 subtyp clarify the principal dimensions and components of UDV levels in different grids (T 7).ArcGis 10.7 was used to calculate the mean and standard deviation of UDV values.Visualisation was conducted by highlighting the UDV values above and below the mean using the colour bands (Figure 16).In this study, the similarity statistics of exposure, sensitivity and adaptive capacity were calculated for each grid and merged into the categories closest to the squared Euclidean distance.They were defined as the Macau vulnerability UDV types: habitable, safe, generally safe, average, highly vulnerable, very vulnerable and hazardous.In addition, the two subtypes with the highest frequency in each UDV type were selected as samples.On this basis, the grid was divided into 19 subtypes to clarify the principal dimensions and components of UDV levels in different grids (Table 7).has 105 grids, of which 27 are divided into three subtypes.The average type has 185 grids, of which 25 are divided into three subtypes.The highly vulnerable type has 118 grids, of which 23 are divided into three subtypes.The very high vulnerability type has 27 grids, of which five are divided into five subtypes.The hazardous type has two grids with On basis of the main dimensions of each UDV type, the planning strategy is recommended as follows.
Habitable, safe: This grid has fewer exposure and sensitivity factors.Although the adaptive capacity is not high in this type of grid, it causes lower exposure and sensitivity.These areas are mainly characterised by high sensitivity caused by a high proportion of the elderly population, a high proportion of the female population, a high proportion of the population under 14 years old and a high proportion of vulnerable people in the total population.Urban planners should establish support organisations for vulnerable people in specific areas.
General safety type, general type: Some grids are composed of moderate exposure and sensitivity.In this grid type, exposure sensitivity is affected by some adverse factors.Theoretically, the government should adopt various urban planning measures to reduce UDV, such as controlling the expansion of construction land, reducing the building footprint and planning urban wind corridors.However, considering the actual situation of Macau as the most densely populated region in the world, the population density and buildings in large amounts cannot be controlled.Therefore, this type of planning should focus on managing land use in Macau, reducing the expansion of building sites and planning for managing illegal sites and new sites.When the plan effectively reduces exposure and sensitivity, the safe grid strategy can be used.
Highly vulnerable, extremely highly vulnerable and hazardous: Involves areas with high or very high exposure and sensitivity but low adaptability.This grid type is affected by several unfavourable factors in terms of exposure, sensitivity and adaptability.Firstly, on top of the previously mentioned control for the land use of Macau, the elevation issue must be considered.Secondly, a survey should be conducted every year prior to the typhoon season in May for buildings with high service life.A decision should be made to demolish, rebuild or reinforce the buildings according to the situation.Dangerous situations due to old buildings, such as sliding, falling off and building collapse, must be prevented.Finally, the vegetation cover must be strengthened because vegetation plays a significant role in regional typhoons.Plans to construct medical resources in the city, in a specific range, must be developed to ensure a certain number of medical institutions.

Conclusions
Case studies on spatial vulnerability assessment systems are common worldwide [9][10][11]49].However, data from the previous literature cannot sufficiently demonstrate all the three dimensions [8,10].In this study, the 12 cumulative indicators were completed with categorisation using three dimensions.This condition allows city managers to tailor the research results to different attributes and provides a basis for selecting indicators for future studies.Furthermore, the innovative addition of 'average wind speed' as an exposure indicator provides a critical indicator reference for future typhoon hazard vulnerability studies.
In the past, few studies have been conducted to analyse the spatial vulnerability of individual cities because another study has conducted macro-or micro-level investigations [8,9,50].However, local governments globally tend to approach disaster policymaking from an urban perspective [51].Therefore, this study not only filled the gaps in identifying urban wind hazard spatial vulnerability areas in Macau but also provides a reliable reference for future urban studies.
The data selected for this study are those generated from CFD simulations of typhoons in Macau.The vulnerability zones generated from these simulations are informative based on the final results.Therefore, an essential theoretical basis for whether CFD simulation data can be used to identify future wind hazard spatial vulnerability zones is provided.
In other words, with the appropriate use of CFD simulation software, the generated data are expected to be an essential data reference or supplement for identifying wind hazard spatial vulnerability zones.
Generally, the UDV evaluation index system for typhoon vulnerability in Macau is used to identify the spatial characteristics of 470 grid samples and the potential disaster vulnerability relationships in the process.A UDV evaluation pointer system is proposed to explore its principal components, and a measurement model using a series of factors is established in three dimensions: exposure, sensitivity and self-adjustment capacity.This approach aims to reveal the factors that strongly influence the typhoon vulnerability of Macau and provide targeted suggestions for future prevention and adaptation to extreme disaster weather in Macau.
The results show that 31.27% of the grid belongs to the severe vulnerable zone (SVZ) (Figure 17).UDV has six principal components with high positive and negative correlations, mainly influenced by population density, building occupancy ratio, elevation, wind speed and other indicators.Moreover, on this basis, countermeasures are proposed, as follows: ISPRS Int.J. Geo-Inf.2024, 13, x FOR PEER REVIEW  1.
Urban planners should optimise the population structure of Macau and set the building volume ratio standards in the planning and design to avoid the increase in disaster risk in this area due to the over-concentration of the population.The floor area ratio is proposed for high-rise and super-high-rise buildings based on the detailed control plans for various types of residential land prepared under the current urban planning regulation system and considering the situation in Macau.The proposed floor area ratio for high-rise buildings is 2.5-3.5, and the floor area ratio for super high-rise buildings above 19 stories is 3.6-5.

2.
Support organisations must be established for disadvantaged individuals in a specific range to avoid increasing the number of injuries during typhoons.

3.
A plan must be developed for controlling land use in Macau, reducing construction land expansion and managing illegal and unused land.4.
We must strengthen the vegetation cover, significantly affecting the geographical impact of typhoons.Medical resources in the city must be established in a specific range to ensure a certain number of medical institutions.

Figure 1 .
Figure 1.Location of Macau in China.

Figure 1 .
Figure 1.Location of Macau in China.

Figure 2 .
Figure 2. Macau's geographical composition and grid division.Macau is located in the central part of the Pearl River Estuary, where tropical cyclones frequently make landfall.The meteorological disasters in Macau are mainly typhoons and rainstorms.May-November represents the season of frequent typhoons in Macau, and their activities are characterised by long duration and high frequency.At the same time, to deeply analyse the specific situation of typhoon disasters in Macau, the method of using Macau as the base map is adopted, and detailed fishing net divisions are conducted in ArcGIS10.8.We used grids as the basic unit of analysis, eliminated invalid grids through precise data processing and finally determined the number of effective grids to be 470, which provided a solid foundation for subsequent typhoon disaster analysis (Figure2).
(a) Statistical data: They mainly reflect Macau's construction status and population data, which are obtained from the statistical geographic information system of the Statistics and Census Bureau of Macau.(b) Spatial data: They describe objects to represent the shape, size, location and distribution characteristics of spatial entities, including vegetation coverage, elevation and population density.The vegetation coverage database is calculated from the 2020 Landsat image of Macau according to the vegetation coverage formula.The Cartography and Cadastre Bureau of the Macau Special Administrative Region provides elevation data.Population density data are obtained from the World Pop 1 km dataset (https://www.worldpop.org/accessed on 10 February 2022).(c) Social information: It mainly reflects the status of medical resources in Macau.In this article, it is represented by the points of interest (POIs) of medical resources, the source of which is Baidu Map.
ISPRS Int.J. Geo-Inf.2024, 13, x FOR PEER REVIEW 8 of 22 (a) Statistical data: They mainly reflect Macau's construction status and population data, which are obtained from the statistical geographic information system of the Statistics and Census Bureau of Macau.(b) Spatial data: They describe objects to represent the shape, size, location and distribution characteristics of spatial entities, including vegetation coverage, elevation and population density.The vegetation coverage database is calculated from the 2020 Landsat image of Macau according to the vegetation coverage formula.The Cartography and Cadastre Bureau of the Macau Special Administrative Region provides elevation data.Population density data are obtained from the World Pop 1 km dataset (https://www.worldpop.org/accessed on 10 February 2022).(c) Social information: It mainly reflects the status of medical resources in Macau.In this article, it is represented by the points of interest (POIs) of medical resources, the source of which is Baidu Map.

4 . Result 4 . 1 .
Distribution of Typhoons from 2000 to 2020 4.1.1.Spatial Distribution of Typhoons in Macau This study provides 100 typhoon disasters that have affected Macau in the past 20 years through the list of tropical cyclone records.Moreover, relevant typhoon disasters from 2000 to 2020 are selected from the typhoon track website (http://typhoon.zjwater.gov.cn/default.aspx)(accessed on 19 January 2022).Furthermore, TXT files are used to intercept and record the latitude and longitude coordinates and wind speed.Typhoons transmit data, such as wind pressure and wind direction.Then, they are imported into ArcGIS and converted into a shapefile layer for density analysis to find the spatial law and action path of typhoon disasters affecting Macau.The typhoon track map (Figure 5) shows that most of the typhoon tracks in the past 20 years have been relatively regular.Red cycle in the figure indicates Macau location in different administrative level (e.g. of national level).The typhoons' operating density and directional trends are summarised on the basis of the particularity and randomness of typhoon track activities.The statistics indicate that the typhoons affecting Macau in the past 20 years can be roughly divided into three different typhoon transit paths.Firstly, the trajectory of Path 1 moves from east to west, moving westward after landing in the eastern coastal area of Macau.The second path represents the track with the highest frequency of occurrence amongst the typhoon tracks affecting Macau.This path made landfall on the southeast coast of Macau and moved northwestward.According to statistics, most of them affected St. Francis Parish.Finally, the third path moves from the southern coast of Macau to the north and affects Guangdong Province and other regions.

4 . 1 .
Distribution of Typhoons from 2000 to 2020 4.1.1.Spatial Distribution of Typhoons in Macau This study provides 100 typhoon disasters that have affected Macau in the past 20 years through the list of tropical cyclone records.Moreover, relevant typhoon disasters from 2000 to 2020 are selected from the typhoon track website (http://typhoon.zjwater.gov.cn/default.aspx) (accessed on 19 January 2022).Furthermore, TXT files are used to intercept and record the latitude and longitude coordinates and wind speed.Typhoons transmit data, such as wind pressure and wind direction.Then, they are imported into ArcGIS and converted into a shapefile layer for density analysis to find the spatial law and action path of typhoon disasters affecting Macau.The typhoon track map (Figure 5) shows that most of the typhoon tracks in the past 20 years have been relatively regular.Red cycle in the figure indicates Macau location in different administrative level (e.g. of national level).The typhoons' operating density and directional trends are summarised on the basis of the particularity and randomness of typhoon track activities.The statistics indicate that the typhoons affecting Macau in the past 20 years can be roughly divided into three different typhoon transit paths.Firstly, the trajectory of Path 1 moves from east to west, moving westward after landing in the eastern coastal area of Macau.The second path represents the track with the highest frequency of occurrence amongst the typhoon tracks affecting Macau.This path made landfall on the southeast coast of Macau and moved northwestward.According to statistics, most of them affected St. Francis Parish.Finally, the third path moves from the southern coast of Macau to the north and affects Guangdong Province and other regions.

Figure 6 .
Figure 6.Simplified verification diagram of wind simulation building.

Figure 6 .
Figure 6.Simplified verification diagram of wind simulation building.
SPRS Int.J. Geo-Inf.2024, 13, x FOR PEER were collected from the Macau Special Administrative Region Governm on 18 October 2022.

Figure 7 .
Figure 7. Simulated wind speed map of a typhoon in Macau.

Figure 8 .
Figure 8. Partial map 1 of simulated wind speed of a typhoon in Macau.

Figure 7 .
Figure 7. Simulated wind speed map of a typhoon in Macau.

Figure 7 .
Figure 7. Simulated wind speed map of a typhoon in Macau.

Figure 8 .
Figure 8. Partial map 1 of simulated wind speed of a typhoon in Macau.

Figure 8 .
Figure 8. Partial map 1 of simulated wind speed of a typhoon in Macau.

Figure 8 .
Figure 8. Partial map 1 of simulated wind speed of a typhoon in Macau.

Figure 9 .
Figure 9. Partial map 2 of simulated wind speed of a typhoon in Macau.

Figure 9 .Figure 10 .
Figure 9. Partial map 2 of simulated wind speed of a typhoon in Macau.

Figure 11 .
Figure 11.Panorama of simulated wind speed of a typhoon in Macau.

Figure 10 .Figure 10 .
Figure 10.Partial map 3 of simulated wind speed of a typhoon in Macau.

Figure 11 .
Figure 11.Panorama of simulated wind speed of a typhoon in Macau.

Figure 11 .
Figure 11.Panorama of simulated wind speed of a typhoon in Macau.

Figure 13 .
Figure 13.Spatial distribution of sensitivity.4.3.3.Adaptive Analysis In this study, PCA was used for two adaptation indicators: vegetation cover and number of medical resources.The principal component explained about 50.87% of the variance in the dataset (Table6).This component contained two factors for vegetation

Table 1 .
Instruments and observing stations for typhoon data collection.

Table 2 .
Comprehensive calculation table of indicators.

Table 2 .
Comprehensive calculation table of indicators.

Table 3 .
Actual wind speed at the Macau Meteorological Monitoring Station.
Note: Data were collected from the Macau Special Administrative Region Government at 9.30 pm on 18 October 2022.

Table 3 .
Actual wind speed at the Macau Meteorological Monitoring Station.

Table 4 .
Eigenvalues and variances of each principal component of exposure ability.

Table 5 .
Eigenvalues and variances of each principal component of sensitivity ability.