Urban Coastal Flood-Prone Mapping under the Combined Impact of Tidal Wave and Heavy Rainfall: A Proposal to the Existing National Standard

: The drivers for coastal ﬂooding may vary from extremely high intensity and persistent rainfall, morphological factors of the coastal area, to extreme waves from the ocean. This means that the ﬂood vulnerability of a coastal area does not solely depend on a single driver but can be a combination with others. A national standard for coastal ﬂooding based on rainfall drivers has been developed. As an evaluation, this study aimed to develop a method for coastal ﬂood-prone mapping by combining rainfall with tidal waves. The steps included the assessment of the coastal ﬂood-prone areas driven by rainfall (CFR) and the coastal ﬂood-prone areas by combined drivers (CFC), which was developed by employing the analytic hierarchy process (AHP), spatial-overlaid, weighted-scored, and logical tests. The coastal area of Mataram City on the Island of Lombok in Indonesia was selected as the study area, since it is frequently a ﬀ ected by ﬂooding. The ﬁndings determined the essentiality of the CFC method for identifying ﬂood vulnerability areas. Thus, the minimum standard for CFC parameters can be deﬁned with climatic and land characteristic factors. Further, the ﬁndings also identiﬁed the need for expert judgment in the development of the CFC weighted score-based method.


Introduction
One of the impacts of climate change is that extreme currents and tidal wave activity has increased flooding rates in coastal and riverine flood plains [1,2]. The intensity of flooding events will extend the degradation of coastal habitats and the devastation of coastal settlements. The occurrence of flooding may be caused by the combination of natural and anthropogenic factors. Natural factors include extremely high intensity and persistent rainfall, morphological factors of the coastal area, and extreme waves from the ocean. Meanwhile, anthropogenic factors could be related to human behaviors in utilizing the coastal environment, such as the conversion of the coastal ecosystem to the unsustainable development of residential areas, tourism, industries, ports, and other exploitations [3].
Coastal flooding usually occurs in low-lying areas, which are flooded by seawater [4]. According to [4], the extent of coastal flooding is a result of the penetration of inland floodwaters controlled by the topography of the coastal land exposed to the flooding. Within coastal areas, flooding can occur not only from overflowing streams caused by heavy rainfall but, also, by inflowing water from the ocean. There are various kinds of flooding, such as tidal flooding, which typically occurs during extremely high tides [5]. Other types of floods can be caused by heavy rains, storm surges, and high waves that arise during coastal storms [6]. The tidal flood in the coastal area of the city of Semarang, Indonesia that led to a coastal recession is one example of the extreme impact of coastal flooding [7,8]. Another example was the West Cliff Drive area in Santa Cruz, CA, USA, which has been moved landward at about 0. 5 m/ year due to the tidal wave activities [9].
Coastal flooding is a global issue that has caused the degradation of the coastal environment and affected the lives of 20 to 200 people who were exposed to flooding [10]. Therefore, the study of coastal flooding and mapping methodology needs to be evaluated. Flooding has occurred frequently in Indonesia; over the past century, 108 times were recorded, accounting for 33.3% of all disaster events [11]. As a consequence, the mapping of flood-prone areas needs to be regularly updated, so that proper planning can be done regarding the inshore urban areas in Indonesia. Flood-prone area mapping within a coastal area should adequately refer to existing standards to support the development of precise, accurate, and standardized spatial information. A national standard, SNI 8197:2015, concerning the mapping method for flood-prone measures at 1:50,000 and 1:25,000 scales with the primary issues of rainfall and overflowing streams have been published and nationally implemented [12]. However, the implementation of the standard is not mandatory. Consequently, the method of flood-prone mapping may vary among scientists, potentially generating different results that may lead to inaccuracy. For example, the standard and [11] used three parameters of flooding-prone criteria: rainfall, slope, and land use. Despite using the same parameters, the methods applied were dissimilar. Ref. [11] used a pairwise comparison approach to specify the weight of each parameter, while the national standard document set the weighted value into 30% and 70% for the rainfall and slope parameters, respectively. Ref. [11] only considered rainfall as the driver of flooding and ignored the existence of ocean drivers such as tidal waves that have been studied and implemented internationally by [1,2,13]. This is not in-line with The General Guidelines for Disaster Risk Assessment that states the importance of wave height, current vegetation cover, coastline shape, and coastal typology for determining a flood-prone assessment [14]. Ref. [15] also states the need for an urban-flooding assessment by considering multiple process drivers. Therefore, the objective of this study is to develop a method for coastal flood-prone mapping by combining the rainfall and tidal waves as the drivers and the land characteristics of the coastal area. Before the development of the combined method, an evaluation of the coastal flooding method caused by heavy rainfall was assessed to develop a better equation as an update to the existing flood-prone mapping method.

Study Area
The coastal area of Mataram City, on the Island of Lombok, Indonesia, was selected as the study area ( Figure 1). Located at 8 • 33 -8 • 38 S and 116 • 04 -116 • 10 E, Mataram City directly faces the Lombok Strait, which separates Lombok Island with the Island of Bali. The Lombok Strait is designated as a channel for the implementation of the right of Archipelagic Sea lane passage based on international maritime law conventions, which connects the Pacific Ocean in the north with the Indian Ocean in the south. The Lombok Strait has a vital role controlling the dynamics of the movement of water masses from the Pacific Ocean to the Indian Ocean [16], which can have an impact on tidal waves on the coast of the city of Mataram.
Lombok Island is also famous as a tourist destination as Bali, but the island is vulnerable to floods, earthquakes, tsunamis, droughts, extreme waves, abrasion, and extreme weather [17]. Based on the land system map of Mataram from Badan Informasi Geospasial, the coastal geomorphological type of the study area is an alluvial plain, with flat topography and sandy clay soil. The two rivers that have the impact of flooding emptying into the coast of Mataram are the Jangkik River and the Unus River. Mataram City is inundated every year, and the impact is exacerbated during monsoons [18]. Tidal waves are the major factor that caused inundation in the coastal city of Mataram [19][20][21]. A major flood occurred in January 2017 where tidal waves inundated 1161 houses on the coast of Mataram [18]. Research by [18] has projected that, within the next 50 years, the coast of Mataram City will likely be inundated without any mitigating actions.

Materials and Methods
This study assessed two methods of flood-prone mapping: coastal flood-prone areas driven by rainfall (CFR) and coastal flood-prone areas by combined drivers (CFC). The latter method combined a series of coastal flood-prone methods in the existing national standard while adding tidal waves as a driver from the ocean. The development of both methods employed multi-analysis criteria, including pairwise comparison using the analytic hierarchy process (AHP) method, overlay method, weighted-scored method, and logical tests to validate the results.
The data for the study consisted of: (a) Tidal data obtained from the website of Tide Observation and Prediction-Geospatial Information Agency (we referred to BIG: "Badan Informasi Geospasial" in Indonesian) (http://tides.big.go.id). The average of highest tidal data was assessed after data acquisition. The coastal alluvial plain was used as the object of the mapping to examine the extent of flooding, as morphological descriptors are good indicators of flood hazard exposure [24][25][26]. In this case, the alluvial plain of the Mataram Coast is a critical geomorphological unit due to its vulnerability to flooding.
The workflow of this research is summarized in Figure 2.

Materials and Methods
This study assessed two methods of flood-prone mapping: coastal flood-prone areas driven by rainfall (CFR) and coastal flood-prone areas by combined drivers (CFC). The latter method combined a series of coastal flood-prone methods in the existing national standard while adding tidal waves as a driver from the ocean. The development of both methods employed multi-analysis criteria, including pairwise comparison using the analytic hierarchy process (AHP) method, overlay method, weighted-scored method, and logical tests to validate the results.
The data for the study consisted of: (a) Tidal data obtained from the website of Tide Observation and Prediction-Geospatial Information Agency (we referred to BIG: "Badan Informasi Geospasial" in Indonesian) (http://tides.big.go.id). The average of highest tidal data was assessed after data acquisition. The coastal alluvial plain was used as the object of the mapping to examine the extent of flooding, as morphological descriptors are good indicators of flood hazard exposure [24][25][26]. In this case, the alluvial plain of the Mataram Coast is a critical geomorphological unit due to its vulnerability to flooding.
The workflow of this research is summarized in Figure 2. The stages are explained as follows:

Evaluation of the Existing Method
Prior to the assessment of the CFC, this study needed to examine the existing national standard (CFR-2) as the basis of the combined driver model (CFC). Previously, the national standard [12] used rainfall to represent the climate and slope to represent the land characteristics but assumed the land cover parameter along the coastal city was homogeneous. However, [27][28][29][30] highlighted the importance of urban sprawls as an increasing driver of urban flood risk and indicated the relevance of land use for assessing the potential of flood hazard impacts in coastal areas. Ref. [31] supported the role of land cover/use for a flooding assessment in Mataram City, because 8.59% of the coastal built-up area was a conversion from agricultural land (rice fields). Besides, [32] stated that inundation adaptation should become a land-use issue rather than an engineering perspective. Therefore, the importance of the land-use cover parameter cannot be discounted. CFR-2 was modified to CFR-3 by adding the land cover/use parameter (see the matrix of multi-criteria evaluation (MCE) in Table 1) Table 1. Shows the matrix of MCE for coastal flood-prone areas driven by rainfall (CFR)-2 and CFR-3 based on [12]. Notes: w = weight and S = score. The stages are explained as follows:

Evaluation of the Existing Method
Prior to the assessment of the CFC, this study needed to examine the existing national standard (CFR-2) as the basis of the combined driver model (CFC). Previously, the national standard [12] used rainfall to represent the climate and slope to represent the land characteristics but assumed the land cover parameter along the coastal city was homogeneous. However, [27][28][29][30] highlighted the importance of urban sprawls as an increasing driver of urban flood risk and indicated the relevance of land use for assessing the potential of flood hazard impacts in coastal areas. Ref. [31] supported the role of land cover/use for a flooding assessment in Mataram City, because 8.59% of the coastal built-up area was a conversion from agricultural land (rice fields). Besides, [32] stated that inundation adaptation should become a land-use issue rather than an engineering perspective. Therefore, the importance of the land-use cover parameter cannot be discounted. CFR-2 was modified to CFR-3 by adding the land cover/use parameter (see the matrix of multi-criteria evaluation (MCE) in Table 1) Table 1. Shows the matrix of MCE for coastal flood-prone areas driven by rainfall (CFR)-2 and CFR-3 based on [12]. Notes: w = weight and S = score.
The equation of the matrix is described as follows: where F i is the flood-prone area of class i, w i is the weighting for the flood-prone area of class i, and s i is the scoring for the flood-prone area of class i.

CFR-New method
A spatial analysis that employs the weighted-score method has a critical issue in the weighting variable, because the values are utilized to calculate the level of priority of an analysis involving multiple attributes [33][34][35]. Thus, to assess the object classification based on the weighted score method, it is essential to learn the weighting factors, especially when multiple features are used to distinguish variable object classes. The AHP method was used to specify the weight of the three parameters in the flood-prone area criteria (CFR-3) to develop the new CFR method. Firstly, the Saaty method, as quoted in [36], was used to identify the hierarchy structure of used parameters with a scale of 1 to 9. Secondly, a pairwise comparison matrix was developed and normalized to synthesize the priority level by calculating the eigenvector to obtain the weight values of the used criteria. The equations of the methods can be expressed as follows .. a 1n a 21 where A 1 to n is the average of (a 1n +..+a mn ) as the eigenvectors; a 11 to a mn is the value of the performance rating towards the parameters.
Then, the weighted decision was calculated based on the normalized matrix and eigenvectors as follows: where w is the weighted value of parameters 1 to n.
Thirdly, the consistency of the weighting value assessment-Saaty's consistency ratio (CR), as quoted in [37]-was assessed. The CR is a comparison of the consistency index (CI) with a random index (RI) value, as follows: where CI is determined by λ max (the largest eigenvalue) and the number of parameters, n, If a pairwise comparison matrix has a CR ≤0.100, then the inconsistency in the decision-makers' opinions is still acceptable; otherwise, the assessment needs to be repeated.
If CR ≤ 0.100, then the spatial analysis will be employed by using the following equations: where F i is the flood-prone area of class i, w 1..n are the weighting factors for the parameters p 1 to p n , and S i is the scoring value for the flood-prone area class i for the parameters p 1 to p n .

CFC Method Development
The most crucial issues of this study were to develop the method of CFC mapping and to assess how much the potential impact is of tidal waves on coastal flooding. The AHP method was also used to specify the weights of the parameters, assessed according to Equations (2)- (6).
An important thing that needed to be considered in the development of the CFC was the extent to which a tidal wave would inundate the coastal area. The spatial analysis for the tidal wave over wash or landward inundation (x r ) can be defined as a function of the height of the tidal wave (h c ), the depth at the still-water line (h o ), and the difference between the still-water level and the normal tide x p , while θ is the angle between the water level when tidal wave landward inundation occurs with the beach slope (see Figure 3). Next, the CFC method can be determined following Equation (7).

CFC Method Development
The most crucial issues of this study were to develop the method of CFC mapping and to assess how much the potential impact is of tidal waves on coastal flooding. The AHP method was also used to specify the weights of the parameters, assessed according to Equations (2)- (6).
An important thing that needed to be considered in the development of the CFC was the extent to which a tidal wave would inundate the coastal area. The spatial analysis for the tidal wave over wash or landward inundation ( ) can be defined as a function of the height of the tidal wave ℎ ), the depth at the still-water line ℎ ), and the difference between the still-water level and the normal tide , while θ is the angle between the water level when tidal wave landward inundation occurs with the beach slope (see Figure 3). Next, the CFC method can be determined following Equation (7).

Implementation of the Methods in Mapping
Equation (6) was used to determine the interval of the vulnerability classes, which can be explained in Table 2.

Vulnerability
Classes

Implementation of the Methods in Mapping
Equation (6) was used to determine the interval of the vulnerability classes, which can be explained in Table 2.

Flood-Prone Area
Class Interval Vulnerability Classes Flood-Prone Classes F i =w 1 s ip1 +· · · +w m s ipn Notes: C _low = low vulnerability, C _mod = moderate vulnerability, and C _high = high vulnerability.
Afterward, the map of the coastal flood-prone area for CRF-2, CFR-3, CFR-new, and CFC were developed by using the overlay analysis using GIS software.

Field data inventory
The field data inventory was carried out as a reference to evaluate the results of the study due to the unavailability of a recent detailed flood disaster map for the study area. The field data consisted of interview data from local people in the identified flood-prone areas and the point data of the flood locations. The respondents for the in-depth interviews were randomly selected from four clusters of flood-prone areas, which were identified based on the coastal village administration, as suggested by The Regional Board for The Disaster Management Office. Three respondents were selected from each cluster, followed by a field-mapping of the flood-prone area and the tracking of the flood-prone boundaries. The contents of the interview dealt with the discussion of flooding events and at least consisted of when, where, how long, how significant the flooding, the frequency of flooding/how often per year, when was the worst flooding, how much it has affected the people's lives, and how is the emergency response and mitigation works.

Evaluation
The evaluation was carried out by assessing the findings with the data obtained from the field survey using the logical test method. Both the references and the study derived map were analyzed according to the arranged matrix; then, the test was carried out following the method described by [38] as follows: i f (x, y) 1 = (x, y) n , then (x, y) 1 is true (8) i f (x, y) 1 (x, y) n , then (x, y) 1 is f alse (9) where (x, y) n is a reference point data of the flood location. The total evaluation value (Ev) can then be determined as where Ev is the result of the evaluation (%) based on t = a true value and f = a false value, 1 . . . n are the flood-prone points from the study, and P ref is the field reference data. The logical test analysis was applied for both the classified CFC map into three vulnerable classes and the not-classified derived map.

Coastal Flood-Prone Area Method Based on Rainfall
The results that compared the standard with two parameters (CFR-2) and the CFR-3 show that, by adding the land cover parameter to the method, different levels and extents of the potential disaster can be identified (see Figure 4). The highly vulnerable area increased with the CFR-3 method (b) compared with the CFR-2 method (a) (Figure 4).
where is the flood-prone area of class i, is the scoring for parameter = land use at class i, is the scoring for parameter = slope at class i, and is the scoring for parameter = rainfall at class i. However, because the consistency of the original weighting remained unclear, we updated the weighting variable and reassessed the CFR-3 to become the CFR-new. The CFR-new method was developed based on an AHP analysis, which used a pairwise comparison matrix to determine the level of importance (priority scale) of the used parameters. A judgment matrix was constructed for CFR-new, as shown in Table 3, which summarizes the normalized weights. The matrix (λ max = 3.0387; CI = 0.0194; RI = 0.58; CR = 0.0334 < 0.1) passed the consistency evaluation. So, the measured weighting values can be used for the development of the new CFR equations (CFR-new). A multi-criteria evaluation (MCE) matrix, as shown in Table 4, and the prototype of the CRF-new map (in Figure 5) are the results of the implementation of the CFR-new equation (Equation (11)).   The findings show slight differences between the CFR-3 and the CFR-new map. The coverage of the high-vulnerability class in the CFR-new was 2.56 ha larger than the CFR-3. Conversely, the moderate-vulnerability class was 179.98 ha and 182.54 ha, respectively ( Figure 5 and Table 5). Table 5. Flood-prone area of CFR-2, CFR-3, and CFR-new.

The new equation for CFR-new can be defined as follows:
F i = 1.20s li p l + 0.20s si p s + 0.03s ri p r (11) where F i is the flood-prone area of class i, s li is the scoring for parameter p l = land use at class i, s si is the scoring for parameter p s = slope at class i, and s ri is the scoring for parameter p r = rainfall at class i. The findings show slight differences between the CFR-3 and the CFR-new map. The coverage of the high-vulnerability class in the CFR-new was 2.56 ha larger than the CFR-3. Conversely, the moderate-vulnerability class was 179.98 ha and 182.54 ha, respectively ( Figure 5 and Table 5). It can be said that the application of the CFR-new method was more accountable compared to the CFR-2 and CFR-3 methods, because it was applied to the measured and proven weight values of each parameter's criterion. Thus, it was essential to measure the weight before conducting the spatial analysis processes.

Coastal Flood-Prone Area Method Based on Combined Drivers
Following the use of AHP to develop the CFR-new, a judgement matrix was constructed for the CFC, as shown in Table 6, which summarizes the normalized weights. The matrix (λ max = 4.01185; CI = 0.0395; RI = 0.90; CR = 0.0439 < 0.1) meets the CR acceptable criteria described by Saaty [39,40]. Table 6. AHP judgement matrix of the coastal flood-prone areas by combined drivers (CFC)-affecting parameters and consistency ratio assessment. Then, the CFC equation was developed as follows:

Slope
F i = 1.24s si p s + 0.28s wi p w + 0.06s li p l + 0.01s ri p r (12) where F i is the flood-prone area of class i, s li is the scoring for parameter p l = land use at class i, s si is the scoring for parameter p s = slope at class i, s wi is the scoring for parameter p w = tidal wave at class i, and s ri is the scoring for parameter p r = rainfall at class i. The average tidal wave height (h c ) on Lombok island is 1.00 to 1.25 m, with the extreme tidal waves reaching >2 m. By considering the extreme and average heights, the terrestrial inundation (x r ) caused by the tidal wave ranged from 14.93 m to >47.68 m. An MCE matrix, as shown in Table 7, and the prototype of the CFC mapping (in Figure 6) was developed based on the CFC equation (Equation (12)).  The findings show considerable differences when compared to the CFR, because the implementation of the CFC method classified almost all of the area as vulnerable (Figure 6). Table 8 shows the differences in the extents of the three models based on the MCE analysis, which significantly expanded the inshore boundary in the CFC model implementation by as much as 244.02 ha for the high-vulnerable class. Table 8. Flood-prone area analyses from the (a) CFR-3 parameters, (b) CFR-new (c), and CFC combined drivers.

Prone Classes
Area (  The findings show considerable differences when compared to the CFR, because the implementation of the CFC method classified almost all of the area as vulnerable (Figure 6). Table 8 shows the differences in the extents of the three models based on the MCE analysis, which significantly expanded the inshore boundary in the CFC model implementation by as much as 244.02 ha for the high-vulnerable class.

Discussion
The methodology established in this study highlights the minimum parameters needed for coastal flood-prone area mapping and the role of AHP for developing an equation that is necessary to specify the vulnerability classes. It shows that climate and land characteristic parameters are important for determining flood-prone area mapping. Several researchers [41][42][43] agree that land characteristics (such as land cover, topography/slope, and geomorphology) and climate (such as rainfall and extreme tidal waves) are influential in determining the occurrence of floods. The findings show how the integration of these parameters can determine the worst impact of the extent of inundation across coastal alluvial plains ( Figure 6 and Table 7). The explanation of the added parameters and the methods are as follows: a.
Land cover is essential for the CFR and CFC methods: [30] shows the empirical link of flooding to land use and land cover conversion that plays a vital role in the occurrence of flooding. Table 5 showed the increase of high and moderate vulnerability classes when the land cover/use parameter was added, with 42.63 ha for the high class and 125.56 ha for the moderate class from CR-2 to CR-3, respectively. b.
The CFC method allows more detailed impacts on the land: [44,45] supports that coastal and deltaic areas are susceptible to flooding from the combined effects of excess rainfall and astronomic tide inundations. Table 8 describes the differences in susceptibility of the moderate and low classes in CFR-3 and CFR-new, with 142.74 ha and 140.18 ha, respectively, to CFC while an extreme tidal wave was added to the model. c.
Flooding depends on geomorphological units: [46] supports the importance of hydro-geomorphological aspects for coastal flood mapping in that an alluvial plain can be used to delineate flood hazard areas based on his research in Constantine City, Algeria. The geomorphological type of the Mataram City coast is an alluvial plain. As can be seen in Figure 6, using either the CFR or CFC method, the alluvial plain area in the study site is vulnerable to flooding and can be used as a mapping unit. The observation area indicates that the highest tidal wave caused a dynamic form of the alluvial plain, which sometimes created a sand embankment along the coast of Mataram City. This condition may intensify the inundation of the area, especially in low areas that are blocked by sand embankments. In his research in Vietnam, [47] indicated that flooding and the dynamic mechanism of alluvial formation were correlated and, in turn, affected the flooding hazard.
Another important finding is the role of expert judgment through AHP to develop an equation for coastal flood-prone mapping. Generated thematic maps by using an overlay analysis can be based on the weighted scores of each factor in the AHP model [48]. Ref. [49] supports the importance of weighting values as a central step in the multi-criteria analysis process. AHP assists in decision-making by organizing quantitative and qualitative data in an integrated, transparent, and structured way [50]. Based on the weights and scores derived using AHP, coastal flooding vulnerability maps can be assessed as low, medium, and high-vulnerability classes (Figures 5 and 6). Therefore, it is essential to measure the weight before conducting spatial analysis processes.
Other issues that should be evaluated is the selection of the drivers of the flooding despite the land characteristic issues and rainfall. Our research selected the tidal wave as the ocean driver, because tidal flood events have been increasing recently in Indonesia [10]. The northern coast of Java Island has been the worst affected by this phenomenon [10]. The implementation of a combined method has also been done by many researchers worldwide. For example, several researchers assessed storm tides and heavy rainfalls as a combined driver in their studies, sequentially, [51] for Crookhaven Head and Shoalhaven Head Estuaries, NSW, Australia, [52] North Folk, VA, USA, and [53] for Saint Louis in the Caribbean Sea. However, the method they used was the hydrological model without any level of flood vulnerability classes. An AHP application for flooding and its vulnerability class research was employed by [35] on the Coast of Lithuania in the South-Eastern Baltic Sea. However, his research did not apply the combined drivers but only the micro-tidal wave drivers. Therefore, the combined drivers supported by the AHP method for the development of the equation-based weighted scored method has become important to enrich the coastal flood-prone area mapping technique.

Evaluation of the CFR-new and CFC methods
The development of the equations for both methods, CFR-new and CFC, were found consistent by using CR analysis. The extent of the accuracy of the prototype with the ground conditions was evaluated by using the logical test method to compare the findings with the field data as a reference. The reference data represents actual flood events occurring in the alluvial plain area. Without distinguishing flood classes, about 97% of the alluvial area is vulnerable to flooding, utilizing either CFR-new = 96.25% or CFC = 97.5%. However, in considering the classification of vulnerability, the results of the evaluation were dissimilar; there was 51.25% conformity obtained for CFR-new and 98.75% for CFC (Figure 7).  Ref. [20] concluded that river discharge was largely responsible for flooding along the river and the river mouth. The [18] study supports that a tidal wave would cause higher flooding along the river mouth of the Ampenan and the Sekarbela subdistricts. Thus, it can be said that the CFR-new value was lower because the high-vulnerability classes of the CFC had a strong relationship with the river mouth. Both CFR and CFC found high vulnerability classes around the Jangkik River in the Ampenan subdistrict and the Unus and Berenyok Rivers in the Sekarbela subdistrict.

Conclusions
Some studies have assessed coastal flooding using tide and rainfall data independently. However, within coastal areas, flooding cannot be caused solely by a single factor. Ref. [52] reported that heavy rainfall, alongside the storm tide, exacerbated the severity of coastal flooding. Indeed, [54] supported that rainfall and tidal waves contribute to coastal inundation. The combination of rainfall, tidal wave, geomorphology, and other land characteristic factors are correlated as the causes of coastal flooding. This study proposed a method and prototype of CFC using the coastal area of Mataram City, Lombok, Indonesia as an example. The findings showed that the combined model increased the identification of flood-prone areas compared to the single-driver model. The minimum parameters for determining the flood-prone areas can be identified from this study; there is the slope, Ref. [20] concluded that river discharge was largely responsible for flooding along the river and the river mouth. The [18] study supports that a tidal wave would cause higher flooding along the river mouth of the Ampenan and the Sekarbela subdistricts. Thus, it can be said that the CFR-new value was lower because the high-vulnerability classes of the CFC had a strong relationship with the river mouth. Both CFR and CFC found high vulnerability classes around the Jangkik River in the Ampenan subdistrict and the Unus and Berenyok Rivers in the Sekarbela subdistrict.

Conclusions
Some studies have assessed coastal flooding using tide and rainfall data independently. However, within coastal areas, flooding cannot be caused solely by a single factor. Ref. [52] reported that heavy rainfall, alongside the storm tide, exacerbated the severity of coastal flooding. Indeed, [54] supported that rainfall and tidal waves contribute to coastal inundation. The combination of rainfall, tidal wave, geomorphology, and other land characteristic factors are correlated as the causes of coastal flooding. This study proposed a method and prototype of CFC using the coastal area of Mataram City, Lombok, Indonesia as an example. The findings showed that the combined model increased the identification of flood-prone areas compared to the single-driver model. The minimum parameters for determining the flood-prone areas can be identified from this study; there is the slope, land cover/use, rainfall, and tidal waves. The findings also indicate the need for expert judgment for the development of equations that involve multi-criteria analysis and so does the assumption of using an alluvial plain as the unit for mapping. Therefore, this study may become a useful input for the existing national standard, which only considers rainfall as the affecting factor for flooding. Thus, the equation derived from the AHP analysis for the mapping technique based on the weighted-score method should also be considered to update the national standard.
Besides that, there are still other factors affecting flooding phenomenon that may enrich the scientific realm in the field of coastal flooding. Tidal waves are not the only factor that can cause coastal flooding. Sea levels rising, storm surges, and other high waves are additional drivers that should be studied due to their impacts on coastal inundation. The characteristics of the drivers are also regionally dependable. Therefore, it is likely that the study of multiple drivers for coastal flooding can be further developed.
Author Contributions: Dewayany Sutrisno and Ati Rahadiati were the main contributors of this article, and both contributed to the conceptualization, methodology, analysis, validation, writing the draft, and article preparation, while Aninda W. Rudiastuti, Ratna Sari Dewi, and Munawaroh were the other contributing members. Aninda W. Rudiastuti contributed to the tidal wave data analysis and field data analysis, Ratna Sari Dewi contributed to the field data analysis, and Munawaroh contributed to the rainfall data analysis. All authors have read and agreed to the published version of the manuscript.

Funding:
This research was funded by the Geospatial Information Agency, State Budget No 083.01.01.3539.967.001.052.C: The Research on Coastal Vulnerability Mapping.