The 2011 Tsunami event in the eastern coastal area of Japan caused a huge amount of damages or devastations on buildings. To this date, several field surveys have been conducted which provide detailed information about inundation areas and building damage characteristics in attacking east coastal areas by this tsunami. In this study, building damage data of Ishinomaki city, with special attention to the plain coast affected area, are classified and analyzed using data surveyed by the Ministry of Lands, Infrastructure and Transportation of Japan (MLIT) for more than 52,000 structures. The classification includes information on six levels of damage, four types of building materials and damages due to tsunami inundation for each building material which are necessary information for an effective hazard mitigation. Notably, damage level percentage distribution of different building materials is plotted for different inundation depth ranges in several sets of figures. This graphic illustration not only shows a better resistant performance of Reinforced Concrete (RC) and steel buildings over wood or other buildings for all inundation depth ranges, but also can explain clearly the inundation-induced damage behavior for each building material as well as the threshold depth for each damage level. Moreover, this research contains an analysis of vulnerable areas due to the coastal topography and the geographical factors. Surveyed data provided by Geospatial information authority of Japan (GSI) that classifies Ishinomaki plain coast area into three classes are compared with the damage map produced using an Analytical Hierarchy Process (AHP) methodology in ArcGIS 10.2 environment. The influence of key geographical features on tsunami-induced building damage, notably Kitakami river and water canals flooding, is taken into account with respect to the weighting of factors. A good agreement produced building damage map with surveyed GSI data shows the power of a GIS tool based on the AHP approach for tsunami damage assessment. The results of this study are useful to understand the damage behavior of buildings with different structural materials located in coastal areas vulnerable to the tsunami disaster.
A natural disaster event that occurs in a populated area is a hazardous event that causes a huge amount of destruction and fatalities. A tsunami is one of these natural hazards which can harm human and properties vastly since it happens in an inhabited area. The 2011 Great east Japan, tsunami that attacked the east coast of Japan causes more than 400,000 buildings damaged or destroyed [
Unfortunately, in different human societies that are in danger of such natural disasters, people are engaging in developing infrastructure and building constructions seemingly without vulnerability assessment of hazards impact. Therefore, it is essential to improve our understanding of building damage characteristics and coastal topography due to a tsunami event. Then, it will be possible to prepare an initial vulnerability and damage assessments of a disaster event as a tool to reduce its damage to the community for an effective hazard mitigation [
Early studies provide discussions about the influence of building constructions material and the number of stories on the building resistance to tsunami inundation in several case studies using surveyed building damage data [
It is shown by these studies that the RC or steel based constructions are stronger than wood or masonry based constructions for different inundation depth ranges through plotting of fragility curves in many case studies [
In addition to this field survey data classification, the effect of coastal topography, geographical features and buildings locations has been treated by means of an analysis of satellite imagery in the plain coast of Ishinomaki city by creating damage map. Previous studies have developed and analyzed integrative remote sensing techniques in assessing building vulnerability to tsunami hazard. High resolution satellite imagery, which has become available in the last few years, makes satellite remote sensing more useful in disaster management since even the damage status of individual buildings can be identified without visiting the sites of disasters [
In our case study, the satellite image of Ishinomaki plain coast, which is shown in
The experience of historical hazardous events through field surveys together with obtaining information
through GIS, management authorities pre-plans appropriate response activities in order to minimize the damages of a future similar event. These disaster risk managements include activities such as providing better evacuation routes and safe places or shelters, educations for public, construction of tsunami barriers and new regulations for building construction.
The first part of this study dedicated to evaluate the building materials behavior and damage because of tsunami inundation according to the field survey conducted by the MLIT provided in database format. In the second part the costal topography impact is treated through damage map produced through GIS incorporated with AHP methodology for the plain coast area. The coastal topography effect on inundation based building damage is investigated by defining shore line, elevation, slope, river and vegetation by means of GIS as a tool to analyze Ishinomaki plain coast satellite image. The analysis results aim to reflect potential differences in damage probabilities related to the tsunami inundation, construction materials and the buildings’ location according to the damage map.
According to MLIT survey the damage levels were categorized in six classes: 1) Low damage, (some part of building damage and water between ground and first floor), 2) Slightly low damage (half damage and water between ground and first floor), 3) Moderate damage (half damage), 4) High damage (water over than first floor and completely destroyed), 5) Collapse and 6) Washed away. Moreover, the main building characteristics, including building construction material, number of stories and buildings locations were included. An example of such building damage data for affected area in Ishinomaki city is represented in
In this table, building damages are classified based on the measure of tsunami inundation depth. Similar files were available for steel, wood and lightweight building materials named as others in our classifications.
Low | Slightly low | Moderate | High | Collapse | Washed away | Total | ||
---|---|---|---|---|---|---|---|---|
~0.5 m | 46 | 29 | 7 | 82 | ||||
0.5 ~ 1.0 | 39 | 160 | 32 | 1 | 1 | 233 | ||
1.0 ~ 1.5 | 8 | 94 | 38 | 4 | 2 | 146 | ||
1.5 ~ 2.0 | 1 | 15 | 23 | 10 | 7 | 2 | 58 | |
2.0 ~ 2.5 | 7 | 14 | 29 | 22 | 2 | 74 | ||
2.5 ~ 3.0 | 2 | 16 | 19 | 20 | 1 | 58 | ||
3.0 ~ 3.5 | 11 | 28 | 24 | 63 | ||||
3.5 ~ 4.0 | 12 | 15 | 22 | 4 | 53 | |||
4.0 ~ 4.5 | 6 | 13 | 13 | 1 | 33 | |||
4.5 ~ 5.0 | 10 | 4 | 14 | |||||
5.0 ~ 5.5 | 11 | 10 | 4 | 25 | ||||
5.5 ~ 6.0 | 7 | 12 | 4 | 1 | 24 | |||
6.0 ~ 6.5 | 1 | 1 | 9 | 11 | ||||
6.5 ~ 7.0 | 1 | 13 | 14 | |||||
7.0 ~ 7.5 | 1 | 1 | ||||||
7.5 ~ 8.0 | 0 | |||||||
8.0 ~ 8.5 | 2 | 1 | 1 | 4 | ||||
8.5 ~ 9.0 | 1 | 5 | 2 | 8 | ||||
9.0 ~ 9.5 | 2 | 6 | 8 | |||||
9.5 ~ 10.0 | 1 | 2 | 5 | 8 | ||||
10.0 ~ 10.5 | 2 | 2 | ||||||
10.5 ~ 11.0 | 4 | 3 | 7 | |||||
11.0 ~ 11.5 | 1 | 1 | ||||||
11.5 ~ 12.0 | 1 | 2 | 2 | 5 | ||||
12.0 ~ 12.5 | 0 | |||||||
12.5 ~ 13.0 | 0 | |||||||
13.0 ~ 13.5 | 1 | 1 | ||||||
13.5 ~ 14.0 | 0 | |||||||
14.0 ~ 14.5 | 0 | |||||||
14.5 ~ 15.0 | 0 | |||||||
15.0 ~ 15.5 | 0 | |||||||
15.5 ~ 16.0 | 0 | |||||||
16.0 ~ 16.5 | 0 | |||||||
16.5 ~ 17.0 | 0 | |||||||
17.0 ~ 17.5 | 0 | |||||||
17.5 ~ 18.0 | 0 | |||||||
18.0 ~ 18.5 | 0 | |||||||
18.5 ~ 19.0 | 0 | |||||||
19.0 ~ 19.5 | 0 | |||||||
19.5 ~ 20.0 | 0 | |||||||
20.0 m~ | 0 | |||||||
Total | 933 | |||||||
The data were provided for 63,157 buildings where the number of buildings in the inundated area of Ishinomaki city and ones in the plain coast area (look to
As illustrated in the
Other important buildings characteristic provided by MLIT was number of stories. It has been explained in some previous studies that more number of stories per building the stronger the building construction because of their structural design to support heavier loads. Since based on MLIT database about 98% of buildings found to be two floors with few multi stories buildings that are located mostly in safe areas, building damage due to the number of stories is neglected in this study.
The MLIT survey analysis method in this study covers an investigation on damage level distribution based on
Building Materials | Ishinomaki City | Ishinomaki-Plain Coast |
---|---|---|
Concrete | 933 | 692 |
Steel | 2450 | 1526 |
Wooden | 45,976 | 24,997 |
Others | 13,798 | 6258 |
Total | 63,157 | 33,473 |
inundation depth, building damage distribution with considering both structural material and inundation depth in Ishinomaki city. For necessary occasions, the similar data analysis is applied to buildings located in the plain coastal area for comparison. Several data classifications were performed to extract valuable charts to assess the impact of construction material in amount of damage. Finally, an analysis of the impact of the coastal topography in the amount of building damage for 33,473 buildings of Ishinomaki city plain coast is provided. For this purpose, Arc GIS software tool based on AHP method is used to analyze the satellite image and classify the vulnerable areas.
In this study, remote sensing data captured by GeoEye-1 a Very High Resolution (VHR) optical satellite with a pixel size about 0.41 m was used to identify critical geographical elements such as buildings, transport infrastructure and the inundation area due to the tsunami event in order to estimate the damage and risk in plain coast area of Ishinomaki City in Miyagi Prefecture as shown in
We extracted the factors of elevation and slope from a Digital Elevation Model (DEM) obtained from the Geospatial Information Authority of Japan while the Normalized Difference Vegetation Index (NDVI) for vegetation density is extracted from Geoeye-1 image and Image Analysis toolbar in ArcGIS 10.2.1 software. The shoreline distance and river distance were measured from the vector maps of the study area.
The damage by tsunami is then estimated by applying AHP to all the mentioned factors. The damage assessment can be displayed via GIS in terms of spatial multi - criteria analysis on a map of the tsunami damage area.
A digital elevation model was created from elevation data obtained from GSI. DEM was downloaded from http://fgd.gsi.go.jp/download/GsiDLSelFileServlet. The mesh elevation data created by interpolating to the elevation point at center point in 0.2 seconds (about 5 meter) mesh that is from the value of elevation (ground data) measured by airborne laser. The height accuracy of the 5 m mesh elevation from the airborne laser is less than 0.3 m and the standard deviation of the altitude acquisition position (from the photogrammetry) is less than 1.0 meter. The data were in JPGIS (format) converted to shapefile in point format using base map viewer converter software version 4.00 (FGDV) provided by GSI.
Finally, this point format was converted to raster for creating the digital elevation model via Arc GIS 10.2.1 software and raster interpolation toolset function.
The slope was determined as the rate of maximum change in z value from each cell of the satellite image. The use of a z-factor is essential for correct slope calculations when the surface z units are expressed in units different from the ground x, y units. The range of values in the output depends on the type of measurement units. The range of slope values is 0 to 90 for degrees and 0 to essentially infinity for percent rise. We created a slope map using the surface creation and analysis tools of the ArcGIS 10.2.1 software to use a third-order finite-difference method for calculating the slope [
The distance from the shoreline was created in a polyline file for buffering the distance from the shoreline to the land. We computed the distance using the proximity and the Euclidean distance analyst tool in the ArcGIS 10.2.1 software. The distance is based on the historical report of the maximum run-up in the area of study. We used the “Equation (1)” to classify coastal proximity and shoreline distance buffering from Bretschneider and Wybro (1976) [
where Xmax is the maximum reach of the tsunami over land, and Y0 is the tsunami height at the shore line.
We classified distance buffers in five classes based on the “Equation (1)” in order to create a tsunami damage map. It explains that 4.55 m to 7.09 m of run-up can reach a distance of 489.94 m from the shoreline, 7.09 m to 9.64 m of run-up can reach 885.76 m, 9.64 m to 12.18 m of run-up can reach 1332.84, 12.18 m to 14.73 m of run-up can reach 1821.46 m and 14.73 m to 17.27 m of run-up can reach more than 2345.53 m.
The Kitakami River, the fourth largest river in Japan and the downstream branches into two water canal, Old- Kitakami river and Kitakami river. Old-Kitakami river flows on a fertile plain which has been highly developed for agriculture and industries and this water canal passes through the study area, while Kitakami river flows through a narrow valley into a small bay faced to the Pacific Ocean. It can be seen that the tsunami had reached 8 to 10 km or more off the coastline towards inland around Kitakami river and Old Kitakami river and the elevation in this region was extremely low, ranging between 0 and 2 m [
Kitakami-unga and Kitakami canals are artificial canals with an average elevation of −9 m below sea level along 15.8 km. Kitakami canal that ran across the city and then inundated the inland area [
The elevation, river distance and shoreline distance were classified into five classes of damage using the Jenks natural breaks method. This classification method indicates by picking the class breaks in best group similar values, maximize the differences between classes and minimize value differences between data within the same class and emphasize the differences between the reclassified classes.
The NDVI is a common tool for identifying and characterizing vegetation and a measure of the difference in reflectance between these wavelength ranges that take values between −1 and 1, while vegetated areas produce values starting around 0.4 and approaching 1.0 and values < 0 indicating no vegetation [
“Equation (2)” was used for NDVI calculation and specification:
In which for GeoEye-1 Band 4 is NIR (near infrared band) and Band 3 is red.
We classified elevation and the distance for river, canal and shoreline based on the values described in
Cell-based modeling in spatial analysis was used to specify the damage area due to tsunami hazard. Cells are classified into five classes of damage in the numbers of 1, 2, 3, 4 and 5, which represent low, slightly low, medium, slightly high, and high damage classes for different inundation depth.
AHP is a Multi-Criteria Decision-making Analysis (MCDA) approach introduced by Saaty (1977, 1980). It relies on the series of pair-wise comparisons between each factor relative to other factors to make a scaled set of preferences in ranges from 1 to 9 to score the importance of each factor (Saaty, 2008).
All of tsunami damage factors are overlaid and weighted based on their dominant influences in determining the class of tsunami induced damage. The relative importance of each factor within the hierarchy is determined
Damage class | Factors name | |||||||
---|---|---|---|---|---|---|---|---|
Elevation (m) | Slope (%) | Shoreline distance (m) | River distance (m) | Vegetation | ||||
Kitakami | Unga canal | Index | Density | |||||
High (5) | −4 - −1 | 0 - 1 | 0 - 489.94 | - | - | −0.998 - 0.151 | Low | |
Slightly high (4) | −1 - 0 | 1 - 3 | 489.94 - 885.76 | - | - | 0.151 - 0.164 | Slightly low | |
Medium (3) | 0 - 2 | 3 - 6 | 885.76 - 1332.84 | 0 - 160 | - | 0.164 - 0.188 | Medium | |
Slightly low (2) | 2 - 8 | 6 - 8 | 1332.84 - 1821.46 | 160 - 680 | 0 - 300 | 0.188 - 0.218 | Slightly high | |
Low (1) | >8 | >8 | 1821.46 - 2345.53 | >680 | >300 | 0.218 - 0.556 | High |
Score | Definition | Explanation |
---|---|---|
1 | Equal importance | Two factors contribute equally to the objective. |
3 | Weak importance of one over another | The judgment is slightly favor one factor over another. |
5 | Essential or strong importance | The judgment is strongly favor one factor over another. |
7 | Demonstrated importance | A factor is strongly favored and its dominance is demonstrated in practice. |
9 | Absolute importance | The evidence favoring one factor over another is of the highest possible order of affirmation. |
2, 4, 6, 8 | Intermediate values between the two adjacent judgments | When compromise is needed. |
by their weights (Saaty, 1977 & 1980) with a pairwise comparison as shown in
The first eigenvector is computed based on the pairwise comparison matrix shown in
AHP is subjective and tolerates inconsistency through the amount of redundancy by providing a measure of inconsistency assessment, which is shown by Consistency Ratio (CR). If the value of CR is smaller or equal to 10%, the inconsistency is acceptable. CR indicates the probability that the matrix judgments were randomly generated and it is defined as the ratio of the Consistency Index (CI), which is the degree of logical consistency among pair-wise comparisons, to the Random consistency Index (RI) which is the average CI value of randomly- generated comparison matrices [
“Equations (3)” and “Equation (4)” describe the algorithms for CR and CI calculation:
In which:
λmax is the maximum eigenvalue of the judgement matrix and calculated from the sum of all factors and is multiplied by its eigen vector which is 42.16%.
n is the size of the comparison matrix. In this study, n = 5.
The RI is based on the random consistency index as shown in
So the consistency index, CI, is calculated 0.093 and CR is 8.37%.
We calculated each raster cell-by-cell basis of the factor to its weight using any GIS system with overlay ca-
Pairwise comparison | Factors name | ||||
---|---|---|---|---|---|
Elevation | Slope | Shoreline distance | River distance | Vegetation | |
Elevation | 1.00 | 2.00 | 0.33 | 0.50 | 2.00 |
Slope | 0.50 | 1.00 | 0.25 | 0.33 | 1.00 |
Shoreline distance | 3.00 | 4.00 | 1.00 | 2.00 | 5.00 |
River distance | 2.00 | 3.00 | 0.50 | 1.00 | 4.00 |
Vegetation | 0.50 | 1.00 | 0.20 | 0.25 | 1.00 |
Random Consistency Index | Matrix Size | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
RI | 0.00 | 0.00 | 0.52 | 0.89 | 1.11 | 1.25 | 1.35 | 1.40 | 1.45 | 1.49 |
pabilities. A Weighted Linear Combination (WLC) analysis is very straightforward in a raster GIS, and each factors are combined by applying a weight, to each followed by a summation of the results to proceed a suitability map. “Equations (5)” describes the suitability calculation [
where:
S = Suitability.
Wi = Weight of factor.
Xi = the criterion score of factor i.
We estimated raster calculator in map algebra menu using the spatial analyst tools of ArcGIS 10.2.1 to produce damaged map by applying “Equations (5)”.
The first part is representing results of building damage for different construction materials. Second part shows the result of GIS for building damage map.
In this section, based on surveyed data several statistical analysis results of building damage for different construction materials are provided and discussed. The damage level percentage distribution of different building materials are plotted for different inundation depth ranges in several sets of figures: 1) Tsunami inundation depth only 2) Construction material only, 3) Construction material and inundation depths together.
The tsunami inundation depth for buildings located in Ishinomaki city was obtained from the MLIT surveyed data and also Geoeye-1 satellite image. According to this data, we classified and summarized the distribution of damage levels based on tsunami inundation depth on an interval of 0.5m for the whole city, regardless of building material type. The result is graphically represented in
According to
The wave forces brought by a tsunami cause great damages to fragile buildings. Standing buildings in the path of the waves tend to block the water, and thus the pressure on them increases. The water force can be calculated by using the wave pressure distribution by Tanimoto formula. For example, water just 2m deep will impose pressure of approximately 20 KN per square meter, much more than any normal construction can remain undamaged. The deeper the water is, the greater the pressure is [
This part is dedicated to comparing statistically different building material damage during the east Japan great
tsunami 2011 for Ishinomaki city with 63,157 buildings for the whole city and with 33,473 buildings for Ishinomaki plain coast as a case study. Each damaged building was measured and classified according to its damage level (six levels were defined by MLIT), construction material (RC, steel, wooden and other light weight construction such as clay or brick mentioned in text as others).
The damage data are classified according to the building material as shown in
The same extracted date is plotted as a percentage distribution of different building materials and shown in
As it is mentioned earlier, the deeper the water, the greater the pressure and the greater the amount of damage is expected. By considering the inundation depth on an interval of 0.5 - 8.5 m, the damage data for the whole Ishinomaki city is classified according to the type of building materials. For a given inundation depth interval, the percentage contribution of each damage level is calculated. The graphical results are shown graphically in
in the inundation depth range less than 2 m only low and slightly low damages are expected. For the inundation depth of 2.5 - 6.0 m RC buildings are mostly high damaged since buildings located in inundation depth over 6.0 m are all collapsed. For steel buildings as depicted in
Wooden buildings shown in
With considering inundation depth, the same data used for
It is also comprehending that for the
Damage Index | Damage Area Classification | |||
---|---|---|---|---|
Damage Classification | Damage Value in Tsunami Damage Map | Area (km2) | Area (%) | |
1 | Low | 1.30 - 1.59 | 2.20 | 18.69 |
2 | Slightly Low | 1.59 - 1.87 | 3.08 | 22.76 |
3 | Medium | 1.87 - 2.34 | 3.02 | 21.35 |
4 | Slightly High | 2.34 - 2.97 | 3.22 | 21.71 |
5 | High | 2.97 - 4.32 | 2.65 | 15.48 |
corded by “GSI” and “2011 Earthquake Tsunami Joint Survey Group” shown in
Based on the result of our study, which is shown in
Building damage data of Ishinomaki city, with special attention to the plain coast affected area, were classified and analyzed using data surveyed by the MLIT for more than 52,000 structures. The classification includes information on six levels of damage, four types of building materials and damage percentage due to tsunami inundation for each building material which are necessary information for an effective hazard mitigation.
The main findings of this study and their applicability are summarized as below:
・ Building damage data of Ishinomaki city, with special attention to the plain coast affected area, were classified and analyzed.
・ Damage level percentage distribution for different building materials and buildings percentage distribution for different damage levels for Ishinomaki city and its plain coast were illustrated.
・ For each building material, damage level percentage distribution of different building materials was plotted under inundation depth range of 0 - 8.5 m. These graphs show the material behavior due to tsunami inundation to assess the damages which are comparable to fragility curves. For each damage level, different inundation threshold depths were also mentioned.
・ RC, steel, wooden and other fragile building materials behavior during a tsunami event was compared under different inundation depth. The results show a better resistant performance of RC and steel buildings over wood or other buildings.
・ Calculation of damage map and its comparison with a damage map based on surveys, allow us to understand the importance of each geographical factor that has significant influence on damage. By means of GIS tool, the damage areas are followed by building material behavior schematics. One can assess the amount of damage for any damage area.
These are important considerations in designing a tsunami evacuation building, material and location, and assessing the damages for other cities in danger of future tsunami attach.
In conclusion, by using multiple geospatial factors of topographic elevation, relation to tsunami direction, coastal proximity, and coastal shape incorporated by the AHP, an appropriate pair-wise comparison of AHP is proposed to construct a weighting scheme for the geospatial factors and assessing tsunami damage. However, the damage is not only controlled by the tsunami height which is dependent highly to coastal topography, but also highly dependent on the building materials. By knowing the material behavior of different materials for each inundation depth, as shown graphically in this study, one can assess with a good precision the amount of damage. We proposed using the illustrated figures of damage for different building materials due to different inundation depths in this paper, together with GIS which produces a damage map as a powerful approach to assess the damages on buildings for an area attacked by tsunami storms.
Our research can be employed to evaluate building damage located in vulnerable areas due to tsunami hazard in future natural disasters.
The technical and financial supports by Research Center for Urban Safety and Security (RCUSS) of Kobe University and Shahid Beheshti University are gratefully acknowledged. Authors are thankful to the GSI for providing DEM data and vector maps of the study area and 2011 Tohoku Earthquake Tsunami Joint Survey for the survey data in the study area.
Mohammad Reza Poursaber,Yasuo Ariki, (2016) Integrated GIS, Remote Sensing and Survey Data for Damage Assessment of Buildings in Tsunami Event, Ishinomaki City, Japan. Journal of Geographic Information System,08,260-281. doi: 10.4236/jgis.2016.82023