Adjustability of exposed elements by updating their capacity for resistance after a damaging event: application to an earthquake–tsunami cascade scenario

The 2011 Great East Japan earthquake and tsunami remind us once again that these types of cascade event can occur and cause considerable damage. The scientific community realizes the need for rapid theoretical and practical progress on cascade events to provide field teams with the necessary tools and information for action during these types of events. The earthquake damage scenario for Martinique and Guadeloupe islands (French West Indies) has already been performed within the framework of French governmental projects, but these areas, in the vicinity of the French West Indies subduction zone, are also subject to tsunami events. In this study, we propose to perform a combined scenario in which an earthquake is followed by a tsunami, as it could arrive one day, considering the seismic characteristics and potential of such a subduction zone. The vulnerability of the buildings is defined considering local specific information based on several years of field inventories and inspections and is later classified into one of the 36 model building types of HAZUS. The calculation of the damages due to tsunamis follows the HAZUS methodologies. The main novelty of our study is the calculation of damage due to the two phenomena occurring one after the other, not in parallel, as is calculated in the existing literature. Therefore, for the calculation of the damages due to the second event (i.e. the tsunami), the vulnerability characteristics of the initial structure are reduced, considering the damage state of the construction after the first event (i.e. the earthquake). Hence, in our case, this calculation approach allows us to update the number of exposed elements and their changed vulnerabilities considering the damages due to the earthquake, since certain structures are already damaged by the earthquake before the arrival of the tsunami wave. The results coming from our study and our manner of treating the cascading hazards are putting into perspective with the Hazus method for combining damages coming from earthquake and the damages coming from consequently tsunami. The results expressed as the sum of the damages in both most damaged states, Extensive and Complete, are more or less in the same range of values for both studies (our study and HAZUS 2017). However, a trend of having more percentage of complete damages (and hence, less the Extensive damages) with our method than the ones obtained with the Hazus combination can be important information for crisis managing. This is a first result for the French West Indies territory, but anyway, more studies should be carried out in order to check this trend and eventually to confirm and validate this issue for others territories with others bathymetries, vulnerabilities and seismological features.


Introduction
Before addressing the question of the damages, due to the several quasi-simultaneous events, the methodologies for calculating the damage of each of these events should be developed and managed well. Next, the treatment of cascade effects from an event to other events and their interactions should be considered to give the most accurate prediction of the damage. The development of the methodologies for calculating the damage due to earthquakes started in the 1970s with the damage probability matrix of Whitman (Whitman et al. 1974) after the 1906 San Francisco Earthquake and currently continued to be improved. Based on this long experience, the methodology for performing the classic damage earthquake scenario is now well known, although each step of the methodology (e.g. shake maps, site effects, vulnerability assessments) can still be improved by, for example, the validation of the numerical results with the post-earthquake observed damage. A number of publications and software have been published and developed with this aim (e.g. GEM Technical Report Crowley et al. 2010a,b;Molina et al. 2010;Hancilar et al. 2010;FEMA 1999;FEMA 2004;Sedan et al. 2013).
Regarding the methodologies for calculating damages due to tsunami events, significant resources have been dedicated worldwide for improving hazard and damage models for tsunamis (Suppasri et al. 2015(Suppasri et al. ,2012 following recent large tsunamis (e.g. Murao and Nakazato 2010;Yamazaki andCheung 2011 andAmakuni andTerazono 2011). More recently, the 2011 Tohoku earthquake tsunami in Japan has provided many data concerning structures (including engineered structures). Hence, the development of a fragility function for tsunamis increased after this event. There were 19 fragility functions derived by the data of this event compared to 11 fragility functions developed for all previous tsunamis combined (Charvet et al. 2017). Recent studies have tried to clarify or develop the different steps necessary to complete the methodology for damage estimations due to tsunami events, including the definitions of (1) the tsunami intensity measurement, (2) the inventory of the exposed elements and of their important characteristics, (3) typologies and their main characteristics and (4) the development or harmonization of the damage description. Park and Cox (2016) and Park et al. (2017) presented the probabilistic tsunami damage assessment and evaluated the influence of five intensity measures of tsunamis applied to Seaside, Oregon.
One of the most recent and complete studies that treat all these aspects was presented by Charvet et al. (2017), but several other previous studies are also very detailed and useful for understanding the main issues related to the tsunami damage scenario (Kircher and Bouabid 2014;Macabuag and Rossetto 2014;Tarbotton et al. 2015;Rehman and Cho 2016). Macabuag et al. (2016) went further in the details of analysis by correlating the damage degree distribution with the inundation depth distribution and building footprint distribution. Jaimes et al. (2016) proposed a new probabilistic tsunami risk assessment approach and exemplified it for public schools in southern Mexico. Browning and Thomas (2016) presented a method for determining an initial assessment of tsunami risk, with application for two coastal areas of Oman.
The work proposed in this paper addresses the cascading effects and interactions in the earthquake-tsunami damage scenarios (Fig. 1, right). With regard to the damage due to the quasi-simultaneous cascade events-in our case, an earthquake followed by a tsunami-not 1 3 many studies are available. When considering the combination of other phenomena, such as earthquakes, landslides or inundations, the SYNER-G (Systemic Seismic Vulnerability and Risk Analysis for Buildings, Lifeline Networks and Infrastructures Safety Gain, www. vce.at/SYNER -G) Project, funded by European Union's Seventh Framework Programme (EU FP7), is a forerunner in the treatment of the combined events (Pitilakis et al. 2014). Of course, this topic of the damages due to the combined events has already been studied in the literature; however, the approaches used are more theoretical or propose a framework for how to consider the cascade or domino effects, and are rarely implemented in a software package or applied to a real case study, as is proposed in this paper. The Bayesian treatment of the cascading effects is modish and the advances proposed by this approach are valuable; however, most often they remain theoretical or applied only to network infrastructures. Gasparini and Garcia-Aristizabal (2014) discuss the role of the cascading effects triggered by the earthquake and its importance for a holistic view of the seismic risk assessment process. Mignan et al. (2016) treated the question of the cascading effects from the social point of view for assessing the adaptive capacities of the communities. The most representative studies that consider the damages both from the earthquake and from the induced tsunami are the paper related to the Great East Japan Earthquake (March 11th, 2011 Tohoku Earthquake). Considering the damage analysis and risk assessment, some studies are related to the 2004 tsunami in Phang Nga and Phuket, Thailand (Ruangrassamee et al. 2006;Rossetto et al. 2007;Römer et al. 2012). Nevertheless, these studies are either observation reports or retrospective questions related to prevention and mitigation measures (Epstein 2011). Generally, these studies treat the impact of the tsunami and the impact of the earthquake separately. Garcia-Aristizabal et al. (2015), within the EU FP7 CRISMA project (Modelling crisis management for improved action and preparedness, www.crism aproj ect.eu), implemented a concept model and tool for evaluating cascading effects into scenario-based analyses.
The scientific challenge comes from the fact that the tsunami-triggered earthquake is at the boundary of two specialties: the tsunami propagation is treated by fluid mechanics, and seismic faults and the building damages are treated by solid mechanics. Tsunami propagation is generally studied in the free field, and when the buildings are included, they are generally considered as non-deformable blocs or barriers, while in the case of the damage of buildings due to earthquakes, the buildings are deformable objects. Note that the terminology used to classify fragility functions, intensity measure or damage scale can be different from one study to another for the same phenomenon (e.g. Tarbotton et al. 2015;GEM Crowley et al. 2010a, b;FEMA-HAZUS 1999;FEMA-HAZUS 2004;Prieto et al. 2018), and even more so when the phenomena are different. Tarbotton et al. (2015) reviewed the existing literature on tsunami fragility curves, noting trends and comparing existing fragility curves to highlight variability in the mean function across a range of studies. Charvet et al. (2017) presented a critical review and discussed the key issues in the current literature for each of the model components: building damage data, tsunami intensity data and the statistical model that links the two. They also presented a very complete table with the published empirical fragility functions for the 1993 Hokkaido-Nansei-Oki Tsunami, 2004Indian Ocean Tsunami, 2009 Samoa Tsunami and 2010 Chilean Tsunami. As described in their table, the "Tsunami Intensity Measures" are numerous (h, inundation depth; v, velocity; F, drag force; MF, momentum flux; MMF, moment of momentum flux adn FQS, a new proposed quasi-steady force estimate), and the number of the damage degrees can vary from 2 to 6. Generally, the fragility functions derived from data for the Great East Japan Earthquake and Tsunami use the damage scale proposed by the Japan Cabinet Office (2013). The equivalence or convergence of terminology from one study to another is a research topic in itself that will not be discussed in this paper. To overcome this debate, we chose to use the HAZUS guidance throughout this paper for characterizing the intensity measure, the building typologies and the damage scale, even if the software is conceived as such that, without any additional development, other terminologies can be applied in our research. The HAZUS (2017) choice is explained by the fact that earthquake and tsunami terminologies are self-consistent and by the fact that this document is the most recent one (November 2017), and can be considered as the most up-to-date development ( Fig. 1, left). Kircher and Bouabid (2014) describe new functions for determining the probability of damage to buildings and essential facilities due to tsunami inundation (flood) and tsunami lateral force (flow) hazards, and for combining the damages due to these hazards with that due to earthquake shaking (i.e. for evaluation of damage and loss due to a local tsunami). The new building damage and loss functions were developed as part of a FEMA-funded project to develop a Tsunami Model for incorporation into the HAZUS Loss Estimation Technology. The Tsunami Model uses the same model building types, occupancy classes and building damage states as those of the Earthquake Model, descriptions of which may be found in the HAZUS technical manuals for earthquakes (FEMA-HAZUS 1999) and the Advanced Engineering Building Module (FEMA-HAZUS 2004). This version of HAZUS (HAZUS 2017) is able to run two types of damage analysis: both near-source (earthquake + tsunami) and distant-source (tsunami only). In this paper, we are interested in the combination of the earthquake and tsunami effects, called "near-source" analysis by HAZUS. Inputs to the calculation of tsunami damage to buildings include hazard data, namely, tsunami flood (inundation height) and tsunami flow (momentum flux), and earthquake damage data from the Earthquake Model (when evaluating local tsunami effects).
For the cascade event (i.e. earthquake followed by tsunami), we consider the simulated damage to structures due to the initial earthquake, and we return to the vulnerability function that we modify in order to consider the increase in its vulnerability regarding the damage degree achieved by the structure after the earthquake. In so doing, we consider this new increase vulnerability of the structure for the calculation of the tsunami damage. This is quite innovative for a scenario approach, and it is important to account for it for such near-source events, since there is likely to be increased vulnerability of some buildings to tsunami loading due to ground shaking, induced weakening of structures and the presence of loose debris that may be entrained in the tsunami flow (Fraser et al. 2014). Thus, the exposed elements are updated after the damage seismic simulations, both in number and in the vulnerability characteristics. As in the HAZUS Tsunami Model from FEMA (2017), we combine earthquake and tsunami inventory attributes, and aggregate and analyse at the census block; however, the update of the exposed element is done by each typology in each census block. The census blocks used in our analysis were established after several years of field inventories (Bertil et al. 2009a;Roulle et al. 2010;Bertil et al. 2009b for Guadeloupe andBelvaux et al. 2013 for Martinique).
The Lesser Antilles arc is associated with the south-westward subduction of the North American plate under the Caribbean plate at a low convergence rate (i.e. 2 cm/ year, DeMets et al. 2000). Moderate to large earthquakes can occur, related either to intraplate active faults or to the subduction process. Damaging historical events have been reported (e.g. SISFRANCE, BRGM 2018; SISFRANCE is the database of historical earthquakes felt in French territories). Moreover, the tsunamis of the French West Indies are mainly associated with regional earthquakes (e.g. Lambert and Terrier 2011), even if some teletsunamis (also called a transoceanic tsunami) were reported (e.g. Zahibo and Pelinovsky 2001;Lambert and Terrier 2011), in particular the tsunami induced by the great 1755 Lisbon earthquake. A teletsunami is a tsunami that originates from a faraway source, which is generally thousands of kilometres from the area of interest (Yu et al. 2011). In this study, we are paying attention to tsunamis induced by regional significant earthquakes. Two large destructive thrust interplate earthquakes occurred in the area: the 11 January 1839, east of Martinique Island, and the 8 February 1843, east of Guadeloupe Island (e.g. Bernard and Lambert 1988;Feuillet et al.2011). Whereas these two damaging earthquakes, among the largest events reported in West Indies, are considered as reference events for seismic hazard analysis, their characteristics (i.e. magnitude, hypocentre or rupture description) still generate debates (e.g. Hough 2013). The 1843 earthquake is the largest known historical event associated with the Lesser Antilles subduction, contemporary testimonies gave reports of intensity  (Bernard and Lambert 1988). In this case been proposed for this event, from 7.5 to 8 (Bernard and Lambert 1988), to 8.5 or higher (Feuillet et al. 2011;Hough 2013). Whereas most of Mw ≥ 8.5 megathrust earthquakes are associated with significant tsunamis, depending on source characteristics and on local conditions, some megathrust events were not associated with major tsunamis (e.g. Nias Mw8.6 2005 earthquake, Briggs et al. (2006). Le Roy et al. (2017) had used such hypothesis in order to model the moderate tsunami generated by the rupture of the deeper interface during 1843 earthquake. As the aim of the present study is not to analyse historical events, but to construct a realistic scenario for a cascade earthquake and tsunami damage estimation, sources characteristics for 1843 and 1839 events proposed by Feuillet et al.(2011) were used. They proposed rupture along the subduction interface for the 1839 and 1843 events associated with Mw8 and 8.5, respectively, and the source parameters they proposed generate significant near field tsunamis (see Table 1).
The main outcome of this study is the consideration of the earthquake damage that may precede tsunami damage and compromise building resistance, in case of nearsource events. This is all the more important, especially for EQ-vulnerable buildings, as is the case in the French Lesser Antilles area. From a practical viewpoint, this implies the update of the exposed elements in the scenario, considering the damages due to  the earthquake. Future studies should also assess tsunami risk to vital lifeline components, such as bridges and roads. Debris generated and casualty estimates also need to be investigated in future studies.
2 Co-seismic-tsunami cascade approach with exposure updating due to earthquake damages Compared with the HAZUS methodology for near-source earthquakes, where the total damage occurs as the result of a combination of earthquake and tsunami hazards calculated separately without interaction ( Fig. 1, left), our approach calculates the damages in one event after the other. This approach allows the updating of the exposed elements and especially their vulnerability, based on the damages caused by the earthquake, for the tsunami damage calculation (blue arrows in Fig. 1, right). This aspect is very important in the areas in which the earthquake vulnerability of structures is very low, as in our case study.
Considering the phenomenon as a cascade effect seems more realistic than treating each hazard separately and then proposing methods for combining the probability of building damage due to a tsunami with the probability of building damage due to the earthquake that generated the tsunami. Figure 1 (left) presents the HAZUS (2017) combined earthquake and tsunami methodology, and the alternative method proposed in this study (Fig. 1,right). Note that in the right side of Fig. 1, we propose to make a computation loop regarding the exposure, to update it considering the spatial distribution and to consider its new vulnerability established after the earthquake damage. Hence, this computation loop is applied to each census block, to each typology included in the census block and to each damage degree reached by each typology after the first event. The area of focus of the study is the French West Indies territory.

Seismic hazard: source characteristics and intensity measures
Two earthquakes are considered: interplate earthquakes offshore Martinique and Guadeloupe (Fig. 2). Large historical events were reported in 1839 and 1843 offshore Martinique and Guadeloupe, respectively. Though they are generally associated with megathrust ruptures, no large tsunami were reported (e.g. Bernard andLambert 1988 or Hough 2013). Nevertheless, source hypothesis from Feuillet et al. (2011) for these events are inspired from macroseismic data and active faults analysis. Both earthquakes are localized offshore, in the Atlantic Ocean, to the East of Martinique and Guadeloupe Islands. The epicentre of the 1843 earthquake is located to the west of Guadeloupe and the 1839 earthquake is located to the west of Martinique (Fig. 2). As explained in the lines hereby, when we are speaking about tsunami events, the 1843 earthquake is referred as IOG and the 1839 earthquake is referred as IOM. Such scenarios can be considered as "classical" scenarios for major rupture of the Lesser Antilles megathrust. These two scenarios are chosen better than others that could better explain that lack of significant tsunami in 1839 and 1843 as they represent an ideal case study for modelling near-source earthquake and tsunami damages. Moreover, even if such source characteristics fit macroseismic data and active faults characteristics, but did not fit tsunami data, they remain realistic hypothesis for rupture of 1 3 the Lesser Antilles megathrust offshore Martinique and Guadeloupe. Two hypotheses are considered for the damage assessment: • A large mega thrust event offshore Guadeloupe from Feuillet et al.( 2011), referred as Interplate Offshore Guadeloupe (IOG) scenario; • A large thrust event offshore Martinique from Feuillet et al. (2011), referred as Interplate Offshore Martinique (IOM) scenario.
The localization of the two earthquakes is an important point for the interpretation of the results. Table 1 presents the parameters used for the generation of earthquakes.
The bedrock acceleration is calculated by a new developed python version of Armagedom software (Sedan et al. 2013 now integrated into VIGIRISKS platforme -Tellez-Arenas et al. 2019; Negulescu et al. 2019) using the fault position, the moment magnitude, the fault mechanism and the attenuation laws. The bedrock acceleration is majored with a soil coefficient to take into consideration the site effects that are quite important in this area. The site effects have been established during French Government projects that funded several years of H/V measurements and soil investigations first, for some communes (Bertil et al. 2009a;Roulle et al. 2010) and subsequently, for the entire territory (Sedan et al. 2008 for the methodology, Bertil et al. 2009b for Guadeloupe and Belvaux et al. 2013 for Martinique). Depending on the hazard intensity measure (IM) parameter used for the damage calculation (e.g. acceleration for HAZUS or EMS98 intensity for RISK-UE level 1 2001 -Milutinovic andTrendafiloski 2003), the acceleration can be converted to intensities. The calculation of the intensities for our simulations also provides the ability to check them with the historical information. Figure 3 compares the Armagedom simulation results in terms of intensities with that described by the SISFRANCE historical events database, which is an abundant and reliable database, and is the reference historical database in France (SISFRANCE, BRGM). Note that the range of intensities is quite similar for the Armagedom (Sedan et al. 2013) simulations and the historical data.

Tsunami hazard: propagation and submersion
To calculate the tsunami hazard parameter, which is the maximum moment flux, HAZUS uses the median value of the maximum inundation height, H (also called the water depth or flood depth, that is, the sea surface elevation minus the topography), and the median value of the maximum velocity of the water, v.
We use the inundation height grids and flow velocities grids calculated based on numerical simulations involving modelling the tsunami source, propagation and inundation. The tsunami hazard data include inundation grids provided from an external numerical tsunami hazard model (Le Roy et al. 2017;Poisson et al. 2009) that corresponds to a Level of Analysis 3 according to the HAZUS description. Tsunami generation is classically modelled with (Okada 1985) formulation. FUNWAVE-TVD (Shi et al. 2012) is applied in order to simulate tsunamis, solving the Boussinesq equations. Two chains of four nested ranks are prepared in order to be able to perform tsunami propagation and coastal flooding models for Guadeloupe and Martinique. Topo-bathymetric grids were built from different sources (Litto3D ©IGN-SHOM, HISTOLITT from SHOM and GEBCO). For finest rank, heterogeneous soil friction according to land use is used and obstacle to flood is integrated. In our study, the final grid mesh is in the range of 20 to 30 m. The inundation grids of the numerical tsunami hazard model give the height, H, and the velocity, v, of the water at each time step; thus, for each area there are 100 inundation grids by each intensity measure.
Generally, the maximum values of the height, H, and the velocity, v, of the water are not reached at the same step of the simulation in the same spatial cell. Hence, by looking separately to the evolution of the two parameters, the height, H, and the velocity, v, of the water, and by using their median values, some information could be lost and we could miss the highest moment flux (H*v 2 ), which is the most important parameter for us, since the flux damages the structure. Hence, for the calculation of the tsunami hazard we decided to calculate the moment flux (H*v 2 ) at each time step in each cell of the grid. Next, the hazard parameter that is compared with the building capacity is the maximum value of the moment flux for all time steps in each grid. Contrary to HAZUS, which uses the median value of the maximum height and velocity parameters, we use the real maximum moment flux, which is not reached at the same time in each pixel of the grid, as this clearly arrives in reality. Figure 4 presents the maximum sea surface elevation (SSE), around Pointe-à-Pitre in Guadeloupe and along Martinique coast for IOM-like and IOG-like events.

Earthquake and tsunami inventory, vulnerability characteristics and damage calculation
The inventory of existing buildings was performed first by the GIS approach that determined the census blocks that could have the same vulnerability and occupancy characteristics. These census blocks are named "area of homogeneous vulnerability". Second, field investigations by the teams of civil engineers of BRGM during a period of 5 years were conducted to validate or refine the "areas of homogeneous vulnerability", called out further on the census block. These studies and field investigations were funded by different projects of the French Government, the same projects that characterized the soil characteristics. Bertil et al. (2009b) for Guadeloupe and Belvaux et al. (2013) for Martinique have finally regrouped these several projects. For some areas, extensive field inventories were performed; for others, extrapolations of the obtained data to the other communes were performed. The typologies were classified considering the local construction mode and local vulnerability aspects. Hence, the assessment of the vulnerability is based on local field inventories that took weeks to months to be performed, which corresponded to Level 2 of the HAZUS methodology. The total number of census blocks is 2950 for Martinique and 2054 for Guadeloupe. The spatial repartition of these census blocks is presented in Fig. 5. For each census block, the following information is provided: typologies, range of construction years, number of houses, number of dwellings, surface, density and type of occupation. In Fig. 5, we represented the density of each area classified as a function of the number of houses per hectare. The earthquake vulnerability is governed first by the earthquake-resistant design for both European and American approaches, and second by the construction material (wood, RC, masonry). In addition to the construction material, the assessment of the construction mode (e.g. cast-in-place, prefabricated) is also important. Finally, thanks to these data acquired during the French programmes, the buildings are classified in the 36 model building types of the HAZUS Tsunami Model, which are the same as those of the HAZUS Earthquake Model (FEMA 2011). For a coherence between the studies, we kept the "local name typologies" as a first column of Table 2 (named "Type CDRS") and we observed the corresponding HAZUS Typologies and RISK-UE typologies. The double assessment of the vulnerability allows us to use, for the calculation of the buildings damages, either the RISK-UE or the HAZUS methodologies, or to use one

Earthquake and tsunami damage calculation
For the earthquake damage calculation, we can apply either the RISK-UE methodology, LM1 (empirical method) and LM2 (mechanical method) (Milutinovic and Trendafiloski 2003;Giovinazzi and Lagomarsino 2004), or the HAZUS methodology (FEMA 1999). Both of them are well known and are not detailed in this paper. Considering the historical database of intensities (SISFRANCE, BRGM), the field inventory data that we have and the experience acquired on the validation of the simulated damages with observed damages (Douglas et al. 2015;Boutaraa et al. 2018), we decided to apply RISK-UE methodology LM1 (the one based on relation between the intensity, the vulnerability index and beta damage distribution) for earthquake damage calculation. The decision to use RISK-UE methodology LM1 is also supported by studies made by others teams, for example Lestuzzi et al. (2016) that conclude that "There are significant differences in global results between LM1 (empirical method) and LM2 methods (mechanical method). The mechanical LM2 method is more pessimistic since it predicts damage grades of about one degree higher than LM1 method. However, the main drawback of the empirical LM1 method is that an a priori determination of an adequate value of the macroseismic intensity is required. Nevertheless, LM2 method may lead to a global overestimation of damage prediction". Regarding the tsunami damage calculation, the methodologies are less known, since they are quite recent. Several new approaches have been published lately (Attary et al. 2017;Petrone et al. 2017), and they generally use a Performance-Based Tsunami Engineering (PBTE) approach for tsunami damage assessments. This is also the approach adopted by HAZUS methodology (published in November 2017-FEMA 2017) and implemented in our study. Our research is focused on the structural damage; hence, we use the functions developed for tsunami "flow" hazard since the lateral forces due to tsunami flow are the primary cause of damage to the building structure, including building collapse. The damage due to tsunami inundation affects primarily the non-structural systems, which are not the object of our study and hence, the tsunami "flood" hazard is not considered in our scenario. More recently, Petrone et al. (2017) have compared tsunami pushovers to dynamic analysis in predicting structural responses, and they concluded that "tsunami peak force" is a more efficient intensity measure (IM) than flow velocity and inundation depth. This new parameter for tsunami intensity measure seems to have a good potential; however, its development is not sufficiently advanced nor sufficiently correlated with the other tsunami parameters (typologies, damage scale) to be integrated currently in general damage scenarios.
The damage on structures is due to the lateral forces caused by drag effects and debris carried along by the tsunami flow (Sect. 5.5 HAZUS Tsunami). The lateral force due to tsunami flow that acts on the structure is expressed by the following equation: where F TS is the tsunami force on a building, K d is a coefficient used to modify the basic hydrodynamic force, ρ s is the density of the debris flow materials, C d is the drag coefficient, B is the building plan dimension normal to the flow direction and hv 2 is the maximum momentum flux (flow depth times velocity squared). (1) The lateral capacity of a building is given by where F y and F u are the yield and ultimate force at the base of the building, respectively, α 1 is the modal mass parameter, A y and A u are the yield and ultimate spectral accelerations, respectively (Table 5.7 Hazus EQ Technical Manual), and W is the total weight of the building.
By equating the debris flow force (Eq. 1) with the lateral building capacity (Eq. 2), we obtain the moment flux parameter, hv 2 , for each damage state. The mean and median values of the building's strength, A u and A y , and their variabilities are taken from the HAZUS tables. The use of the lognormal damage probability functions is necessary to reflect the uncertainty. The metric for structural damage is the maximum flux as presented in Sect. 3.2. For the classification of typologies, for which we associated the A y and A u values, the building type and design level are important, as presented in Table 2.
At this stage, for the calculation of the tsunami damages, we do not consider the damages to structures due to the initial earthquake; in other words, all the structures have the initial vulnerability or are undamaged (damage state D0).
The influence of the debris of the collapsed buildings after the earthquake has not be accounted in Eq. 1. This issue is important and can affect the results but the estimation of this issue is complexe and require for an important study in itself, that it is not in the scope of this paper. The debris of the collapsed buildings after the earthquake can have an important impact on the tsunami flooding in the damaged area and hence the risk of the first event can affect not only the vulnerability of structure, but also the extend of the hazard area of the second event. The debris of the collapsed buildings can change the soil roughness taken into consideration the flood extent estimation trough the Manning coefficient. In this case, a better estimation of the tsunami flooding should be considered (up to high-resolution modelling) including an update of the Manning coefficient after earthquake damage calculation.

Earthquake and tsunami building damage states
For earthquakes, the two widely used damage scales in loss estimation studies are: (1) EMS-98 (Grunthal et al. 1998) an intensity (damage) scale developed primarily for postearthquake field investigations; and (2) the US HAZUS (FEMA 1999) damage scale, an earthquake engineering scales from FEMA. These scales are used in both their original form and in a "hybrid" form (named by Hill and Rossetto (2008) "hybrid scales"-scales that are adapted versions of original scales). The RISK-UE project proposes scales for two methodology levels (Milutinovic and Trendafiloski 2003): LM1, based on post-earthquake field investigations and, hence, intensity parameters; and LM2, based on physical parameter values calibrated to European construction that follows similar principles to that used in HAZUS. In our software, both methodologies can be used. However, there is not a unique damage scale that can be used without interpretation and adaptation in the US and European context. This subject is difficult; for more detail the authors recommend the Hill and Rossetto (2008) study that carried out a complete review of damage scales for buildings with a view to assessing their suitability for use in earthquake loss modelling in Europe. They concluded that none of the considered damage scales adequately satisfies all the criteria necessary for their use in European seismic loss estimations, and consequently, they propose equivalence tables, showing the relationship between the damage states of each (2) F y = a 1 A y W and F u = a 1 A u W considered scale. For earthquakes, the scales of damage that we are using in our simulations are RISK-UE (Milutinovic and Trendafiloski 2003) and US HAZUS (FEMA 1999).
For tsunamis, the building damage states are defined differently from one study to another and from one region to another. Among all the studies, three damage state definitions related to the tsunami as being the most adequate, the most used and therefore documented are (1) the Global Earthquake Model (GEM); (2) FEMA-HAZUS methodology; and (3) Japanese Ministry of Land Infrastructure Tourism and Transport. Eguchi et al. (2013) in their study, "HAZUS Tsunami Benchmarking, Validation and Calibration", using the data after the 2011 Japan tsunami, mapped the damage states used in Japan into categories that are used by the HAZUS methodology (Fig. 6).
As our area of study is the French West Indies and the building characteristics are closer to the American structures than those constructed in Japan, we use the damage state definitions given by the HAZUS methodology. As shown in Table 3, three non-nil damage states are used for the tsunami by HAZUS: Moderate, Extensive and Complete. These damage states are the same as those (of the same name) used by the HAZUS Earthquake Model to describe the extent and severity of damage due to ground shaking and ground failure. Although the specific cause and manifestation of tsunami damage can be quite different from that of an earthquake, tsunami and earthquake damage states are considered to be the same when they represent a common extent and severity of damage (HAZUS 2017, 5.1.1).
The probability of occurrence of the complete damage is the value read on the corresponding fragility curve. The probability of occurrence of the extensive damage is the value read on the corresponding fragility curve minus the value read on the fragility curve corresponding to the complete damage state.
For clarity, and to separate the damages that come from the two phenomena, in Sect. 5 we address the damages from earthquakes in terms of the D0-D5 damage degrees (Column 3 of Table 3) and the damages from tsunamis in terms of Moderate, Extensive and Complete (Column 6 of Table 3). The damage coming from the combined earthquake-tsunami event is addressed in terms of Slight, Moderate, Extensive and Complete.

Update exposed elements and their vulnerabilities considering successive earthquake and tsunami aggression
The differences discussed in this section compared to the previous section are that we do consider the likely damage to structures due to the first event, the earthquake. Using the first hazard event and the initial vulnerabilities, we calculate the damage due to the first event. Then, we update the capacity of resistance to the lateral force for each typology in each census block taking into consideration the number of buildings in each typology that reached each damage state (D0 to D5). Finally, we calculate the second event damages, based on the second hazard event and on the reduced vulnerability. The parameters that characterize the vulnerability of the structure were adapted (diminished) in order to consider the earthquake damage. In the calculations, we update two elements of the processes: (1) the number of the exposed element already damaged by the earthquake and (2) their vulnerabilities considering the impact of the earthquake. For the exposed elements, we assume that the buildings that reached the damage degrees D4 and D5 after the earthquake no longer exist as exposed elements for the tsunami hazard, which comes after the earthquake damages. Hence, the first update is the total number of buildings by typology and by census block, considering the damage degrees caused by the earthquake.
Considering the vulnerability adaptability for the structures that are no damage, D0, or in damage degree D1, the tsunami method is applied without changing the vulnerabilities. For the buildings in damage degrees D2 and D3, the initial method presented in Sect. 4.2 is modified to consider that the structure is not in its original state, but is already damaged. That means the modification at the level of the yield and ultimate spectral acceleration (A y and A u ) parameters. Figure 7 shows the initial capacity curve of the structure and the points green, yellow, orange and red that represent the spectral displacement when the corresponding damage state is reached. At the same time in Fig. 7, we represent in the colours green and orange the reduction in the rigidity of the structure, through the inclination/ angle of the elastic branch of the capacity curve, after reaching a certain damage degree. As an example, in the top right part of Fig. 7, we are showing a real case of four identical buildings for which we performed microtremor measurements and analysis after the 2007 Martinique earthquake. The measured frequencies follow the reduction in the rigidity of the structures due to the damages occurred; the damages are different from one structure to another. The red, orange and green colours of this table can be associated, for an exemplification aim, with the red, orange and green points of the capacity curve. In order to apply a general approach for all the structures, to compute the residual capacity of a mainshockdamaged building to withstand an aftershock (or another lateral solicitation as is the case SD is structural damage; N-SD is non-structural damage of the tsunami flow) we used the tool SPO2IDA (Vamvatsikos and Allin Cornell 2006). SPO2IDA provides a direct connection between the static pushover (SPO) curve and the results of incremental dynamic analysis (IDA). Therefore, for the purpose of calculating the A y and A u parameters for the building that already reached a certain limit state, LS i , we use a practical tool called SPO2IDA, which is capable of converting static pushover curves into 16%, 50% and 84% IDA curves using empirical relationships from a large database of incremental dynamic analysis results. Hence, using SPO2IDA, we recreate the seismic behaviour of oscillators with complex quadrilinear backbones and different rigidity inclines to simulate a building that has already reached a certain LS i . From Luco et al. (2004) in Fig. 8 (left), there is a representation of a pushover curve for an undamaged structure (D0); Fig. 8 (right) shows the "reduced" pushover curves considering the limit state, LS, reached by the structure after the first event.
Hence, from a development viewpoint, the initial file that describes the vulnerability is replaced by six files of vulnerability to have the new parameters that characterize the structure that reached each one of the damage states (D0 to D5) after the main shock. To implement this procedure in a more general process of a damage scenario, this procedure is integrated and implemented in the software for each typology and for each census block, as described in Fig. 9.
Of course, each of these two processes can be perfected. The aim here is neither to debate what reference to choose for one or another coefficient, nor to have an exhaustive view of these values, but to apply to the scenario process all these updates, and to evaluate for the chosen values the impact on the total damage. In this case, the aim is to apply the informatics and physical processes, and not to test one or another coefficient.

Earthquake scenario
The main work in this paper is related to the tsunami and cascade damage calculation, and hence, the descriptions of the damages due to earthquakes are not highly detailed. Naturally, a good estimation of the earthquake damage is essential for a reliable treatment of results thereafter. Hence, we here present the results for the entire islands (Figs.10,11) and for the zoomed areas, respectively, corresponding to the area that is most affected by the tsunami in each one of the two islands: Trinity in Martinique and Pointe-à-Pitre in Guadeloupe. Further attention is paid to these zoomed areas. Figure 10 presents the results for the 1839 earthquake in Martinique, where we note that the Eastern part of the island is more affected than the Western part. This is normal considering the epicentre position of the earthquake (Fig. 2, left) added to some digital elevation model (DEM) characteristics leading to local tsunami amplification (e.g. La Trinité to the west of Martinique). Hence, Fort-de-France is not much impacted by the earthquake and tsunami, contrary to the east of the Island where the percentages of the damage degrees D4 and D5 reached 30% in a census block. Hence, Fort-de-France is not much impacted by the earthquake and tsunami, contrary to the east of the Island where the percentages of the damage degrees D4 and D5 reached 30% in a census block. The high rate of damage coming from the earthquake event considerably affected the results on the combined cascade event. Figure 11 presents the earthquake damages in Guadeloupe for the 1843 earthquake, which was more damaged than in 1839 in this island. The percentage of damage degrees D4 and D5 can reach up to 50% in some census blocks. Even though the epicentral position is located in Northeast of Guadeloupe (Fig. 2,  right), the area of the Pointe-à-Pitre Bay is largely affected by both earthquake and tsunami events.
For the More Classical Scenario for 1843 (or IOG-like), Fig. 12 shows a bigger damage degree for Guadeloupe since the damage distribution is centred on D2 (Guadeloupe is closest from the epicentre of 1843 event), while for Martinique the damage distribution is centred on D0. Inversely, the IOM-like earthquake damage is more important in Martinique centred in D2 (closest to the 1839 epicentre), while for Guadeloupe the large majority of buildings are undamaged (damage degree D0). As outlined in Chapter 3, Seismic  and Tsunami Hazards, the historical observations of the intensities confirm the simulated seismic intensities (Fig. 3); hence, when we are speaking about the earthquake events the 1843 and 1839 earthquake can be mentioned like this, since they provoked these intensities. However, the historic data cannot sustain a hypothesis of a significant tsunami with the value modelled in our study; so, in this case we refer rather to IOG-like event than to the 1843 earthquake and, respectively, to IOM-like event than to the 1839 earthquake.

Tsunami scenario and cascade event earthquake-tsunami scenario
The difficulty of combining or comparing the damages coming from different events stems also from the difference of the scale of the impacted areas. The earthquakes, situated in the sea at East of either Guadeloupe or Martinique, affect both entire islands, in a greater or lesser degree of damage. Hence, in the case of earthquakes, the scale of the treatment of the damage is the scale of the entire Islands. The tsunami flooding maps show that for both earthquakes, the areas affected by the generated tsunamis are quite localized. For Guadeloupe, the area most affected by the tsunami coincides with the densest area in the Bay of Pointe-à-Pitre. Figure 13 shows the impact of the tsunami that affected the bay, described in terms of the intensity measure by the momentum flux. The amplification of the tsunami waves is certainly due to bathymetry and to the geographical configuration of the bay, with the coasts relatively close where tsunami wave rushes. In the legend of Fig. 13, we present the inundation height and the momentum flux, which take also in account the water velocity. This explains why, we can observe, for almost the same value of inundation height, different range of values for the moment flux represented at the scale of census blocks. For  Fig. 13, for an inundation height of around 3 m, we have one census block in the moment flux range between around 20 and 30, and the other neighbour census blocks in the moment flux range between around 10 to 20. For a better understanding, in order to discuss the results of tsunamis, we present here the figures zoomed in on the tsunami affected areas: Trinity in Martinique and Pointe-à-Pitre in Guadeloupe. To exemplify the difference of scale of the affected area by the earthquake and tsunami, Figs. 14 and 15 present the damages from the IOG-like tsunami and the damages from the IOG-like earthquake at the area of Pointe-à-Pitre Bay, respectively.
Note that the area impacted by the tsunami is contiguous to the coast line, while for the earthquake, it is spread throughout the entire island. This is an obvious observation from a phenomenological viewpoint; however, this observation alerts us as to how combining the damages coming from the cascade events that do not have the same impact scale. Figure 16 quantifies this localization of damage for the tsunami by comparing the damages for all islands and for the impacted areas by tsunami for the IOG-like tsunami. The bottom illustrations of Fig. 16 show that there is practically no difference of damage between all of Guadeloupe and the Bay of Pointe-à-Pitre, and hence, all the tsunami damage is concentrated in a very small territory. Concerning Martinique, there is a more important difference in the damages for all of Martinique and the Trinity damages, especially for the complete damage, which indicates that some other costal communes were affected by the IOG-like tsunami. Figure 17 shows the distribution of damages for the IOM-like tsunami. Guadeloupe has not been impacted (see Fig. 17 right). Concerning Martinique, as for the IOG-like tsunami, we note that a large amount of damage is concentrated in Trinity. In Trinity, for the IOM-like and IOG-like tsunamis, we observe almost the same damage distribution; certainly, the same area is impacted. For the entire Martinique Island, the number of damaged Keeping in mind Fig. 14 (only tsunami event) and 15 (only earthquake event), Fig. 18 compares to themes, the damage at the scale of Pointe-à-Pitre coming from the combined cascade event-the earthquake and tsunami. We note that the majority of the areas affected by the tsunami in one category of damage (e.g. 0-10%, 10-35%) make a jump in the higher category of damage when we consider the cascade event. This is the case for the central area of Fig. 18, already affected by the tsunami and now almost completely damaged after the cascade event. For the bottom areas, right and left, of Fig. 18, we also note that some areas that seem to not have complete damage state neither after tsunami nor after the earthquake, present finally complete damage state after the cascade event.
Zooming in on Trinity, Martinique, Figs. 19, 20 and 21 show the damages due to the earthquake, tsunami and combined event earthquake, and tsunami for the IOM-like event, respectively. Figure 19 presents the earthquake damage and presents a maximum percentage of complete damage of maximum 30% per census block in a large affected area. In contrast, Fig. 20 presents the tsunami damage and presents a maximum percentage of complete damage going until 100% per census block; however, a smaller area is affected. Figure 19 and Fig. 20 focus on the difference of damages when considering only a tsunami or a combined earthquake and tsunami event. Certain areas on the coast of Trinity, affected by only tsunamis, are in a more important category of damage when considering a combined earthquake and tsunami event. The most impressive change is the larger impacted area in the centre of Trinity, due to fact that the buildings already damaged by the earthquake are

Discussion
Tables 4 and 5 present the damage results for Martinique and Guadeloupe, respectively, for the IOG-like and IOM-like earthquakes and for the three phenomena sequence approaches: earthquake only, tsunami only and combined cascade events earthquake and tsunami. The earthquake damages are expressed in terms of damage degrees, from D0 to D5, and the tsunami and combined event damages are expressed in terms of Slight, Moderate, Extensive and Complete. Table 3 makes the connection between the definition of damage levels following different methodologies. However, the damage state classification is a complicated and tricky issue that can have a big impact on the results. As discussed in Introduction and in Chapter 4.2, we make some hypotheses of the damage states, but future studies done by other researchers for better characterization, and connection between the damage states, can be easily adopted in our methodology. As the estimation of damages is based on a statistical methodology, the numbers in Tables 4 and 5 are rounded to the nearest decimal.
In Martinique, we inventoried a total number of 2950 census blocks that cover the entire island. For this total number of census blocks, only 97 areas are affected by the tsunami induced by an event IOM-like earthquake and only 74 for a tsunami induced by IOG-like earthquake. The affected percentage of census block is around 3%, but the extensive and complete damages increasing much more, arriving until almost doubling in certain combination of event and exposed territory. The same range of values, related to the contrast between the small territory affected by the tsunami and the high increasing of damages for the combined event, are observed for the Guadeloupe Island. For a tsunami induced by IOG-like earthquake, 150 census blocks are impacted on a total of 2054 census blocs. However, the hypotheses used for a tsunami induced by IOM-like earthquake affect only 6 census blocks in Guadeloupe. This is a very insignificant number and is completely in the error range of the methodology; hence, the zoom in Guadeloupe on IOM-like event is not presented in Table 5.
When we study the occurrence of two natural phenomenon, a tricky issue is to evaluate the temporality as well as the spatiality of the two events. As described above, the area affected by the tsunami is, of course, the costal territory, and it is comparatively much smaller that the area affected by the ground motion. The first issue is to correctly identify the area, and the exposed elements affected by the two events. The correct estimation of this area, and of the number of the affected elements, is essential for a good application of the relations used to combine the damages coming from the both events and can reduce the bias on the results. To quantify the effect of the increasing of the vulnerability due to the damages of the first event, we compare the damage results only in the areas affected by the tsunami, since they are a subject of the vulnerability changes. Table 6, and Figs. 22 and 24 present the results only for the area affected by the tsunami. Further, on this paper, we are paying a particular interest to the differences of the results expressed in damaged buildings in this area, due to the different applied approaches.
Contrary to IOM-like tsunami, the IOG-like tsunami, for which the epicentre is situated in offshore of Guadeloupe, affects the both islands. In Martinique, the damages caused only by tsunami are not so different from the damages caused by tsunami and earthquake combined events (Fig. 22, left). This can have two explanations: either the water height is shallow, hence there is not a sur-damage due to tsunami event, or the earthquake did not damaged the structures, and hence, the combined event can be assimilated to an only tsunami event (the only event that has a sufficient intensity to provoke damages).
In order to help to discriminate between the two possibilities mentioned above, we are plotting, in Fig. 23, the moment flux values as a function of the probability of exceedance of a certain damage degree. It can be notice that the maximum MMF for IOG is in the range of 25 m*(m/s) 2 and around 40 m*(m/s) 2 for IOM. For IOG, the damage state the most impacted is the Moderate State for a value of MMF beneath 10 m*(m/s) 2 , while for IOM the three damage states are impacted by the manner and order of treating the events (Moderate, Extensive and Complete). The inversion of the probabilities, between the tsunamis only event and cascade events earthquake and tsunamis, is more apparent for the IOM. For IOM, in between the range from 20 to 30 m*(m/s) 2 we have systematically higher percentage of tsunami damages in Moderate state and, in consequence, higher percentage of damage for the Complete damage state, due to conjoint events, tsunami and earthquake.    In Guadeloupe, the damages coming from the combined events are much more important than the damages coming only from tsunami event (Fig. 24). These should be associated with Figs. 14, 15 and 18, which localize, at the scale of the census block, the range of percentage of damages coming from different sequences of physical events and with Fig. 25 which presents the moment flux values as a function of the probability of exceedance of a certain damage degree, in Guadeloupe.
Using the damage results obtained from separate calculation coming from earthquake and coming from tsunami, we calculated the probabilities of damage to the structures due to the combine earthquake and tsunami hazards using HAZUS (2017) formulas from 5.18 to 5.21. These probabilities of damage from combined hazards calculated using HAZUS formulations are compared with the values calculated numerically in this study for the cascade event simulation, with the degradation of the vulnerability of the exposed elements affected by the first event. The comparison between the probabilities of damage to the structures using HAZUS and this study is presented in Figs. 26 and 27 for Trinite, the zoomed area of Martinique for IOM and, respectively, IOG events.
We can observe that in our study the probabilities of being in Complete damage state are systematically superior to the same probabilities calculated with the HAZUS Formulas. For the Extensive damage state, our study gives less probabilities that the HAZUS Formulas. The results for Slight and Moderate damage states are not so different between the two approaches, which means that the reduction in the capacity curves do not affect significantly these damage states. Also, the equivalence of the damage scales (Table 3) can induce a bias to this issue for the lower damage state (Slight and Moderate), but this was discussed previously. We can also see that the more the earthquake event is damageable (before the impact of the tsunami), the more the differences of the results coming from the two approaches are higher. In Fig. 26, for IOM event, the Complete damage state for HAZUS approach is equal to 5%, and in our study, the Complete damage state is around 13%. This difference is vanishing in Fig. 27, IOG event, where of the both approaches the Complete damage state is around 10%. In both situations, we remember that you are treating the same territory with the same approaches, and only the characteristic of the seismic event are changing. In both cases, the tsunami event is almost the same damageable, but for IOM case, the earthquake damage is more important that for the IOG event, which may cause the differences in the combined results. This observation is also verifiable in the case of Guadeloupe (Fig. 28), where the damages from earthquake event are important and the difference in the Complete damage state from the both approaches is of around 10%. As discussed in the previous sections, the combination of the damages coming from a tsunami and earthquake is quite complicated, since the scale of the impacted area and the damage description are not of the same order. Some formulas for combining the results from the separated hazard calculations are proposed in the literature; however, the results obtained in this paper show that the number of damages is closely related to a specific area, its topography, cartography and its building physical characteristics, and it can be very difficult to generalize this combination formula to others areas. Moreover, for the same affected area, but different characteristics of the seismic event, and hence of the tsunami impact, the combination of the results can be difficult. This is the case for the Trinité area affected by IOM and IOG events. For the IOM event, the results between the two approaches are quite similar, while for the IOG event, in exactly the same area, the results show significant differences especially for the Complete damage state. In addition, in our study, we can noticed that the amount of important damages caused by the earthquake event, may be the cause of the differences between the results, but this should be confirmed by more studies in other geographic areas.

Conclusions
The relationship between the number of damaged structures and earthquakes alone and earthquake-and-tsunami cascade events was investigated in this paper from informatics, practical and realistic viewpoints. The database for buildings and soil responses are based on field investigation campaigns performed during several projects funded by the French Government. Several practical novelties are applied to the combined earthquake-tsunami damage scenario. First, the acquisition of each intensity measure for all the time of the event is primordial, since the maximum sea surface elevation and the maximum water velocity are not necessary in the same time in the same spatial cells. Hence, this information allows us to stick as close to the phenomenon as possible, instead of using median values of the intensity measure parameters. Second, for the cascade scenario, we considered the degradation of the initial vulnerability of a damaged structure, expressed here by a capacity curve for tsunami calculation, considering the reached limit state by the structure after the first event (i.e. the earthquake). This step complicates the direct hazard-vulnerability-damage calculation, since the exposed elements are updated in the number of structures but also consider the new degraded vulnerability, which implies an informatics loop and a multiplication of the vulnerability characteristic equal to the number of the damage degree used to describe the first event. From a general point of view, this procedure can be applied to other areas of study but also to the other combinations of the two natural 1 3 Fig. 25 Relation between the moment flux (H*v ²) and the Damage States (Moderate, Extensive, and Complete) calculated for tsunami event and for combined tsunami and earthquake events, for Guadeloupe Island phenomena. This step has the merit of overcoming a complicated stage of the combination of damage coming from two different hazards, as is currently used. From a specific point of view, the combined scenario provides an opportunity to observe real damages and to quantify them stemming from three different approaches: earthquake only, tsunami only and combined earthquake and tsunami event. This quantified result highlights the complexity of treating two phenomena that have a difference in their scale of impact and of intensity of damage: a large area of impact for an earthquake, and graduating damage and a considerably more localized area of impact for a tsunami, but heavy damages. However, the correspondences between the RISK-UE and the HAZUS typologies (Table 2) and between the RISK-UE and the HAZUS damage degrees (Table 3) are based on experts choice and these can introduce some biases in the results, but even with others choice of the correspondence between these two parameters, the calculation method remains reproducible and stable. Future work will focus on a better integration of the tsunami flooding in the damaged area (up to high-resolution modelling) including an update of the Manning coefficient during Complete (C) damage states resulting from (i) earthquake event, (ii) tsunami event, (iii) HAZUS approach for combining the damages coming from earthquake and tsunami and (iv) this study approach for combining the damages coming from earthquake and tsunami damage calculation. Another important issue to be addressed in the future studies is the influence of the debris of the collapsed buildings on the force applied to the structure but also on the changing of manning roughness coefficient due to the debris. and Complete (C) damage states resulting from (i) earthquake event, (ii) tsunami event, (iii) HAZUS approach for combining the damages coming from earthquake and tsunami and (iv) this study approach for combining the damages coming from earthquake and tsunami