Ex post analysis of engineered tsunami mitigation measures in the town of Dichato, Chile

Due to Chile’s notorious and frequent seismic activity, earthquake- and tsunami-related studies have become a priority in the interest of developing effective countermeasures to mitigate their impacts and to improve the country’s resilience. Mitigation measures are key to accomplish these objectives. Therefore, this investigation adopts a tsunami damage assessment framework to evaluate the direct benefits of tsunami mitigation works implemented by the Chilean government in the town of Dichato in the aftermath of the 2010 tsunami. We perform an ex post analysis of the potential damage reduction produced by these works studying what would have been the consequences on the built environment if they were in place for the tsunami that hit this area after the Maule earthquake in February 27, 2010. We use state-of-the-art tsunami simulation models at high resolution to assess the reduction in tsunami intensity measures, which serve as input to evaluate the benefit from averted damage against the costs of the mitigation measures. The obtained results show a reduction in the flooded area and a delay in the arrival times for the first smaller tsunami waves, but a negligible damage reduction when confronted to the largest waves. In conclusion, the tsunami mitigation measures would not have been effective to reduce the impact of the tsunami generated by the Maule earthquake in the town of Dichato, but could have had a benefit in retarding the inundation of low-land areas for the first smaller tsunami waves. The latter suggests that these works might be useful to mitigate storm waves or tsunamis of much smaller scales than the one that hit central-south Chile in 2010.


Introduction
Chile is a highly seismic country, mainly due to the subduction process of the Nazca Plate underneath the South American plate (Lomnitz 1970;Araya 2007;Fritz et al. 2011;Centro Sismológico Nacional 2013), where the trench is located very close to the coastline. Many megathrust earthquakes have generated tsunamis, with very short arrival times. Over the past 10 years, four large-magnitude events (Mw > 7.5) have occurred in the country, which generated tsunamis with varying levels of destructiveness. These include the Mw 8.8 Maule earthquake (Fritz et al. 2011) in 2010, the Mw 8.2 Pisagua earthquake (Catalán et al. 2015) in 2014, the Mw 8.3 Illapel earthquake (Aránguiz et al. 2016) in 2015, and the Mw 7.6 Melinka earthquake (Xu 2017) in 2016. The 2010 event caused large-scale economical, structural, and human impacts in central and southern Chile (Gobierno Regional del Bío-Bío 2010; Fritz et al. 2011). This event also triggered an increase in basic and applied research in the country. In addition, the Chilean government invested in several tsunami mitigation projects and reconstruction plans (Siembieda et al. 2012;Herrmann Lunecke 2015;Khew et al. 2015), including hard structures.
Among the most damaged coastal towns in 2010 was Dichato, a relatively small settlement with three thousand inhabitants (according to the last census) (Koshimura et al. 2010), located north of the epicenter, that suffered the destruction of houses and infrastructures near the beach and estuary (Koshimura et al. 2010;Yamazaki et al. 2010;Martínez et al. 2016). After the event, the government proposed a tsunami mitigation plan, which included hard structural countermeasures to reduce the impact of tsunamis. The setting of this new infrastructure provides the opportunity to compare the effect of the tsunami with and without the presence of such mitigation works.
The significant impacts associated with recent tsunamigenic events worldwide have motivated the development of numerous studies devoted to hazard assessment and tsunami modeling (e.g., Arcas and Titov 2006;González et al. 2009;Montenegro-Romero and Peña-Cortés 2010;Nandasena et al. 2012;Imamura et al. 2012;Suppasri et al. 2013;Melgar and Bock 2013;Adriano et al. 2014Adriano et al. , 2016Ozer Sozdinler et al. 2015;Catalán et al. 2015;Goda et al. 2015;Park and Cox 2016;Santa María et al. 2016;Aránguiz et al. 2016;Martínez et al. 2016). For example, the Great Eastern Japan Earthquake and tsunami offered the opportunity to evaluate different types of mitigation measures. Nandasena et al. (2012) made an ex post assessment of the effectiveness of tsunami mitigation works, including vegetated dunes, coastal forest, and seawalls. These are tested using numerical modeling to determine how they influence inundation extent, and other hydrodynamic properties such as momentum flux and velocities (Nandasena et al. 2012). Adriano et al. (2016) evaluated the effect of breakwaters in the coast through similar tsunami modeling approaches, concluding that the breakwater present in Onagawa at the time of the tsunami attack reduced the tsunami impact by diminishing the maximum inundation depth in 2 m, even if it was completely destroyed by the flow (Adriano et al. 2016). The performance of a breakwater located in Kamaishi is assessed by Ozer Sozdinler et al. (2015), who state that even though the presence of a breakwater may generate higher flow velocities, these hard countermeasures were beneficial in decreasing inundation depths and retarding the arrival of peak tsunami inundation. Moreover, it was found that a damaged breakwater could still provide some degree of protection.
In Chile, no studies on the effectiveness of tsunami mitigation works have been published to date, but some efforts have produced data and relevant information that are worth mentioning. Santa María et al. (2016) produced an exposure model for residential structures that classifies them by characteristics such as construction materials, geo-localization and estimated replacement costs. Martínez et al. (2016) presented an assessment of the vulnerability of Dichato considering physical, socioeconomic, and educational dimensions, comparing pre-and post-event conditions. It is claimed that vulnerability remains large in this town, despite the existence of new tsunami mitigation measures.
In the present work, we aim to assess the benefits that the presence of new mitigation works in Dichato may have generated, if they had been present during the tsunami of February 27, 2010. To this end, we implement a methodology to estimate the direct damage from the tsunami inundation and compare its results with and without the new infrastructure in place. The methodology used is intended to provide ex post evaluation, for the first time in Chile, on the tsunami impact reduction that such works may produce.
The manuscript is organized as follows: Sect. 2 shows a brief overview of the important definitions and concepts being considered here, along with a characterization of the case under study. Section 3 describes the methodology, while in Sects. 4 and 5, we present the main results of the study and a discussion on their potential implications. The conclusions are summarized in Sect. 6.

Risk, hazard, exposure, and vulnerability
In what follows, risk is considered to be composed of three factors: hazard, exposure, and vulnerability (Weichselgartner 2001;González et al. 2009;Venegas San Martín 2012). The hazard is quantified by intensity variables that reflect the capability of a certain geophysical phenomenon to inflict harm (González et al. 2009). While in the risk analysis framework it is necessary to assess all the possible events that may occur in the area of interest to yield a probabilistic hazard assessment (Cutter et al. 2000;Pelling et al. 2004), here we limit ourselves to a scenario-based analysis using the actual 2010 tsunami, following similar lines as Suppasri et al. (2013).
The characterization of exposure requires an inventory of all the elements that could be affected by the hazard, such as communities and physical infrastructure (Pelling et al. 2004), identifying those located in the flooded area by the tsunami (Penning-Rowsell et al. 2005). Lastly, in the context of this study, vulnerability is accounted for structural elements only (i.e., buildings), by means of fragility curves which allow to link hazard intensity variables with the damage probability. Linking hazard, exposure, and vulnerability provides means of assessing the risk level and gives the opportunity to quantify damage reductions from the benefits associated with mitigation measures (e.g., urban planning and relocation, engineered mitigation measures, educational programs, etc.).

Estimation of the benefit from averted damage
The damage reduction of the mitigation works in Dichato can be quantified using a cost-benefit analysis that considers the investment of such works and the benefits they produce against a tsunami. It must be noted that these benefits are exclusively evaluated over direct damage assessment, as other less tangible benefits are beyond the scope of this research. The benefit of the project (B) is estimated as (Ministerio de Desarrollo Social 2013): where P(i) is the probability of occurrence of the event i , C NP refers to the costs with no project in event i , and C WP are the costs with project in event i . Since a scenario-based analysis is used in this study, the consideration of the probability of its occurrence is beyond the scope of the investigation. Therefore, Eq. (1) will be adapted, as explained in Sect. 3.2.3.

Estimation of damage
A tsunami fragility function or fragility curve constitutes a direct relation between a hydrodynamic feature of the tsunami flow and the probability of damage of a structure (Koshimura et al. 2009a, b;Mas et al. 2012;Suppasri et al. 2012;Adriano et al. 2014;Favier et al. 2014a;Urra Espinoza 2015). They are usually built empirically from different data sources (Koshimura et al. 2009a). As a result, the level of refinement in the data and its categorization (such as construction material and number of stories of the structure) may affect the outcome of the analysis (Goda et al. 2015). Despite this potential shortcoming, fragility functions are essential in the methodology to estimate tsunami damage. The characteristics of the ones used here are described in Sect. 3.2.3.

Description of Dichato and the local impact of the 2010 tsunami
Dichato is located at the south end of the Coliumo Bay, almost 40 km north of Concepción (see Fig. 1). It consists of over 3000 inhabitants according to the last official census in 2002 (Instituto Nacional de Estadísticas) and about 2000 homes. The most common construction materials of buildings are wood and masonry (Servicio de Impuestos Internos). The Coliumo Bay has low wave energy, owing to the protection provided by the Tumbes Peninsula to SW incoming waves. On the other hand, its horseshoe shape may induce long wave resonance (Gobierno Regional del Bío-Bío 2010). The 2.4-km-long beach of Dichato makes the area a touristic attraction, especially during summer, when the population can increase up to 5000 inhabitants (Koshimura et al. 2010;Venegas San Martín 2012).
At 3:34 a.m. (local time) in February 27, 2010, a violent earthquake shook the central region of Chile. The epicenter was located off the Chilean coast, about 105 km northeast from the city of Concepción (Saito et al. 2010;Yamazaki et al. 2010) with a reported moment magnitude (Mw) of 8.8 (Pulido et al. 2010;Saito et al. 2010;Mas et al. 2012;Robertson et al. 2012). The tsunami generated free surface disturbances observed all around the Pacific Ocean basin, including the coast of Japan (Koshimura et al. 2010;Robertson et al. 2012). Five hundred and twenty-one people died as a direct result of either the earthquake or tsunami. The tsunami itself claimed 124 victims and 46 missing, mainly in the coastal regions between latitudes 34.5° and 38° south (Nahuelpan López and Varas Insunza 2010;Fritz et al. 2011). Nearly 370,000 homes were damaged and the economic loss was estimated in USD 30 billion, which was roughly 15% of Chile's gross domestic product in 2010 (Yamazaki et al. 2010).
Dichato was one of the most affected locations. Inundation depths in Dichato were estimated by post-tsunami surveys to be around 8 m, and water penetrated as far as 1.3 km inland (Martínez et al. 2016). Though most people were able to evacuate, there were 66 fatalities in Dichato, mostly tourists and elders who underestimated the intensity of the phenomenon ). More than 1200 families reported damages to their properties (Koshimura et al. 2010), and approximately 80% of the exposed built structures were washed away (Gobierno Regional del Bío-Bío 2010), evidencing the high degree of destruction from the tsunami (Venegas San Martín 2012; Martínez et al. 2016). This level of damage can be explained by the existence of wood houses, which were unable to uphold the large inundation depths and velocities. Damage was also induced by floating debris, including loose boats that collided against structures.

Mitigation works and relocation
Due to the tsunami impact in 2010, the Chilean government financed a mitigation project aimed to provide protection for Dichato and the community against future tsunamis (Martínez et al. 2016). This mitigation project consists of a low-height seawall approximately 4 m tall and 1.7 km long, made of reinforced concrete along the coastline, and a channeling of the Dichato Estuary at its mouth (see Fig. 1c). With the information provided by the Division of Port Infrastructure of the Minister of Public Works (Dirección de Obras Portuarias-DOP) and reference prices from a catalogue (Portal Ondac Construcción), the estimated construction costs of the wall and channel are US$ 6.8 million (details of this estimation are given in "Appendix 1," Fig. 8 and Tables 4,5,6,7,8). On the other hand, five neighborhoods were relocated either to higher areas or to piloted houses north of the Bay. Building restrictions were established in the town, allowing only commerce-related infrastructure along the coastline of the town. In this area, there is also a mitigation park, consisting of trees planted along the coastline, but its construction only started in 2015 (Martínez et al. 2016).

Methodology
The tsunami hazard associated with the 2010 tsunami is estimated by means of numerical modeling. Exposure is assessed from building inventories, and the vulnerability is estimated through fragility curves.
The benefit of tsunami damage reduction from engineered mitigation measures on the coastal urban area is evaluated by a cost-benefit analysis of direct costs associated with the inflicted damage to buildings by the tsunami, and the investment costs for the mitigation measures. Hence, by taking the difference between results considering the situation with and without mitigation works, we can estimate the potential improvement that the project would have produced. We are not attempting to give a comprehensive economic evaluation of the future benefits, but to provide a quantification of the potential improvement that this works could produce if the town was subjected to a similar event on a ceteris paribus framework.

Tsunami modeling
We employ the open source code GeoClaw to perform the tsunami modeling. GeoClaw solves the nonlinear shallow water equations using finite volume methods and adaptive mesh refinement. The code has been used to model several historical tsunamis using bathymetric and topographic data (MacInnes et al. 2013;Melgar and Bock 2013;Arcos and LeVeque 2015). Tsunami waves along the coast of Chile have been known to last several hours (Catalán et al. 2015), including the 2010 event (Venegas San Martín 2012); hence, the simulation time is set to 8 h, which is considered sufficient to model the largest latearriving waves (Yamazaki and Cheung 2011).
The information available did not allow us to include buildings in the topography, so a spatially constant Manning roughness coefficient of 0.2 s/m 1/3 is considered, as it produces the best agreement with in situ measurements (see Sect. 4.1). The spatial grid for the tsunami simulation is composed of eight rectangular nested grids, with a threefold increase in resolution between each one. The smallest computational grid, enclosing the Coliumo Bay, is built from bathymetric data collected before the event (2009) and has a resolution of 2 m (see Fig. 1c), while the spatial discretization in the largest grid is 1 km. The time steps in the simulations are adjusted dynamically to enforce the fulfillment of the Courant-Friedrichs-Lewy (CFL) stability conditions. The reference level used for the bathymetry and topography was the mean sea level, corrected by the tide level at the time of the event.
To incorporate the mitigation measures in the tsunami modeling, an in situ campaign was carried out to obtain the updated topography with high resolution (wall, park, channels) using differential GPS. Data were corrected by means of the Ntrip network. As a result, the topography for the tsunami simulation considered the aggregate of the data collected before and after the event. The difference between the digital elevation model (DEM) including mitigation works and the original can be appreciated in Fig. 1c. The most important changes are around the wall and along the estuary.
The finite fault seismic source model proposed by Delouis et al. (2010) is used and validated. The initial condition of the displacement of the ocean's free surface is obtained from the Okada formulation (Okada 1985). Figure 1b shows this displacement surface where it can be seen that Dichato is located slightly north of one of the two maxima.

Exposed buildings
The exposed residential buildings in Dichato and the tsunami effect on them are considered in this study. For simplicity, only the physical damage induced by the tsunami waves to these buildings is considered. The social impacts or other economic losses are beyond the scope of this work because of the difficulty in quantifying their vulnerability, especially when it comes to social vulnerability (Cutter et al. 2003). The characteristics of houses, construction materials, price values, and year of construction were estimated from a database provided by the Revenue Service (Servicio de Impuestos Internos). This database consists of 1878 houses; however, during the fieldwork, only 452 of them could be properly identified and geo-localized of the database. Hence, we use herein only the identified buildings which are considered as a sample of nearly 40% of the houses directly exposed to the tsunami waves. The characteristics and the spatial distribution of the inventory can be seen in Fig. 2. Most of the houses in Dichato are built with wood and masonry, while the rest are built with concrete, steel or other less common materials (see Fig. 2b). On the other hand, the majority of houses were valued in approximately US$7000 (Servicio de Impuestos Internos) (see Fig. 2a).

Tsunami intensity measures
The hazard intensity is estimated by computing inundation depths, velocity, energy, and arrival times (Macabuag et al. 2016;Park and Cox 2016). Maximum inundation depths are considered as the most relevant parameter because when used in conjunction with fragility curves, they allow to evaluate the probability of damage of individual buildings (Cançado et al. 2008). Nevertheless, the other intensity measures are also of interest. Flow velocities, for instance, have been studied and utilized for the development of fragility curves for structural damage in buildings Park and Cox 2016;De Risi et al. 2017), as they may contribute with information about the interactions between the flow and the topography that otherwise cannot be attained Arcos and LeVeque 2015). In addition, both maximum inundation depths and velocities have been found to be applicable for the estimation of building damage due to tsunamis (Park and Cox 2016). Thus, maximum flow velocities and maximum energy-as an extension of the latter-are also estimated from the tsunami simulation as reference.

Damage estimation
Fragility curves are used to estimate the probability of damage generated by the tsunami to the physical infrastructures. Direct damage is thus considered to be the result of the interaction between hydrodynamic forces and the structures (Penning-Rowsell et al. 2005). The damage depends on the construction material of buildings, and several damage levels can be considered.  Suppasri et al. (2012) performed an analysis of all available data delivered by the Ministry of Land, Infrastructure and Transportation of Japan (MLIT) regarding the 2011 Great East Japan tsunami impact zone. These authors developed fragility curves with six damage levels for wood, masonry, reinforced concrete, and steel frame houses, with the damage levels shown in Table 1  ). These fragility curves were chosen since they involved the same four main construction materials observed in Dichato, with the possibility of including different damage levels for each materiality.
Considering that a common type of building in Dichato is best represented as a concrete house in the first story (2 m in height), and as a wooden house in the second story, an additional fragility curve was arranged for the purpose of this investigation, by combining those of wooden and concrete fragility curves. More details are provided in "Appendix 2.1." With the probability of having the different damage levels for each building, the expected loss is quantified in monetary terms. The cost of repair (C r ) is assumed to be a percentage of the construction cost (C t ), as the more expensive the construction cost of a house would mean a more costly reparation as well. Likewise, if a building is constructed with cheap materials, the costs of repairing are assumed to be low. The relation between the two costs is a function of the damage level, which determines how much of the total construction cost the repairing cost is. Table 1 shows this relation. For example, if the building suffers minor damage, the repair cost is a fraction of the construction cost, assumed as C r = 0.2C t . If the house is washed away, the repair cost is larger than the construction cost (C r = 1.2C t ). This estimated larger cost is attributed to the management and disposal of debris.
Finally, the total costs from the loss induced by the tsunami for each scenario are estimated as the aggregate of costs from all buildings. This is simply the sum of the product between the cost ratios and their probability of occurrence for every damage level and every type of material, i.e., where ds i stands for damage state i , c i for the corresponding cost ratio Cr/Ct for the damage state i , and P ds i for the probability of damage for damage state i . Finally, the expected damage ∑ i P � ds i � * c i for each building k is summed up to give the total expected damage. The cost-benefit analysis considers the comparison between costs that are averted by the application of mitigation measures and those without mitigation (Penning-Rowsell et al. 2005;Ministerio de Desarrollo Social 2013;Iwata et al. 2014). Equation (1) is adapted to take into consideration only one realization:

2010 Tsunami impact without engineered mitigation measures
First, the tsunami modeling is carried out without mitigation works in order to estimate a reference scenario, and to validate the model. For the latter, results are compared to field data from in situ surveys summarizing observations from five campaigns, resulting in a set of 37 measurements (Imamura et al. 2010;Matsutomi et al. 2010;Fritz et al. 2011;Mikami et al. 2011).
In 2002, the Japan Society of Civil Engineers (JSCE) proposed that K and values in the range 0.95 < K < 1.05 and < 1.45 indicate a good estimation (Tsunami Evaluation Subcommittee and Nuclear Civil Engineering Committee 2002).
The observed and numerical sets of inundation heights are used to evaluate the K i coefficients from five post-tsunami surveys, resulting in a set of 37 observations (Imamura et al. 2010;Matsutomi et al. 2010;Fritz et al. 2011;Mikami et al. 2011). These 37 measurements are then compared with simulated data to obtain K i , and by the use of Eqs. (4) and (5), the values of the parameter K and are found. The Dichato simulation yields K = 0.97 and = 1.38 , which is within the recommended limits.
The tsunami simulation results are also qualitatively compared with inundation maps obtained by Mas et al. (2012) (Fig. 3a). Relative sea levels during the event are compared with the information registered at the Talcahuano port inside the Bay of Concepción, as seen in Fig. 3b and with data provided by a Deep-ocean Assessment and Reporting of Tsunamis (DART) Station off the coast of northern Chile and southern Peru (Fig. 3c).
Peak values of inundation depths, flow velocities, energy, and arrival times estimated with the tsunami model are shown in Fig. 4. The highest inundation depth was nearly 10.8 m. The maximum run-up reached 16 m. For flow velocities, the maximum value, as shown in Fig. 4b, is nearly 11.4 m/s. Local energy values are calculated using Bernoulli's expression H = h + v 2 2g , where h stands for inundation depth, v for flow depth-averaged velocity, g for the gravitational acceleration, and H for the hydraulic head. The results obtained for maximum energy are shown in Fig. 4c. Figure 4d also shows that there is an inundation of low areas within 50 min, but that the larger inundation extent occurs after 100 min. The total costs for this scenario are of nearly USD $ 2.9 million.

2010 Tsunami impact considering the engineered mitigation measures in place
Next, simulations including the mitigation works are carried out considering the updated topography. Figure 5 shows the differences between this modified scenario and the baseline scenario. Negative values denote the baseline scenario values exceed those of the modified case. Values for inundation depths and run-ups with the incorporation of mitigation works have maxima of 10.6 m and 16 m, respectively, which are very similar to the base scenario. Near the estuary, depths are larger owing to the channeling effect of the works. However, for the maximum hydraulic head, the modeled values in the modified situation are lower than the baseline scenario, which highlights the potential in reducing the energy of the flow. The total costs estimated for the direct damage in this case are $USD 2.7 million. A summary of these results is presented in Table 2. Considering inundation depths are similar and that the damage estimation depends directly and solely on this variable, the benefit obtained with the presence of the mitigation works is estimated to be only of USD$ 211,000, i.e., a reduction of 7%.

Discussion
The Dichato reconstruction plan incorporated hard mitigation measures to minimize future tsunami impacts in the bay, but our results show that when confronted to the 2010 tsunami, the level of direct damage reduction is only marginal. The most relevant mitigation aspect corresponds to a delay in arrival times of tsunami waves, as the inundation area of the second wave is reduced to approximately 40% in the presence of the seawall and channel. The time available to carry out an evacuation in the area located in the vicinity of the river mouth is thus increased to approximately an hour, as shown in Fig. 5d. This aspect is relevant for evacuation purposes, but it is a benefit that is not granted for other tsunamis since it is highly dependent on the characteristics of the seismic source; indeed for the 2010 tsunami, the largest tsunami waves that arrived to Dichato were not the first to come, but this situation could be different for other events. Indeed, the largest wave arrived 120 min after the earthquake in 2010 and inundated roughly the same area as it did in the case we simulated with the engineered structures in place. Hence, the conducted mitigation works might be efficient against more frequent and less destructive events, such as storm waves or surges or even minor tsunamis, but are not able to withstand the largest tsunami waves generated by an event similar to the one that occurred in 2010.
This finding motivates us to evaluate alternatives for the wall height and to produce a sensitivity analysis on this design parameter. The rationale behind this analysis is that one would expect that a higher wall height should produce larger direct benefits than the costs of construction, and eventually, and optimum constrained to specific objectives or criteria could be found. Thus, we consider different wall heights from + 0.5 m up to + 4 m from its actual configuration, to assess potential mitigation effects and damage reductions in view of determining an optimal wall height (Favier et al. 2014b;Shimozono and Sato 2016).
We thus perform additional tsunami simulations and re-assess the cost-benefit analysis considering walls 0.5, 1.0, 2.0, 3.0, and 4.0 m taller than the actual one of 5 m on average. Special structural designs or protection elements that would increase the construction costs of higher walls to ensure its correct performance in this situation are not included, since the costs are only assumed to be proportional to the height of the wall. It is also important to emphasize that a more comprehensive analysis would require the inclusion of a time window for life operation of the engineered works, the probability of exceedance of certain tsunami hazard level each year, and a cash flow for costs and benefits. Although difficult to quantify, these indirect impacts are also important and should be considered in the design and decision process.

3
The resulting inundation maps for the same 2010 event are shown in Fig. 6. As expected, as the wall height increases, the inundated area gradually decreases and the damage cost reduction increases, as summarized in Table 3. Indeed, when the height of the seawall is increased in 3 m, the flooded area is greatly reduced (43.3% of area reduction and an 88.9% reduction in damage costs). The different scenarios are summarized  Fig. 7, where it can be seen that the application of the methodology brings valuable technical information for decision making; nonetheless, this dimension constitutes only one element among many others that come into play before deciding investing in mitigation works. Indeed, the final decision has more to do with political issues, as it seems to have been the case in Dichato since we could not find any technical justification for the selected wall height in the official documents we could gather.
In addition to hard mitigation structures, which are not foolproof, other measures can also be considered. With this regard, sustaining education and awareness programs to ensure an effective evacuation response in case of a disaster should not be undermined (Favier et al. 2014b). Mitigation and risk communication strategies unbalanced toward hard measures may induce a false sense of security (Suppasri et al. 2013).
Urban planning is another important risk reduction tool that can be informed by the type hazard and damage assessment analysis presented here. The latter could be useful to improve the preparedness of Dichato against an event similar to the one that hit Chile in 2010. For example, if additional studies of potential evacuation are considered, restriction areas based on this information could be defined to ensure a prompt and safe evacuation (e.g., Imamura et al. 2012;León et al. 2019).

Conclusion
The present investigation provides an ex post analysis of the effectiveness of the mitigation measures that were built in the town of Dichato in the aftermath of the 2010 tsunami. This is the first study of this kind in Chile, but follows similar lines than research published after the 2011 Tohoku tsunami in Japan. Our results suggest that the seawall and the channeling of the river account for a direct damage cost reduction of nearly 7% if faced to a similar tsunami than the one that hit this town in 2010, with an estimated cost of implementation of USD $6.8 million. However, damage reduction could be increased if taller walls were considered, as a reduction of over 85% in damage costs is estimated with a wall 3 m taller than the reference one. Indeed, the direct benefits of the current measures are minor (approximately 200 thousand US dollars) and would deserve to be confronted with other indirect costs which should also include the impact on tourism and beach erosion over the operation lifetime of the mitigation works.
The most important effect of the conducted mitigation works is an increase in evacuation times for the area close to the river, which should bring additional indirect benefits that are not accounted for in the present study. The latter is explained by the particularities of the tsunami hydrodynamics at this location since the biggest wave is not much affected by the mitigation works, and it arrives 2 h after the nucleation of the earthquake. This can be a situation particular to this place and event, so these results should be treated as a first estimate of the benefit.
Despite this potential shortcoming, our simulations suggest that the mitigation works in place could be effective against minor tsunamis, storm waves or surges, but not against the biggest waves produced by a similar event than the one that occurred in 2010.
While the present research is important to provide an ex post assessment of the mitigation work effectiveness, a more comprehensive analysis for the definition of optimal measures should consider additional aspects that are not easy to evaluate such as probabilistic simulation of the tsunami hazard, risk perception, and cultural issues. Indeed, mitigation plans should be evaluated from more complex perspectives than direct damage reductions attributed to hard works, including community risk awareness, land use, and evacuation plans.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Appendix 1: Valuation of mitigation works
To perform an approximate valuation of the mitigation works, a simplified version of the wall and channels' structures was considered, as composed of fewer materials. This simplification was carried out because of practical matters. All the information used to perform this valuation was provided by the Dirección de Obras Portuarias (DOP) and the Ondac catalogue, and in Table 4 is a detail of the materials considered in both the wall and channel. The prices are according to the Chilean market and include workforce and other minor materials required.
The mitigation works can be better appreciated in Fig. 8.

Appendix 1.1: Wall
The wall's construction was divided in three sectors, each with a different set of profiles. The profiles each have different sizes and material proportions that define their construction costs, so to obtain the total cost of the wall's construction, it is necessary to calculate the amount of each material required for each profile along every set. Every profile's length together with the materiality information was obtained from a set of drawings provided by the DOP. The following tables specify the amount of each material used per profile and their length for each sector. The excavations along each sector were calculated using a mean area, which in the case of Litril Sector is 22.50 m 2 , for Etapa 1 Sector is 13.05 m 2 and for Estero Sector is 9.88 m 2 (Tables 5, 6, 7). It should be noted that the geometry of every profile was simplified into basic geometric elements, such as squares and triangles.

Appendix 1.2: Channel
In the case of the channel, there was one average profile used along the 286 m of its length. Thus, the material summary of it is shown in Table 8. 1 3

Appendix 2: Selection of fragility curves
The fragility curves used correspond to a work by Suppasri et al. (2012) where six levels of damage were defined as 1: minor damage, 2: moderate damage, 3: major damage, 4: complete damage, 5: collapse, and 6: washed away. All six levels were conceived for wood, masonry, reinforced concrete, and steel frames, as shown in Fig. 9.

Appendix 2.1: Fragility curve for mixed materials
A fifth typology is proposed in the present investigation, which corresponds to a mix between reinforced concrete and wood for a specific kind of two-story houses. To this end, we consider the fragility curves for both reinforced concrete and wood, where the first story's materiality is reinforced concrete, while the second is made from wood. Thus, the idea is to apply an inundation height to each curve according to the relative (with respect of the base of each story) level of inundation of the first and the second stories. Therefore, using Fig. 10 as a reference, an inundation level equal to h ′ would be considered for the RC's fragility curve, and in the case of the wood's curve, an inundation height equal to h * is used. The final percentage of damage becomes the average between both curves. For instance, if a house built with mixed materials is exposed to an inundation height of 2.5 m, the total damage would be an average of the damage on the RC first story with an inundation of h ′ = 2 m and the damage on the wooden second story with an inundation of h * = 0.5 m. The RC curve shows that with a 2 m inundation, the approximated damage probabilities are: 100% chance of minor damage, 95% chance of moderate damage, 80% chance of major damage, 40% chance of complete damage, 10% chance of collapse, and 0% of washed away. On the other hand, the wood curve shows that with a 0.5 m inundation, the approximated damage probabilities are 85%, 60%, 20%, 5%, 0%, 0%, respectively. With these numbers, it is possible to apply the methodology explained in Sect. 3.2.3, with Eq.
(2) and the information presented in Table 1.

Appendix 3: Valuation of buildings
All the information utilized for the valuation of the buildings' sample was obtained from the Servicio de Impuestos Internos (SII) and is presented in Table 9. The actual price of a building was considered to include the value of the terrain because of limitations in the available information.