Assessment of urbanization impact on cultural heritage based on a risk-based cumulative impact assessment method

Urbanization is inevitable in both developing and developed countries. However, this growth and transformation of the urban area can pose a significant threat to urban cultural heritage, which is a sensitive component of the urban environment. As cities modernize and change, a risk of irreparable loss of cultural heritage exists. Therefore, taking steps to protect and preserve these sites for posterity is crucial. To better protect urban cultural heritage, decision-makers must rapidly assess the impact of urbanization on cultural heritage while maintaining a balance between cultural heritage preservation and urban growth. This study developed a risk-based cumulative impact assessment (CIA) method that integrates a set of quantifiable indicators to address these issues. This method generates standardized assessment results based on game theory and set pair analysis theory. In contrast to other CIA methods, this risk-based approach effectively manages scientific complexity and uncertainties, thus enhancing the quality of the assessment results. This method was applied to 21 classical gardens in Suzhou, China. The results show that Suzhou classical gardens are impacted by urbanization. This study confirmed that adopting a risk-based CIA method that considers the impact of adverse urbanization on cultural heritage sites is an efficient approach. Successful implementation of the proposed method can provide decision-making support for different types of cultural heritage in other areas.


Introduction
In the current context of rapid economic and social development, urbanization is occurring at a never-before-seen rate, and the surroundings of urban settings are changing as a result of the speedy and usually unchecked development.However, urban features are constantly changing causing challenges in the conservation of urban cultural heritage.
It is universally acknowledged that urbanization is the primary factor contributing to the decline in value and integrity of cultural heritage particularly in historical areas where preparation is still lacking [1][2][3].As an integral component of the urban area, cultural heritage may be threatened by urbanization process.This may result in the loss or distortion of material properties that have cultural heritage value, thus preventing sustainable development and conservation.Notably, in recent years, many investigations and reports have demonstrated increasing levels of threats from urbanization on urban cultural heritage sites, such as building and development projects [3][4][5], transportation infrastructure [6], and environmental pollution [7].
Urban cultural heritage is an essential component of urban areas and can reflect the historical development of cities. Maintaining the connection between the past and the future is greatly dependent on the preservation and sustainable development of cultural heritage.Therefore, supporting the sustainable development of urban cultural heritage and addressing the threats faced by it requires a suitable and applicable assessment tool.In this context, the Heritage Impact Assessment (HIA) was introduced to identify and analyze the possible adverse effects of anthropogenic activities on cultural heritage sites [3].However, the use of an assessment tool cannot be limited to only supporting the management of urban cultural heritage, but it must also have utility as a method to help protect the critical values embodied in cultural heritage.The HIA has been successfully applied and improved upon in several studies [3,4,[8][9][10].However, this approach has a drawback because the current HIA method only considers the impact of a single development project and fails to account for the combined and cumulative impacts of multiple projects.Urbanization impact on urban cultural heritage is a complex problem that requires evaluation of cumulative impacts, which can be derived from different projects.To address this issue, further research is required to investigate cumulative impact assessment (CIA) techniques that are applicable to cultural heritage sites.
CIA has long been acknowledged as a crucial component of impact assessments [11].It is defined as a comprehensive evaluation of the combined effects of human actions and natural processes on the evaluation subject across space and time [12][13][14].The assessment of cumulative impacts has been used in Canada and the United States for many years, where its procedures are incorporated into the Environmental Impact Assessment and Strategic Environmental Assessment procedures [15,16].CIA has been explored by researchers in many fields of science, including oceanography and marine area protection [13,[17][18][19][20], toxicology and ecotoxicology [21,22], human health risk assessment [23], forestry [24,25], fish and wildlife management [26][27][28], freshwater and watershed management [29,30] and environmental justice [31,32].Compared to the HIA, CIA is more comprehensive because it considers all the potential environmental factors across time.However, no previous study has considered the application of a CIA to cultural heritage.
Characteristics of an urban environment, in particular its vast openness, high interconnection of economic and social factors, and large variability and uncertainty, generate additional complexities and challenges for CIA.However, a thorough and open framework is frequently absent from the operationalization and integration of CIA into decision-making processes.A risk-based CIA framework that divides the CIA method into risk identification, risk analysis, and risk evaluation can better analyze such complexities [13].Moreover, cumulative impacts are defined as the combined effects of impacts and growth over time, including past, present, and foreseeable future pressures [33].Currently, existing studies lack a prediction of possible future risks.Therefore, the risk-based CIA method is a promising approach for aligning the cumulative impact assessment with the future risks [34].It also helps develop preventive measures, adaptive surveillance, and impact reduction measures during the urbanization projects [35].
Given the distinctive characteristics and unparalleled value of urban cultural heritage sites, this study aimed to present a risk-based CIA framework that incorporates a set of indexes to assess the impact of urbanization on urban cultural heritage.The proposed method is designed to provide scientific, quantitative, and standardized methods for assessing the impact of urbanization on cultural heritage sites based on objective and quantitative data.

Overview of the method and design
Cumulative impacts are related to the complexity of human endeavors and initiatives that result in numerous pressures that build up and impact cultural heritage sites [36].To address this complexity, the CIA was formulated in a risk-based setting.
Risk is the consequence of the interaction between the impacts of an event or hazard and the associated probability that it will occur [37].The consequences are negative impacts caused by social, political, economic, and environmental factors.In risk analysis, the impacts are often expressed in terms of exposure and vulnerability.Exposure is the totality of objects present in hazardous areas that are subject to potential losses.Vulnerability refers to the characteristics, environments, systems, or resources that make an individual susceptible to the negative impact of a hazard [38].Hence, the CIA method should reveal the probability of cumulative impact, exposure of cumulative impact, and vulnerability of cumulative impact on cultural heritage.
where CP is the cumulative impact probability, CE is the cumulative impact exposure and CI is the cumulative impact vulnerability.
This study presents a risk-based CIA framework that integrates a set of indicators to assess urbanization's impact on cultural heritage sites (Fig. 1).This framework is embedded in the International Organization for Standardization (ISO) risk management process, which (1) Cumulative impact = f (CP, CE, CV ) includes three parts: risk identification, risk analysis, and risk evaluation.
First, risk identification is the process of identifying and assessing risk sources and their potential impacts.Creating cause-and-effect links or risk pathways to describe the vulnerability of cultural heritage to the effects of urbanization is a crucial challenge in risk identification [11].In this step, an index system including cumulative impact probability, cumulative impact exposure, and cumulative impact vulnerability is constructed, and the evaluation criteria of each index are determined by Jenks classification.
Second, risk analysis comprises comprehending the nature of risk and determining the weight share between risk sources [39].In this stage, the analytic hierarchy process (AHP) and entropy weight (EW) methods are used to calculate the subjective and objective weights of each index, and game theory is used to determine the combined weight.Third, risk evaluation identifies specified risk criteria and determines the risk level.The decision concerning risk treatment is aided by risk evaluation, which also necessitates a performance review of new measures [40].This step introduced set pair analysis (SPA) theory to determine the cumulative impact risk level.

Cumulative impact assessment index system
A cumulative impact assessment index system was constructed using the risk-based CIA framework (Table 1).Cumulative impact probability refers to the probability of future urbanization having an impact on the cultural heritage site.Cumulative impact vulnerability refers to the current impact of urbanization on the cultural heritage site and its vulnerability.
The cumulative impact probability indices are measured as follows: A) Population density (D1) has helped examine population growth as it indicates how cultural heritage sites are prone to urban growth.The distance from cultural heritage sites to densely populated areas has been found to be proportional to the possibility of interference from human activities [41].The index refer to the population density of the subdistrict where the heritage sites are located.B) The new construction land area (D2) is used to examine future construction projects around the cultural heritage sites.Construction projects can generate mechanical shock, pollute the environment, and damage heritage sites.This index is measured by the ratio of the new construction land area within the buffer zone to the buffer zone area.C) The percentage of surrounding commercial land use (D3) is used to measure the possibility of interference from human activities.The higher the proportion of commercial land around the heritage site, the more human activities and the more intensive the potential urban construction activity.This index is measured by the ratio of the commercial land use area within the buffer zone to the buffer zone area.
The cumulative impact exposure indices are measured as follows: A) The area of cultural heritage (D4) is an index that reflects a heritage site's level of exposure to the impact of urbanization.The larger the area of the cultural heritage site, the higher the exposure.B) The percentage of open space (D5) is an index that reflects a heritage site's level of exposure to the impact of urbanization.Open spaces are exposed to urban environments and are more sensitive to the visual impact and noise caused by urban construction.C) Accessibility (D6) is used to measure the accessibility of cultural heritage sites through transportation networks.Mechanical shock and pollutants from heavy traffic on streets often exceed this limit, posing a risk to the cultural heritage site.Noise can also degrade the tranquility of the site [41].This index can be calculated using the ArcGIS software based on the origin destination (OD) cost matrix analysis tool.D) Distance from the subway station (D7) is an index used to measure the distance between a cultural heritage site and subway stations.Typically, more commercial and urban construction activities occur around metro stations.E) The distance from subway tracks (D8) is an index used to measure the distance between cultural heritage sites and subway tracks.Both the construction Table 1 Cumulative impact assessment index system of urban cultural heritage and operation of the metro system affect cultural heritage sites [42].
The cumulative impact vulnerability indexes are measured as follows: A) Harmony index of neighboring building forms (D9) is used to indicate the degree of harmony reflected in the relationship between human-made facilities, buildings, infrastructure, and heritage sites.High harmony rarely reduces landscape integrity [43].This indicator is obtained from the ratio of the area of incongruous building forms within the construction control area to the total construction control area.B) Heritage site history (D10)-that is, the antiquity of a cultural heritage site-is used as a proxy for heritage value.The longer the heritage site's history, the higher its value and vulnerability to disasters.C) The heritage protection level (D11) reflects the value of the heritage site, the higher the protection level, the higher the value.Five classes of heritage protection levels exist, namely heritage on the World Heritage List, national-level protected heritage, provinciallevel protected heritage, city-level protected heritage, and non-protected heritage units, which are assigned the values of 9, 7, 5, 3, and 1, respectively.

Grading standard of assessment indexes
The risk standards connected to the indices were divided into five classes per the Jenks classification (Table 2).It is a data categorization method that is designed to optimize the arrangement of a set of values into natural classes.This approach clusters data into groups that minimize the within-group variance and maximize the betweengroup variance [44,45].Natural groupings, inherent in the data, serve as the foundation for natural breaks.The best group of similar values and those that emphasize class distinctions are used to determine class breaks [46].

Risk analysis Combination weight calculation
Determining the appropriate weights for the cumulative impact assessment index system is an essential step in the CIA method.There are two weighting methods, the subjective and the objective, both with advantages and drawbacks [47].The use of the subjective weighting method can reveal experts' subjective opinions and experiences.
Objective weights based on mathematical theory can avoid subjectivity.The objective weight method-named EW-has been applied in several studies [50,51].However, because the objective weight approach can only analyze the value change between indices, the information relevance between indices cannot be assessed.Therefore, game theory was applied in this study to combine the advantages of each of the two weighting methods.

Subjective weight based on the AHP
The AHP is a decision-making strategy that uses a systematic method to evaluate and integrate the effects of many factors [52].By creating a pairwise comparison judgment matrix between the parameters, this wellknown decision-analysis technique enables an analyst to generate weights.Ten experts were invited to participate in this study and complete a questionnaire to construct the judgment matrix.

Objective weight based on EW
Based on the degree of dispersion of each index value, EW may objectively calculate the index weight [53].
Four main steps exist in the entropy method calculation

Combination weight based on game theory
Game theory is an essential subfield of operations research that focuses on the interaction between competing alternatives and the rational behavior of various decision-makers, including their decision equilibrium when their actions affect each other.In the game, all parties make decisions to minimize their losses or maximize their interests; this requires all parties to find a balanced combination to maximize common interests.Consequently, it is possible to make an individual yet collaborative decision that optimizes all stakeholders' projected returns [47].The procedure is as follows: 1) The basic combination of weights is carried out by collecting and using results from different weighting methods W = {w 1 , w 2 , . . ., w m } , where wm is the weight vector of weights for m indices, and W is obtained by the linear combination coefficient ai as follows: 2) By optimizing the weight coefficient ak, the most satisfactory W can be found, which could minimize the deviation between W and wk.
According to the differentiation properties of the matrix, the condition for the optimal first-order derivative in Eq. ( 3) is as follows: Then, to facilitate analysis, the matrix form of Eq. ( 4) is expressed as follows: (2) 3) To obtain a*, the weight coefficient is calculated and normalized as follows: 4) Finally, the most satisfied combination weight can be reached as follows:

Risk evaluation Set pair analysis theory
In 1989, Zhao created the SPA theory as an innovative system analysis approach that aims to explore the certainty and uncertainty of a particular system [54].SPA theory is used to study the connection between certainty and uncertainty by treating two sets with certainty and uncertainty as a set pair and describing the sameness, difference, and opposition of the set pair in terms of the degree of connection [55].

Calculation of single correlation degree
Suppose the set of evaluation factors, A, and the set of evaluation criteria, B, constitute the set pair M, which is analyzed to obtain N characteristics; the set pair is described by the degree of association u: where S/N denotes the homogeneity of the set pair; F/N denotes the degree of difference of the set pair; P/N denotes the opposition degree of the set pair, denoted as a, b, c; a + b + c = 1; i denotes the difference degree coefficient, taking values between [-1,1]; and j is the opposition degree coefficient, generally assumed to be -1.
Depending on the study object, the degree of connection can be diversified.The expression for the k degree of connection is as follows: When studying the risk level of the cumulative impact, the risk is divided into five levels, whose single-indicator correlation can be expressed as follows: The indexes were divided into positive and negative indicators depending on their characteristics, and their correlation degrees were calculated as follows: Positive index correlation degree calculation:

Negative index correlation degree calculation:
where S1 to Sk are the upper thresholds for the respective risk levels and xm is the actual value of the evaluation indexes.

Calculation of the integrated correlation degree
After calculating the single-index correlation degree of each index by Eqs.11 and 12, the integrated correlation degree of the set pair H (A, B) can be calculated as follows: where w m is the weight of the index and µ AB is the integrated correlation degree of the risk index on the corresponding risk level.

Determining the risk level of the cumulative impact
Because the value of the coefficient of variation varies, this study introduced a confidence criterion to determine the risk level of the evaluation unit to avoid a discussion of the coefficient of variation and reduce the complexity of the calculation process [56].In this study, the cumulative impact risk was classified into five risk levels: low-risk, lower-risk, medium-risk, higher-risk, and high-risk, labeled as I, II, III, IV, and V, respectively.(11) In this study, where λ is the confidence degree, which is generally taken between [0.5-0.7].In this study, λ was taken as 0.6.If f 1 > λ, the risk of the cumulative impact of the heritage site was level I.If f 1 + f 2 > λ, it was level II.If f 1 + f 2 + f 3 > λ, the risk was level III; and so on.The risk level of each heritage site can thus be calculated.

Study area
Suzhou is located in the Jiangsu Province of China and has a history of over 2500 years [57].It is the prefecturelevel city in China with the highest economic development and one of the first 24 historical and cultural cities in China.As one of the most important cultural heritages of Suzhou, Suzhou Classical Gardens have a long history that extends back to the Spring and Autumn Period in the sixth century BC.Suzhou Classical Gardens are gardens built in ancient times or inherited from the typical local gardening techniques in Suzhou, China.They contain various landscape elements and complex spaces, are the essence of Chinese gardening art and valuable world cultural heritage sites [58].Nine of the Suzhou Classical Gardens are on the World Heritage List designated by UNESCO.However, the urbanization process of Suzhou has threatened the protection of the classical gardens.The urbanization outbreak in Suzhou is strongly linked to the Chinese strategy for reform and opening up in 1978.Rapid economic development has promoted the process (17) w m c m of urbanization.Additionally, their recognition as World Heritage Sites by UNESCO in 1997 has caused a wave of tourists to the city of Suzhou.Suzhou has had to reshape and develop its urban infrastructure to meet the requirements of the tourism industry.Therefore, it is necessary to assess the impact of urbanization on the Suzhou Classical Gardens.
Gusu District belongs to Suzhou City and is located in the central part of the city (Fig. 2).With more than 4000 years since the first written records were published, Gusu District is the oldest central city of Suzhou and has the first National Historical and Cultural City Reserve in China.According to the statistics on the official website of the Suzhou City Landscape and Greening Administration, more than half of Suzhou Classical Gardens are located in the Gusu District.
This study focused on 21 Suzhou Classical Gardens in Gusu District (Fig. 3), analyzing the impact of rapid population growth and economic expansion in recent years.The name and protection level of the 21 Suzhou classical gardens can be seen in Table 3.The aim was to provide theoretical support for cultural heritage protection in Suzhou and a reference for analyzing urbanization effects on cultural heritage sites in other cities.

Data sources
The data used in this study included urban planning, heritage statistics, social and economic data, and street map data.
2) Heritage statistical data: Statistical data on Suzhou Classical Gardens were obtained from the official website of the Suzhou City Landscape and Greening Administration (http:// ylj.suzhou.gov.cn/ szsylj/ ylml/ nav_ list.shtml).
3) Social and economic data: A statistical report on the national economic and social development of Suzhou's subordinate districts and counties in 2022 was used.4) Street map data: We collected road data from the Open Street Map (https:// www.opens treet map.org/).Street map data were used to calculate garden accessibility.

Weight analysis
The subjective weights and the objective weights of each index were obtained following the AHP and EW methods, respectively.Subsequently, the combination weights were obtained based on game theory.Table 4 lists the obtained weights.As shown in Table 3, D2, D4, and D9 were significant.This result shows that construction and commercial activities have a critical cumulative impact on the cultural heritage sites.Additionally, D6, D7, and D8 played a key role in the cumulative impact assessment.This demonstrates that transportation activities, such as subways and heavy traffic, often exceed limits and potentially damage the heritage sites.However, other indexes, such as D1, D5, D10, and D11, accounted for small weights.These results demonstrate that the cumulative impact is influenced more by the external urban environment of the heritage site, whereas indexes related to information about the heritage site itself are less important.

Cumulative impact assessment of Suzhou Classical Gardens in the Gusu District
Based on Eqs.11 and 12, the single correlation degree of each index of Suzhou Classical Gardens can be calculated.Taking the Master of Nets Garden as an example, the single correlation degree for each index of The Master of Nets Garden is listed in Table 5.The integrated degree of correlation can be obtained using Eqs.14, 15, 16, 17, 18 through each index's weight: f 1 = 0.1378, f 2 = 0.2667, f 3 = 0.4786, f 4 = 0.0398, f 5 = 0.0771, f 1 + f 2 + f 3 > 0.6.Therefore, the risk of the cumulative impact of the Master of Nets Garden can be determined as Level III.
By repeating the above process, the integrated degree of correlation and risk levels of the 21 Suzhou Classical Gardens were obtained (Table 6).The results are listed in Table 7. Fifty-seven percent of the classical gardens have cumulative impact risk levels at or below "medium" risk.Overall, five classical gardens were classified as having a cumulative impact risk level of "highest, " four as having a cumulative impact risk level of "higher, " nine as having a "medium" level, two as having a "lower" level, and one was classified as having a cumulative impact risk level of "lowest." The assessment results confirmed that the urbanization process in Suzhou City has had an impact on the Suzhou Classical Gardens because 85.72% of the gardens have a cumulative impact risk level at or higher than the medium risk level.Moreover, among the nine high-risk Suzhou Classical Gardens, seven are located at the edge of the protection boundary of the historical and cultural city, and only two are located in the center (Fig. 4).This finding underscores the significance of historical and cultural cities in preserving Suzhou Classical Gardens.

Discussion
This study developed a risk-based CIA framework that incorporates a set of indexes to assess the impact of urbanization on urban cultural heritage.In contrast to the other CIA method, this study built a risk-based CIA method that relies on describing the probability of adverse impact, which is useful for managing excessive scientific complexity and uncertainties.Different from other CIA methods, future adverse impacts were considered in this study and index like Areas of new construction land (D2) and percentage of surrounding commercial land (D3) were added in this study.Moreover, unlike other research, game theory was introduced in this study to obtain a more reasonable weight of each index.
Although the 21 gardens selected in this study are all located in the Suzhou Historical and Cultural City Protected Area, the results of the study show that Suzhou's urbanization has inevitably affected these historic gardens.According to the combination weight obtained by game theory (Table 4), the area of new construction land, the area of the cultural heritage and the harmony index of neighboring building form are three of the most important indexes in this study.All three indexes are related to the surrounding urban environment.Therefore, urbanization impact on 21 gardens is probably due to the increase in construction projects, subway construction, and commercial land use in the vicinity of the historic gardens as a result of the city's tourism development in recent years.Therefore, more restrictive management of buffer zone construction are needed for Suzhou Classical Gardens to guarantee their heritage integrity.Meanwhile, among the nine high-risk Suzhou Classical Gardens, seven are located at the edge of the protection boundary of the historical and cultural city.This demonstrates that the establishment of Suzhou National Historical and Cultural City has a positive impact on the protection of Suzhou Classical Gardens.Moreover, six out of the nine high-risk gardens are below the national-level-protected, which indicates  that the current protection measures for gardens with a low level of protection are weak.While the proposed risk-based cumulative impact assessment method yields clear results on the impact of urbanization on cultural heritage, some limitations must be clarified.First, this method was developed using a reduced set of indexes, and a more thorough assessment result can be obtained by including more indices and data.This method aimed to provide a standardized riskbased CIA method for cultural heritage sites that can be modified and improved by including additional indicators.Second, the study focused on one type of cultural heritage, the Suzhou Classical Gardens.The assessment index system should be modified when applied to other types of cultural heritage based on their characteristics.When adopting the suggested approach to analyze the cumulative impact of urbanization on urban cultural heritage, it is vital to account for these limitations.Third, urbanization is a complex activity that can be affected by social and economic factors.Building future risk scenarios for all anthropogenic activities and pressures that could impact cultural heritage sites is a significant challenge.This index-based CIA method is a first-level analysis of cultural heritage sites.To obtain more accurate assessment results, the next stage of the risk-based CIA

Conclusions
This study proposed a risk-based CIA method to evaluate the impact of urbanization on urban cultural heritage sites at a regional level.This method is based upon the SPA theory and uses 11 indexes as assessment factors, and the combined weight of the indexes was rationally obtained using the game theory.The effectiveness of this method was demonstrated through its successful implementation in 21 classical gardens in Suzhou, China, to assess the impact of urbanization on urban cultural heritage sites.The assessment results show that the Suzhou Classical Gardens, which are important cultural heritage sites in China, are inevitably affected by urbanization.This is because 85.72% of the case studies have a cumulative impact risk level at or higher than the "medium" risk level.The proposed method is based on objective and unbiased data and can provide standardized outcomes for decision-makers, which allows them to prioritize cultural heritage sites with a higher cumulative impact.
The CIA method aims to provide an overall understanding of the conservation status of cultural heritage sites, identify potential sources of impact such as development projects or transportation infrastructure, and then provide conservation agencies with a picture of the impact of the current value of cultural heritage sites and guidance by proposing targeted conservation strategies.The CIA method proposed can help fill the gap between city development demands and heritage protection by facilitating the determination and efficient implementation of mitigation strategies in the decision-making process.Meanwhile, the universality of the suggested risk-based CIA method is beneficial for cultural heritage protection, as input datasets are easily obtainable and analysis processes are reproducible for other kinds of cultural heritage sites.Furthermore, the index system can be adapted to different regions and different types of cultural heritage.Therefore, we argue that the risk-based CIA method can be a strategic approach used before the urban planning process, making it easier to communicate key aspects of heritage protection recommendations to decision makers.

Fig. 1
Fig. 1 Flowchart of the risk-based CIA method employed in this study from subway station (D7) Negative Distance from subway tracks (D8) Negative Cumulative vulnerability Harmony index of neighboring building form (D9) Negative Heritage site history (D10) Positive Heritage protection level (D11) Positive

Fig. 2
Fig. 2 Photographs that depict the urban environment of Gusu District; a, b: Photograph of the contemporary style buildings; c: Photograph of the Suzhou classical style buildings

Fig. 3
Fig. 3 Distribution of 21 Suzhou classical gardens in the Gusu District

Fig. 4
Fig. 4 Spatial distribution of the cumulative impact level of Suzhou Classical Gardens in the Gusu district

Table 2
Grading standard of cumulative impact assessment indexes

Table 3
Name and protection level of the 21 Suzhou Classical Gardens

Table 4
Weight of each index obtained by the combination weighting method

Table 5
Single correlation degree of each index of the Master of Nets Garden

Table 6
Integrated correlation degree and risk levels of the cumulative impact of 21 Suzhou Classical Gardens

Table 7
Risk levels of 21 Suzhou Classical Gardens