Impacts of Climate Change and Land Subsidence on Inundation Risk

In this study, a physiographic drainage-inundation model was used to analyse the impacts of land subsidence and climate change on inundation disaster and risk in a land subsidence area. The results indicated that for land subsidence and land subsidence combined with climate change, inundation area, and volume increased under oneand two-day storms for 2-, 5-, 10-, 25-, 50-, 100-, and 200-year return periods. Moreover, locations that originally had high inundation depth showed even greater inundated areas and volumes in the presence of land subsidence. The inundation phenomenon under the combination of land subsidence and climate change proved to be severe, compared to that of land subsidence alone. Land subsidence increased not only inundation depth but also inundation duration. Given land subsidence and climate change, the average inundation duration for each return period increased. The average flooding duration for each return period post land subsidence was found to be 1.05–1.1 times greater than that preceding it. Under the combination of land subsidence and climate change, the average flooding duration for each return period post land subsidence was about 1.13–1.27 times greater than that before it. Furthermore, by assessing inundation risk with inundation depth index, inundation duration index, and damage index from different land uses, it was found that after land subsidence, inundation risk showed an increase, which was amplified in the presence of land subsidence combined with climate change.


Introduction
Owing to an inadequate supply of surface water or economic considerations, people in certain regions of Taiwan have extracted considerable groundwater, causing serious problems, such as land subsidence, seawater encroachment, and soil salinization.According to a Taiwan Water Resources Agency report [1], land subsidence, caused by excessive extraction of groundwater is a concern for approximately one-tenth of Taiwan's plains.In addition, most of the land subsidence areas are located in coastal plains.When such coastal plains experience land subsidence, flooding and seawater encroachment complicate water drainage, leading to severe inundation disasters and jeopardizing national safety.
On average, Taiwan is struck by 3-4 typhoons annually, which lead to seawater intrusion, difficulty in draining flow discharge in rivers, and increased water levels, eventually resulting in flood disasters.According to [2], the maximum inundated depth, area, and time increase with precipitation volume and tidal surge level, leading to severe inundation disasters.Recently, global warming and climate change have intensified, and, according to a report from the Intergovernmental Panel on Climate Change (IPCC) [3], global temperatures are continually rising.Such rising temperatures Water 2018, 10, 157 2 of 21 must inevitably lead to abnormal global climate and extreme hydrological phenomena, which will subsequently result in the frequent arrival of natural disasters, such as floods and droughts.Taiwan is situated in a region that is highly vulnerable to climate change.According to statistical data [4], climate change has been increasing the frequency and volume of precipitation and has caused the sea level to rise.For instance, the average rate of sea-level rise for oceans near Taiwan has reached 2.4 mm/year [5], which is approximately 1.3 times higher than the corresponding rate worldwide (1.8 mm/year).In Taiwan's Southwestern coastal areas, the rate has surged to a maximum of 7.89 mm/year.Therefore, it is apparent that the combined effects of climate change and land subsidence are among the greatest environmental hazards faced by Taiwan, leading to the worst ever inundation and sea-level rise, causing serious harm to the sustainable development of the economy and society [4].
Situations with similar inundation potentials do not necessarily generate the same inundation disaster.The occurrence of a disaster and its degree are typically influenced by the target, time, space of concern, and land utilization conditions.For example, several inundations of identical depth, occurring at different times, locations, and land utilization conditions, cause disasters of different severities.Uitto [6] proposed that the occurrence and severity of a disaster is determined by three basic variables: (1) hazard (natural phenomena such as earthquakes and typhoons), ( 2) exposure (structures, buildings, humans, and other entities at risk), and (3) vulnerability (propensity to suffer loss).Among these, the concept of disaster risk is the concept that the occurrence of a natural disaster involves a certain degree of uncertainty.As a result, the risk concept is introduced for quantifying disaster loss-and disaster management-related studies [7].The degree or severity of a disaster depends on the target, the environment exposed to the disaster, and the amount of loss induced by the disaster.Therefore, researchers define disaster risk differently, depending on different research conditions.In the field of natural disaster risk assessment, several assessment equations have been established.Factors such as probability, disaster damage (consequence), vulnerability, and hazard have all been considered by these equations.For example, Hammer [8] defined disaster risk as (1) a combination of probability for the disaster to occur and the damage caused by the disaster (Risk = Probability × Consequence); (2) a combination of probability for the disaster to occur and the number of entities exposed to the area affected by the disaster (Risk = Probability × Exposure); (3) Risk = Probability × Vulnerability; and finally, (4) Risk = Probability × Exposure × Vulnerability.
The severity of a flooding disaster is related to the number of people and entities exposed to the disaster as well as the resulting loss.Therefore, how to quantify the damage caused by an inundation disaster is an important factor for planning disaster preventive measures, adjustment strategies, and disaster risk management.Inundation risk assessment is broadly categorized into two types: (1) risk = flood hazard × fragility (exposure × social and regional fragility) (e.g., Satu [9][10][11]) and (2) flood risk = flood occurrence probability × disaster damage (e.g., [11,12]).For the inundation hazard factor, besides inundation depth, the probability of inundation is considered.Hashimoto, Stedinger and Loucks [7] defined inundation probability by accounting for the duration of inundation, which means that inundation probability is the ratio of inundation duration during a rain incident to the total duration of the rain incident.Vulnerability is categorized according to the land use of the studied area, and is allocated a score based on the monetary value.Finally, the degree of damage is used to determine the vulnerability score.
This study aimed to understand the impacts of land subsidence and climate change on inundation disaster and risk.In addition, an inundation risk assessment method was established.The serious land subsidence area at the Linbian River Basin was selected as the study area.A physiographic inundation-drainage model was applied and combined with a geographic information system (GIS), to simulate the potential inundation, inundated depth, inundated area, inundated volume, and inundation duration under a normal hydrological situation or under a hydrological situation changed by land subsidence and climate change.Further, inundation depth, duration, and loss (vulnerability) were considered, to conduct an inundation risk analysis and to discuss the impacts of land subsidence and climate change on inundation disaster and risk.

Study Area
The study area was Linbian River Basin, located in the Southern part of Taiwan.It has a surface area of 336.3 km 2 and an approximately 41.3 km long mainstream.The mainstream has an average slope of 1/88 and an average annual precipitation of approximately 2115 mm.Linbian River, whose map is shown in Figure 1, is one of the major rivers in Pingtung County.

Study Area
The study area was Linbian River Basin, located in the Southern part of Taiwan.It has a surface area of 336.3 km 2 and an approximately 41.3 km long mainstream.The mainstream has an average slope of 1/88 and an average annual precipitation of approximately 2115 mm.Linbian River, whose map is shown in Figure 1, is one of the major rivers in Pingtung County

Method and Data
Figure 2 shows the flowchart of inundation risk analysis in this paper.The Physiographic Drainage-Inundation model was used to simulate inundation depth, inundation duration, inundation area and volume for various scenarios, to calculate the impact of climate change and land subsidence on inundation risk.

Physiographic Drainage-Inundation (PDI) Model
Based on the topography, landform and river system of the catchment area, the PDI model utilized GIS to partition the catchment area into either land or river computational cells.Hydrological and physiographical data within each cell were obtained through spatial and 3D analyst tools embedded in ESRI ArcMap 10.3.The PDI model then calculated the inundation potential according to the partitioned cell and its attributes.In the migration process of water for the flooding zone, water level change of each cell should then satisfy the continuity equation, as expressed in (1) [13]: where t is time; Asi is the area of the i cell; hi and hk represent water stages of the i and k cells, at time t, respectively; and Qi,k denotes the discharge from the k cell into its neighbouring i cell.Discharge is positive when flowing into the i cell from its adjacent k cell and is negative when flowing out of the i cell.Pe,i expresses the effective rainfall volume per second in the i cell, which is equal to the effective rainfall per second in the i cell, multiplied by its area.

Method and Data
Figure 2 shows the flowchart of inundation risk analysis in this paper.The Physiographic Drainage-Inundation model was used to simulate inundation depth, inundation duration, inundation area and volume for various scenarios, to calculate the impact of climate change and land subsidence on inundation risk.

Physiographic Drainage-Inundation (PDI) Model
Based on the topography, landform and river system of the catchment area, the PDI model utilized GIS to partition the catchment area into either land or river computational cells.Hydrological and physiographical data within each cell were obtained through spatial and 3D analyst tools embedded in ESRI ArcMap 10.3.The PDI model then calculated the inundation potential according to the partitioned cell and its attributes.In the migration process of water for the flooding zone, water level change of each cell should then satisfy the continuity equation, as expressed in Equation (1) [13]: where t is time; A si is the area of the i cell; h i and h k represent water stages of the i and k cells, at time t, respectively; and Q i,k denotes the discharge from the k cell into its neighbouring i cell.Discharge is positive when flowing into the i cell from its adjacent k cell and is negative when flowing out of the i cell.P e,i expresses the effective rainfall volume per second in the i cell, which is equal to the effective rainfall per second in the i cell, multiplied by its area.

Hazard
Exposure Vulnerability

Monetary value investigated
Vulnerability score index  Flow discharge between two adjacent cells is based on the connections between these cells, and can be flow type and box culvert flow type, as described below:

Inundation loss analysis
(1) River flow type If no obvious obstacle exists between neighbouring cells where the flow take place, or there is no flow exchange between neighbouring cells, then it is regarded as overland flow and the average resistance equation (i.e., Manning formula) is used to calculate the flow from a cell to its neighbouring cell.Accordingly, the flow of water from the k cell into its neighbouring i cell is denoted as follows: where ∆x is the distance between the centre of the i and k cells; n is the Manning roughness coefficient of overland flow between the two neighbouring cells; and A and R are the hydraulic area and radius, respectively, at the border between the two neighbouring cells.
(2) Weir flow type If the areas are divided by hydraulic or artificial structures, such as roadways, embankments, field ridges, or banks, then the border may be treated as broad-crested weir, and the weir flow formula is used to obtain flow from one cell to its neighbouring cell.Such flow exchange between cells is regarded as the weir flow type.If h k > h i , then there are two possible cases-the free weir and the submerged weir-as described below: a.
Free weir where Z w is the weir height, which is the roadway, bank or ground height; b is the effective width of the weir top, which is equivalent to the intersection length of two adjacent cells; g is gravity; µ f and µ s are the weir coefficients of the free and submerged weirs, respectively.
If the flooding area were partitioned into N cells, then there will be N differential equations as Equation ( 1) with N unknown cell water stages.The correlation between discharge from the k cell into its neighbouring i cell, Q i,k , and water stages, h i and h k , for cells i and k may be derived from Equations ( 2)-( 4) based on the cell connection type.Using the explicit finite difference method and taking i cell as an example, Equation (1) can be expressed as: where h m+1 i represents the water stage in the i cell at time m + 1; h m i represents the water stage in the i cell at time m; ∆t signifies the time increment between m and m + 1.The correlation between discharge from the k cell into its neighbouring i cell, Q i,k , and water stages, h i and h k , for cells, i and k, may be derived from Equations ( 2)-( 4).The water stage hydrographs can be obtained from Equation (5) and the corresponding hydrological and physiographic data.

Land Subsidence Depth
Linbian River Basin is regarded as an area with serious land subsidence (photos in Figure 3).Data were accumulated to quantify the land subsidence of previous years; estimates were calculated based on the land subsidence data collected at the land subsidence inspection points, set up by the Water Resources Agency, Ministry of Economic Planning and Agriculture (MOEA), in the Linbian River catchment area and in nearby areas.The amount of land subsidence accumulated for each inspection point was used together with GIS to construct the accumulated land subsidence numerical diagram shown in Figure 3.
By overlapping the computed cells and the accumulated land subsidence numerical diagram, the amount of land subsidence in each cell can be determined.By subtracting this amount from the original topographical height of each cell, the height of each computed cell after land subsidence can be obtained.This height can then be used to conduct inundation simulation analysis.
Water 2018, 10, x FOR PEER REVIEW 6 of 21

Land Subsidence Depth
Linbian River Basin is regarded as an area with serious land subsidence (photos in Figure 3).Data were accumulated to quantify the land subsidence of previous years; estimates were calculated based on the land subsidence data collected at the land subsidence inspection points, set up by the Water Resources Agency, Ministry of Economic Planning and Agriculture (MOEA), in the Linbian River catchment area and in nearby areas.The amount of land subsidence accumulated for each inspection point was used together with GIS to construct the accumulated land subsidence numerical diagram shown in Figure 3.
By overlapping the computed cells and the accumulated land subsidence numerical diagram, the amount of land subsidence in each cell can be determined.By subtracting this amount from the original topographical height of each cell, the height of each computed cell after land subsidence can be obtained.This height can then be used to conduct inundation simulation analysis.

Inundation Risk Assessment
In this study, inundation risk assessment was performed using the three variables proposed by [6], which affect the occurrence and degree of a disaster: (1) hazard; (2) exposure; and (3) vulnerability.For instance, the concept of Risk = Hazard × Exposure × Vulnerability was applied, where the exposure variable not only accounted for the effect of land utilization but also that of inundation duration.The method used for inundation risk assessment is described as follows: (1) Inundation depth score index The degree of damage caused by flooding depends on inundation depth.According to the results of the inundation depth and monetary loss curve for different land uses, inundation depths induce different degrees of damage-the higher the inundation depth, the larger the inundation damage.For this reason, the inundation risk assessment performed in this study considered the effects of different inundation depths.The inundation depths were categorized into five levels, with each level assigned a score index, as shown in Table 1.
(2) Inundation duration score index Besides inundation depth, inundation damage is also dependent on flooding duration-the longer the flooding lasts, the greater the damage.Therefore, inundation duration was also

Inundation Risk Assessment
In this study, inundation risk assessment was performed using the three variables proposed by [6], which affect the occurrence and degree of a disaster: (1) hazard; (2) exposure; and (3) vulnerability.For instance, the concept of Risk = Hazard × Exposure × Vulnerability was applied, where the exposure variable not only accounted for the effect of land utilization but also that of inundation duration.The method used for inundation risk assessment is described as follows: (1) Inundation depth score index The degree of damage caused by flooding depends on inundation depth.According to the results of the inundation depth and monetary loss curve for different land uses, inundation depths induce different degrees of damage-the higher the inundation depth, the larger the inundation damage.For this reason, the inundation risk assessment performed in this study considered the effects of different inundation depths.The inundation depths were categorized into five levels, with each level assigned a score index, as shown in Table 1.
(2) Inundation duration score index Besides inundation depth, inundation damage is also dependent on flooding durationthe longer the flooding lasts, the greater the damage.Therefore, inundation duration was also Water 2018, 10, 157 7 of 21 considered while performing the inundation risk assessment in this study.Here, D cr represents the critical inundation depth in cell i and was used to determine whether cell i was undergoing flooding.For example, D i ≥ D cr indicates inundation.Following the guidelines of the Taiwan Water Resources Agency [1], we set D cr as 30 cm.If the rainfall duration is r m hour and the inundation duration is r d hour in cell i, the inundation probability is r d /r m .Based on the level of flooding probability, an inundation duration score index was assigned, as shown in Table 1.(3) Vulnerability score index The study area covers a wide range of land uses (Figure 4), such as buildings, agricultural land, forest land, and water resources (rivers and drainage); the flooding zone of the study area constitutes low-altitude coastal plains, and the coastal area is mainly used for the construction of populated buildings and agricultural land.Agricultural land includes rice farms, rain-fed crop farms, fruit farms, and fish farms.The agricultural land of the investigated area in this study is mostly occupied by fish farms, at which high-value species, such as grouper and shrimp, particularly the Epinephelus lanceolatus (common name: Giant grouper, Brindle bass, and Queensland grouper) are cultivated.For this reason, the vulnerability index of fish farming was separated from that of normal agricultural land, and was assigned according to the monetary values shown in Table 3.
Water 2018, 10, x FOR PEER REVIEW 7 of 21 considered while performing the inundation risk assessment in this study.Here, Dcr represents the critical inundation depth in cell i and was used to determine whether cell i was undergoing flooding.For example, Di ≥ Dcr indicates inundation.Following the guidelines of the Taiwan Water Resources Agency [1], we set Dcr as 30 cm.If the rainfall duration is rm hour and the inundation duration is rd hour in cell i, the inundation probability is rd/rm.Based on the level of flooding probability, an inundation duration score index was assigned, as shown in Table 1.(3) Vulnerability score index The study area covers a wide range of land uses (Figure 4), such as buildings, agricultural land, forest land, and water resources (rivers and drainage); the flooding zone of the study area constitutes lowaltitude coastal plains, and the coastal area is mainly used for the construction of populated buildings and agricultural land.Agricultural land includes rice farms, rain-fed crop farms, fruit farms, and fish farms.Table 2 shows the monetary value of each land use according to the investigation conducted by the Sixth River Bureau of the Water Resources Agency [1].The monetary value of output for various land use per unit area and the consumer price index [14] of Taiwan were used to compute the monetary values presented in Table 2.According to the table, construction land or buildings have the highest monetary values and fishing farms have the highest monetary values among other agricultural activities.Land uses influence not only monetary values, but also the damage induced by disasters, even multiple floods have the same inundation depth.Hence, classification of the investigated area is based on the different land uses and vulnerability score indices, as determined by the monetary values of these land uses.The agricultural land of the investigated area in this study is mostly occupied by fish farms, at which high-value species, such as grouper and shrimp, particularly the Epinephelus lanceolatus (common name: Giant grouper, Brindle bass, and Queensland grouper) are cultivated.For this reason, the vulnerability index of fish farming was separated from that of normal agricultural land, and was assigned according to the monetary values shown in Table 3.    (4) Inundation risk Inundation risk represents the possibility of human injury or death and the loss of property due to flooding.Inundation damage is correlated with inundation depth, inundation duration, and land use of the inundated area.Therefore, the inundation risk assessment conducted in this study mainly considered inundation depth, inundation duration, and the inundation loss of land areas with different uses in the flooding area.This index was calculated by multiplying the inundation depth index and inundation probability index (inundation duration) with the inundation disaster loss index, and then dividing it by the product of the highest loss index, as shown below: Inundation Risk (%) =

Hazard×Exposure×Vulnerability
The product of the highest score index × 100 =

Depth index×Duration index×Loss index
The product of the highest score index × 100 (6)

Construction of the Computed Cells
Water runoff from rainfall is influenced by catchment topography, landform, vegetation distribution, and land use.Therefore, to simulate the water runoff process for the investigated catchment area, we partitioned the area into several computed cells based on the aforementioned factors.The complex topography and landforms generated considerable hydrological and physiographical data, and to handle these data ESRI ArcMap 10.3 was used, along with specific spatial analyst, hydrology, 3D analyst, and network analyst extension modules.In addition, object-oriented programming language, flow direction analysis, and formation of river lines were used to automatically partition the catchment area into several sub-catchment areas.
The digital elevation model (DEM) forms the basis of the automatic partitioning of the sub-catchment area.The adopted DEM (Figure 5) had a spatial resolution of 5 m × 5 m in a regular cell.River width and cells for the automatically partitioned catchment area were assigned based on the recent river measurement data provided by the Linbian River Bureau.To aid in catchment partitioning, the DEM was overlaid with additional maps, such as digitalized soil, land use, and road system maps.Extremely small or narrow cells were combined, and large cells were divided into appropriate cells.Finally, 4290 cells were defined; they are shown in Figure 5. Attributes of each cell served as input to the PDI model.

Hydrological Data of Linbian River Basin
Downscaled monthly precipitation data of climate change scenarios, at a resolution 5 km × 5 km grid, was provided by the Taiwan Climate Change Projection and Information Platform Project (TCCIP), which is coordinated by the National Science and Technology Centre for Disaster Reduction (NCDR).TCCIP's processing is based on 24 general circulation models (GCMs) [16], reported by the Fourth Assessment Report (AR4) of Intergovernmental Panel on Climate Change (IPCC) [3].According to the AR4 of IPCC, A2, A1B, and B1 are the major scenarios.Among these, A1B is believed to be the most probable scenario in the future and is expected to receive the most widespread attention.Henceforth, this study adopted the "worst-case scenario" of A1B-S, which is similar to the A1B scenario to evaluate the potential impacts of climate change.The worst-case scenario is mainly determined by adding or subtracting one standard deviation between the estimated values of GCMs from the multi-model ensemble of all GCMs.The future period in this study was set to 2020-2039 and the baseline period set to 1980-1999.The precipitation data of baseline and A1B-S for one and two-day storms under various return periods are shown in Table 4.

Hydrological Data of Linbian River Basin
Downscaled monthly precipitation data of climate change scenarios, at a resolution 5 km × 5 km grid, was provided by the Taiwan Climate Change Projection and Information Platform Project (TCCIP), which is coordinated by the National Science and Technology Centre for Disaster Reduction (NCDR).TCCIP's processing is based on 24 general circulation models (GCMs) [16], reported by the Fourth Assessment Report (AR4) of Intergovernmental Panel on Climate Change (IPCC) [3].According to the AR4 of IPCC, A2, A1B, and B1 are the major scenarios.Among these, A1B is believed to be the most probable scenario in the future and is expected to receive the most widespread attention.Henceforth, this study adopted the "worst-case scenario" of A1B-S, which is similar to the A1B scenario to evaluate the potential impacts of climate change.The worst-case scenario is mainly determined by adding or subtracting one standard deviation between the estimated values of GCMs from the multi-model ensemble of all GCMs.The future period in this study was set to 2020-2039 and the baseline period set to 1980-1999.The precipitation data of baseline and A1B-S for one and two-day storms under various return periods are shown in Table 4.The downstream flooding phenomenon is closely related to tidal surge levels.Therefore, in this study, the tidal levels at the Dong Gang Tidal Monitoring Station in 2011 (measured by Tainan Hydraulic Laboratory, National Cheng Kung University) and the tidal levels near Linbian River were considered in the tidal level analysis.These data were also used to determine the downstream boundary conditions in our simulations.The existing tidal level and the increased sea level were used to determine the tidal levels in the climate change scenario.

Model Verification
Typhoon Morakot (2009), which caused severe flooding disaster, was selected to validate the PDI model.We simulated flooding for Typhoon Morakot, from 6 to 10 August 2009, for a total of 120 h.Precipitation data were collected from three rainfall stations, and Thiessen's polygon method was used to partition the control area of each rainfall monitoring station, to serve as an upstream boundary input for the PDI model.The downstream boundary conditions for the simulation were based on data from the Dong Gang Tidal Level Monitoring Station.
The investigated and simulated inundation areas for Typhoon Morakot are shown in Figure 6.The maximum possible inundation depth and inundation area simulated by the PDI model agreed well with the inundation determined on-site.
Water 2018, 10, x FOR PEER REVIEW 10 of 21 The downstream flooding phenomenon is closely related to tidal surge levels.Therefore, in this study, the tidal levels at the Dong Gang Tidal Monitoring Station in 2011 (measured by Tainan Hydraulic Laboratory, National Cheng Kung University) and the tidal levels near Linbian River were considered in the tidal level analysis.These data were also used to determine the downstream boundary conditions in our simulations.The existing tidal level and the increased sea level were used to determine the tidal levels in the climate change scenario.

Model Verification
Typhoon Morakot (2009), which caused severe flooding disaster, was selected to validate the PDI model.We simulated flooding for Typhoon Morakot, from 6 to 10 August 2009, for a total of 120 h.Precipitation data were collected from three rainfall stations, and Thiessen's polygon method was used to partition the control area of each rainfall monitoring station, to serve as an upstream boundary input for the PDI model.The downstream boundary conditions for the simulation were based on data from the Dong Gang Tidal Level Monitoring Station.
The investigated and simulated inundation areas for Typhoon Morakot are shown in Figure 6.The maximum possible inundation depth and inundation area simulated by the PDI model agreed well with the inundation determined on-site.

Inundation Disaster Analysis
To understand the effects of land subsidence and climate change on inundation disaster in the land subsidence area, the PDI model was adopted to simulate the impact of land subsidence and land subsidence combined with climate change on the inundation disaster with rainfall durations of 24 h and 48 h, under the hydrological condition of normal (baseline) and climate change (A1B-S) scenarios, for the return periods of 2, 5, 25, 50, 100, and 200 years.One-and two-day storms for the 5-and 50year return periods were selected to demonstrate the increase in inundation depth and distribution after land subsidence and in land subsidence area, under the impact of climate change, as shown in

Inundation Disaster Analysis
To understand the effects of land subsidence and climate change on inundation disaster in the land subsidence area, the PDI model was adopted to simulate the impact of land subsidence and land subsidence combined with climate change on the inundation disaster with rainfall durations of 24 h and 48 h, under the hydrological condition of normal (baseline) and climate change (A1B-S) scenarios, for the return periods of 2, 5, 25, 50, 100, and 200 years.One-and two-day storms for the 5-and 50-year return periods were selected to demonstrate the increase in inundation depth and distribution after land subsidence and in land subsidence area, under the impact of climate change, as shown in Figures 7-10.Furthermore, the percentages of inundation area and volume increase under different inundation depths in land subsidence and climate change with land subsidence areas are listed in Tables 5-8.5-8 indicate that under the influence of land subsidence and land subsidence combined with climate change, the percentage of inundation area and volume with oneand two-day storms, for different inundation depths in each return period, may increase or decrease.Nevertheless, the overall results showed greater percentages of increase for inundation area and volume.Furthermore, the places with greater inundation depths originally tended to show increased inundation areas and volumes after land subsidence.

Inundation Risk Analysis
By applying inundation depth and duration derived from the inundation potential analysis and combining the land uses as well as the corresponding monetary value of the inundation area, the inundation loss index was assigned.This index can be used to calculate the inundation vulnerability for different return periods.Finally, inundation risk was determined.By subtracting the inundation risk before and after land subsidence, maps of the increase and distributions of inundation risk were prepared.The one-and two-day storms for the return periods of 5 years and 50 years were selected as examples, to demonstrate the influence of land subsidence on inundation risk, as shown in Figures

Inundation Risk Analysis
By applying inundation depth and duration derived from the inundation potential analysis and combining the land uses as well as the corresponding monetary value of the inundation area, the inundation loss index was assigned.This index can be used to calculate the inundation vulnerability for different return periods.Finally, inundation risk was determined.By subtracting the inundation risk before and after land subsidence, maps of the increase and distributions of inundation risk were prepared.The one-and two-day storms for the return periods of 5 years and 50 years were selected as examples, to demonstrate the influence of land subsidence on inundation risk, as shown in Figures

Inundation Risk Analysis
By applying inundation depth and duration derived from the inundation potential analysis and combining the land uses as well as the corresponding monetary value of the inundation area, the inundation loss index was assigned.This index can be used to calculate the inundation vulnerability for different return periods.Finally, inundation risk was determined.By subtracting the inundation risk before and after land subsidence, maps of the increase and distributions of inundation risk were prepared.The one-and two-day storms for the return periods of 5 years and 50 years were selected as examples, to demonstrate the influence of land subsidence on inundation risk, as shown in Figures 13  and 14.These figures show that locations such as the region between the railway and provincial road No. 17, as well as the boundary between the coastal bank and catchment area had higher increases in inundation risk since they were the most affected by land subsidence.Therefore, land subsidence worsens inundation hazard and disaster loss.
13 and 14.These figures show that locations such as the region between the railway and provincial road No. 17, as well as the boundary between the coastal bank and catchment area had higher increases in inundation risk since they were the most affected by land subsidence.Therefore, land subsidence worsens inundation hazard and disaster loss.13 and 14.These figures show that locations such as the region between the railway and provincial road No. 17, as well as the boundary between the coastal bank and catchment area had higher increases in inundation risk since they were the most affected by land subsidence.Therefore, land subsidence worsens inundation hazard and disaster loss.Due to the influence of global climate change in recent years, abnormal weather conditions, such as extreme rainfall, have been observed more frequently in Taiwan.To understand the impacts of climate change and land subsidence on inundation risk in the study area, separate scenarios of inundation depth and duration for land subsidence with one-and two-day storms for each return period in the A1B-S scenario were simulated.By further applying the land use as well as the corresponding monetary value of the inundation area, the inundation loss index was assigned, which can be used to calculate the inundation vulnerability for different return periods and further determine the inundation risk.By subtracting the determined inundation risk of the A1B-S scenario from the current inundation risk, the increase and distribution maps of inundation risk were prepared.The one-and two-day storms for the return periods of 5 years and 50 years were selected as examples, to demonstrate the influence of climate change on inundation risk, as shown in Figures 15 and 16.These figures show that the inundation risk with one-and two-day storms under different return periods in the A1B-S scenario exhibited an increasing trend, and the increase in inundation risk was significantly higher than that with land subsidence alone.This implies that the inundation risk for the A1B-S scenario is higher than that of the existing condition, particularly for a storm with a larger return period.Due to the influence of global climate change in recent years, abnormal weather conditions, such as extreme rainfall, have been observed more frequently in Taiwan.To understand the impacts of climate change and land subsidence on inundation risk in the study area, separate scenarios of inundation depth and duration for land subsidence with one-and two-day storms for each return period in the A1B-S scenario were simulated.By further applying the land use as well as the corresponding monetary value of the inundation area, the inundation loss index was assigned, which can be used to calculate the inundation vulnerability for different return periods and further determine the inundation risk.By subtracting the determined inundation risk of the A1B-S scenario from the current inundation risk, the increase and distribution maps of inundation risk were prepared.The one-and two-day storms for the return periods of 5 years and 50 years were selected as examples, to demonstrate the influence of climate change on inundation risk, as shown in Figures 15 and 16.These figures show that the inundation risk with one-and two-day storms under different return periods in the A1B-S scenario exhibited an increasing trend, and the increase in inundation risk was significantly higher than that with land subsidence alone.This implies that the inundation risk for the A1B-S scenario is higher than that of the existing condition, particularly for a storm with a larger return period.Applying the inundation depth and duration derived from the inundation potential analysis and combining the land uses and the corresponding production value of the inundation area, the inundation loss index was assigned.Subsequently, the inundation risk was determined.The average inundation risk (Ravg) was calculated by the area-ratio-weighted average method, and can be expressed as: where Ri represents the inundation risk in the i cell and Ravg represents the average inundation risk.Applying the inundation depth and duration derived from the inundation potential analysis and combining the land uses and the corresponding production value of the inundation area, the inundation loss index was assigned.Subsequently, the inundation risk was determined.The average inundation risk (R avg ) was calculated by the area-ratio-weighted average method, and can be expressed as: A si (8) where R i represents the inundation risk in the i cell and R avg represents the average inundation risk.The risk value of each computational cell in the land subsidence area (Provincial Highway No. 1 to river outlet boundary) was multiplied by the area of the corresponding computational cell.The above-mentioned products for the entire land subsidence area were then accumulated and the sum was divided by the area of the land subsidence area to obtain the average risk value.This value was then used to derive the average inundation risk with one-and two-day storms for each return period.The inundation risks in the presence of land subsidence and land subsidence combined with climate change for each return period are shown in Figures 17 and 18.These figures show that the average inundation risk in the presence of land subsidence alone, or land subsidence combined with climate change with one-and two-day storms, for each return period, is higher than the inundation risk before land subsidence.Moreover, the average inundation risk of land subsidence combined with climate change is higher than that of land subsidence alone.
After land subsidence, the average inundation risk for each return period was 1.08-1.19times greater than that before land subsidence.With the influence of both land subsidence and climate change, the average inundation risk for each return period was 1.14-1.73times greater than that before land subsidence.Therefore, land subsidence combined with climate change seriously increases inundation risk, especially in any scenario with a two-day rainfall event.The risk value of each computational cell in the land subsidence area (Provincial Highway No. 1 to river outlet boundary) was multiplied by the area of the corresponding computational cell.The above-mentioned products for the entire land subsidence area were then accumulated and the sum was divided by the area of the land subsidence area to obtain the average risk value.This value was then used to derive the average inundation risk with one-and two-day storms for each return period.The inundation risks in the presence of land subsidence and land subsidence combined with climate change for each return period are shown in Figures 17 and 18.These figures show that the average inundation risk in the presence of land subsidence alone, or land subsidence combined with climate change with one-and two-day storms, for each return period, is higher than the inundation risk before land subsidence.Moreover, the average inundation risk of land subsidence combined with climate change is higher than that of land subsidence alone.

Conclusions
An inundation risk evaluation model was developed using a PDI model, integrated with inundation depth, duration, and vulnerability.The inundation risk evaluation model analysed the impacts of land subsidence and the combination of land subsidence with climate change on inundation depth, area, volume, and inundation risk in a land subsidence-prone area.The PDI model could effectively simulate the inundation phenomenon after it was successfully validated by typhoon Morakot.The model was further applied to assess the effects of land subsidence and climate change on inundation disaster and risk.The results indicated that the inundation depth, area, and volume in the investigated area increased significantly after land subsidence; the inundation phenomenon was even worse when both land subsidence and climate change were present.Land subsidence increased not only inundation depth but also inundation duration.The average inundation duration for each return period after land subsidence was 1.05-1.1 times greater than that before land subsidence.With After land subsidence, the average inundation risk for each return period was 1.08-1.19times greater than that before land With the influence of both land subsidence and climate change, the average inundation risk for each return period was 1.14-1.73times greater than that before land subsidence.Therefore, land subsidence combined with climate change seriously increases inundation risk, especially in any scenario with a two-day rainfall event.

Conclusions
An inundation risk evaluation model was developed using a PDI model, integrated with inundation depth, duration, and vulnerability.The inundation risk evaluation model analysed the impacts of land subsidence and the combination of land subsidence with climate change on inundation depth, area, volume, and inundation risk in a land subsidence-prone area.The PDI model could effectively simulate the inundation phenomenon after it was successfully validated by typhoon Morakot.The model was further applied to assess the effects of land subsidence and climate change on inundation disaster and risk.The results indicated that the inundation depth, area, and volume in the investigated area increased significantly after land subsidence; the inundation phenomenon was even worse when both land subsidence and climate change were present.Land subsidence increased not only inundation depth but also inundation duration.The average inundation duration for each return period after land subsidence was 1.05-1.1 times greater than that before land subsidence.With the influence of both land subsidence and climate change, the average inundation duration for each return period was 1.13-1.27times greater than that before land subsidence.By further applying inundation depth and duration indices and combining the land uses as well as the corresponding monetary value of the inundation area, the inundation vulnerability index was assigned, which can be used to calculate the inundation risk for different return periods.The average inundation risk for each return period after land subsidence was 1.08-1.19times greater than that before land subsidence.With the influence of both land subsidence and climate change, the average inundation risk for each return period was 1.15-1.73times greater than that before land subsidence.Therefore, it is concluded that land subsidence and climate change significantly change the inundation potential and increase inundation depth, volume, and risk, making inundation disasters more severe.

Figure 3 .
Figure 3. Linbian River land subsidence inspection point distribution and accumulated land subsidence numerical diagram.

Figure 3 .
Figure 3. Linbian River land subsidence inspection point distribution and accumulated land subsidence numerical diagram.
also the damage induced by disasters, even multiple floods have the same inundation depth.Hence, classification of the investigated area is based on the different land uses and vulnerability score indices, as determined by the monetary values of these land uses.

Figures 7 -
indicate that under the influence of land subsidence and land subsidence combined with climate change, the percentage of inundation area and volume with oneand two-day storms, for different inundation depths in each return period, may increase or decrease.Nevertheless, the overall results showed greater percentages of increase for inundation area and volume.Furthermore, the places with greater inundation depths originally tended to show increased inundation areas and volumes after land subsidence.

Figure 8 .Figure 9 .
Figure 8.(a) Spatial distribution of inundation depth increase after land subsidence for a one-day storm under a 50-year return period; (b) Spatial distribution of inundation depth increase after land subsidence combined with climate change for a one-day storm under a 5-year return period.

Figure 8 .Figure 8 .Figure 9 .
Figure 8.(a) Spatial distribution of inundation depth increase after land subsidence for a one-day storm under a 50-year return period; (b) Spatial distribution of inundation depth increase after land subsidence combined with climate change for a one-day storm under a 5-year return period.

Figure 9 .
Figure 9. (a) Spatial distribution of inundation depth increase after land subsidence for a two-day storm under a 5-year return period; (b) Spatial distribution of inundation depth increase after land subsidence combined with climate change for a two-day storm under a 5-year return period.

Figure 11 .
Figure 11.Average flooding duration for a one-day storm under different return periods and scenarios.

Figure 12 .
Figure 12.Average flooding duration for a two-day storm under different return periods and scenarios.

Figure 11 . 21 Figure 11 .
Figure 11.Average flooding duration for a one-day storm under different return periods and scenarios.

Figure 12 .
Figure 12.Average flooding duration for a two-day storm under different return periods and scenarios.

Figure 12 .
Figure 12.Average flooding duration for a two-day storm under different return periods and scenarios.

Figure 13 .
Figure 13.(a) Distribution map of inundation risk increase after land subsidence for a one-day storm under a 5-year return period; (b) Distribution map of inundation risk increase after land subsidence for a two-day storm under a 5-year return period.

Figure 13 .
Figure 13.(a) Distribution map of inundation risk increase after land subsidence for a one-day storm under a 5-year return period; (b) Distribution map of inundation risk increase after land subsidence for a two-day storm under a 5-year return period.

Figure 13 .
Figure 13.(a) Distribution map of inundation risk increase after land subsidence for a one-day storm under a 5-year return period; (b) Distribution map of inundation risk increase after land subsidence for a two-day storm under a 5-year return period.

Figure 14 .
Figure 14.(a) Distribution map of inundation risk increase after land subsidence for a one-day storm under a 50-year return period; (b) Distribution map of inundation risk increase after l and subsidence for a two-day storm under a 50-year return period.

Figure 15 .
Figure 15.(a) Distribution map of inundation risk increase after land subsidence and A1B-S scenario for a one-day storm under a 5-year return period; (b) Distribution map of inundation risk increase after land subsidence and the A1B-S scenario for a two-day storm under a 5-year return period.

Figure 15 .
Figure 15.(a) Distribution map of inundation risk increase after land subsidence and A1B-S scenario for a one-day storm under a 5-year return period; (b) Distribution map of inundation risk increase after land subsidence and the A1B-S scenario for a two-day storm under a 5-year return period.

Figure 16 .
Figure 16.(a) Distribution map of inundation risk increase after land subsidence and A1B-S scenario for a one-day storm under a 50-year return period; (b) Distribution map of inundation risk increase after land subsidence and the A1B-S scenario for a two-day storm under a 50-year return period.

Figure 16 .
Figure 16.(a) Distribution map of inundation risk increase after land subsidence and A1B-S scenario for a one-day storm under a 50-year return period; (b) Distribution map of inundation risk increase after land subsidence and the A1B-S scenario for a two-day storm under a 50-year return period.

Water 2018 ,Figure 17 .Figure 17 .
Figure 17.Average inundation risk for a one-day storm under different return periods and scenarios.

Figure 17 .
Figure 17.Average inundation risk for a one-day storm under different return periods and scenarios.

Figure 18 .
Figure 18.Average inundation risk for a two-day storm under different return periods and scenarios.

Figure 18 .
Figure 18.Average inundation risk for a two-day storm under different return periods and scenarios.

Table 1 .
Inundation depth and inundation duration index.

Table 2
[1]ws the monetary value of each land use according to the investigation conducted by the Sixth River Bureau of the Water Resources Agency[1].The monetary value of output for various land use per unit area and the consumer price index [14] of Taiwan were used to compute the monetary values presented in Table2.According to the table, construction land or buildings have the highest monetary values and fishing farms have the highest monetary values among other agricultural activities.Land uses influence not only monetary values, but

Table 1 .
Inundation depth and inundation duration index.

Table 2 .
[15]e of output for various land use per unit area[15].Value of output for land use per unit area multiplied by Consumer Price Index.NTD is New Taiwan Dollar. ()

Table 3 .
Vulnerability index for different land uses.

Table 4 .
Precipitation data of baseline and scenario A1B-S for one-and two-day storms under various return periods.

Table 4 .
Precipitation data of baseline and scenario A1B-S for one-and two-day storms under various return periods.