The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Chapter 4 Seismic Risk Assessment
Seismic risk assessment was carried out as the update of the risk assessment of 2002 JICA project of “The Study on Earthquake Disaster Mitigation in the Kathmandu Valley” for the purpose of providing basic information for the formulation of disaster risk reduction and management plans for pilot municipalities. The assessment was performed based on the latest available data and taking into account the new research results on risk assessment as well as the characteristics of ground motion, building damage and human casualties caused by the Gorkha Earthquake, which occurred on 25 April 2015, right before the commencement of the project. The main features of risk assessment of this project are summarized as below.
No building inventory, which is important for risk assessment, exists in KV. In order to establish building inventory, building survey for all buildings was made by the project for former Lalitpur Sub-metropolitan City and Bhaktapur municipality. In addition, building information for all buildings in Budhanilkantha and former Karyabinayak municipalities was collected. The building inventory for the other municipalities were created by estimation based on satellite image and sample structure type survey.
The damage function of buildings was developed based on the building damage data of the Gorkha Earthquake as well as the seismic resistant capacity analysis of typical structures of Nepal.
Nonlinear flexural deformation of RC piers subject to seismic excitation was estimated by a statistical method and was used to classify the damage degree of bridges.
Human casualties were estimated based on the death rate and injured rate of the Gorkha Earthquake for different damage levels and different structure types. Three earthquake occurrence scenes: night, weekday noon and weekend afternoon, were considered for human casualty estimation.
Building damage and human casualties were also estimated for 2030 for the purpose of determination of disaster risk reduction target, assuming different cases of structure type distribution.
The contents of seismic risk assessment cover structural damage of buildings, roads, bridges,
4-1 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report) water supply pipelines, sewage pipelines, power poles and mobile base transceiver station (BTS) towers as well as human casualties and the direct economic loss and impact on the tourism industry. The summary of risk assessment results is shown in Table 4.1
Seismic risk assessment results from one scenario ground motion are shown in the main report as an example. A map book, which includes all of the risk assessment results together with the basic information for seismic hazard and risk assessment as well as seismic hazard results, was compiled for easy reference. The contents of the map book are given in Table 4.2.
REMARKS
The scenario earthquakes are not the prediction of future earthquakes. Risk assessment was carried out based on scientific research and investigation results but with inevitable assumptions. Its results might have uncertainties and are not the guarantee of the future damage of a scenario earthquake. The purpose of risk assessment is to provide basic information for the development of disaster risk reduction and management plan of pilot municipalities in KV.
4-2 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.1 Summary of risk assessment results
Physical Damage Economic Loss (mil. NPR)*1 Human Casualty (Population: 2016: 2,786,929; 2030: 3,805,926) Category Scenario Earthquake Ground Motion Scenario Earthquake Ground Motion Scenario Earthquake Ground Motion WN CNS-1 CNS-2 CNS-3 WN CNS-1 CNS-2 CNS-3 WN CNS-1 CNS-2 CNS-3 Night (Weekday and weekend) 3,034 9,133 22,179 35,726 24,961 65,314 136,060 199,643 Death 0.11% 0.33% 0.80% 1.28% Heavy damage 11,880 35,766 86,861 139,914 (EMS DL4&5) Injured 0.43% 1.28% 3.12% 5.02% 5.6% 14.7% 30.6% 44.9% 279,031 642,743 1,196,080 1,613,314 Evacuee 10.01% 23.06% 42.92% 57.89% Weekday (noon, 12:00) 2,784 8,282 19,959 31,956 21,967 42,940 62,691 67,418 Death 0.10% 0.30% 0.72% 1.15% Building (2016) Moderate damage 132,999 371,003 761,531 1,098,353 10,905 32,435 78,168 125,152 (Total building 444,554) (EMS DL3) Injured 0.39% 1.16% 2.80% 4.49% 4.9% 9.7% 14.1% 15.2% 285,850 652,798 1,206,530 1,619,792 Evacuee 10.26% 23.42% 43.29% 58.12% Weekend (afternoon, 18:00) 2,123 6,393 15,526 25,008 43,564 67,770 77,713 70,462 Death 0.08% 0.23% 0.56% 0.90% Slight damage 8,316 25,036 60,803 97,940 (EMS DL2) Injured 0.30% 0.90% 2.18% 3.51% 9.8% 15.2% 17.5% 15.9% 279,942 645,483 1,202,734 1,624,032 Evacuee 10.04% 23.16% 43.16% 58.27% 33,763 88,681 185,796 273,269 4,121 12,508 30,583 49,381 Case-0 Death 5.6% 14.6% 30.6% 45.1% 0.11% 0.33% 0.80% 1.30% 28,377 79,075 171,977 258,044 3,434 11,017 27,930 46,017 Case-1 Death 4.7% 13.0% 28.4% 42.5% 16.7% 11.9% 8.7% 6.8% 13,627 56,452 146,361 234,477 1,721 8,135 24,356 42,526 Case-2 Death Building (2030) 2.2% 9.3% 24.1% 38.7% 58.2% 35.0% 20.4% 13.9% (Heavy damage, Total building 606,506)*2 12,162 49,970 131,095 213,481 1,438 6,733 20,526 36,715 Case-3 Death 2.0% 8.2% 21.6% 35.2% 65.1% 46.2% 32.9% 25.6% 16,147 52,413 129,904 210,181 2,052 7,887 23,086 41,146 Case-4 Death 2.7% 8.6% 21.4% 34.7% 50.2% 36.9% 24.5% 16.7% 11,138 41,230 111,854 189,357 1,476 6,524 20,842 38,733 Case-5 Death 1.8% 6.8% 18.4% 31.2% 64.2% 47.8% 31.9% 21.6% 237 737 1,654 2,486 Heavy damage 444 1,545 4,002 6,555 4.1% 12.9% 28.9% 43.4% Death
School 253 539 810 875 0.05% 0.18% 0.47% 0.77% Moderate damage 20,462 51,231 98,171 134,932 (Total building 5,731) 4.4% 9.4% 14.1% 15.3% 1,739 6,051 15,673 25,671 568 916 1,057 960 Injured Slight damage 9.9% 16.0% 18.4% 16.8% 0.20% 0.71% 1.84% 3.02% 20 64 153 235 Heavy damage 3.4% 11.0% 26.2% 40.2%
Health facility 24 55 83 94 Moderate damage 27,534 68,588 165,683 232,782 (Total building 584) 4.1% 9.4% 14.2% 16.1% 51 85 105 97 Slight damage 8.7% 14.6% 18.0% 16.6% 20 59 126 186 Heavy damage 4.2% 12.3% 26.4% 38.9%
Government building 20 44 66 73 Moderate damage 2,444 8,669 16,514 22,708 (Total building 478) 4.2% 9.2% 13.8% 15.3% 44 71 85 80 Slight damage 9.2% 14.9% 17.8% 16.7%
Length in landslide 0 6.6 98.5 390.6 area (km) Road*3 0.0% 0.1% 1.7% 6.7% 0 471 1,620 2,878 (Total length 5,811 km) Length in liquefaction 0 76.1 274.9 455.3 area (km) 0.0% 1.3% 4.7% 7.8% 0 1 12 32 Heavy damage 0.0% 2.2% 26.7% 71.1%
Bridge 2212711 Moderate damage 377 898 1,359 1,914 (45 bridges assessed)*4 4.4% 46.7% 60.0% 24.4% 18 17 6 2 Slight damage 40.0% 37.8% 13.3% 4.4%
982 1,921 3,496 5,161 Water supply (Existing) Damage points 36 71 129 191 (Total length 1,167 km) 0.84 1.65 3.00 4.42
124 255 460 676 Water supply (Planned) Damage points 5 9 17 25 (Total length 699 km) 0.18 0.36 0.66 0.97
Sewage 4.81 8.15 11.94 18.21 Damage Length (km) 76 135 200 290 (Total length 1,192 km) 0.4% 0.7% 1.0% 1.5%
1,327 3,991 9,156 13,992 Power distribution Pole broken 19 56 129 197 (Total pole 190,851) 0.7% 2.1% 4.8% 7.3%
Mobile BTS tower 43 143 372 601 Tower damage 82 272 707 1,142 (Total tower 1,043) 4.1% 13.7% 35.7% 57.6% Source: JICA Project Team
4-3 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.2 Contents of Map Book
No. Caption No. Caption A. Basic Conditions C. Earthquake Risk Assessment A-1 Study Area and Pilot Municipality C-1 Distribution of Heavily Damaged Building & Ratio in 2016 (WN) A-2 Distribution of Population in 2016 at Night C-2 Distribution of Heavily Damaged Building & Ratio in 2016 (CNS-1) A-3 Distribution of Estimated Population in 2030 at C-3 Distribution of Heavily Damaged Building & Ratio in 2016 Night (CNS-2) A-4 Building Distribution in 2016 C-4 Distribution of Heavily Damaged Building & Ratio in 2016 (CNS-3) A-5 Estimated Building Distribution in 2030 (without C-5 Distribution of Moderately Damaged Building & Ratio in 2016 BSPS) (WN) A-6 Estimated Building Distribution in 2030 (with BSPS C-6 Distribution of Moderately Damaged Building & Ratio in 2016 Case-1) (CNS-1) A-7 Estimated Building Distribution in 2030 (with BSPS C-7 Distribution of Moderately Damaged Building & Ratio in 2016 Case-2) (CNS-2) A-8 Estimated Building Distribution in 2030 (with BSPS C-8 Distribution of Moderately Damaged Building & Ratio in 2016 Case-3) (CNS-3) A-9 Estimated Building Distribution in 2030 (with BSPS C-9 Distribution of Heavily Damaged Building for 2030 without Case-4) BSPS A-10 Estimated Building Distribution in 2030 (with BSPS C-10 Distribution of Heavily Damaged Building for 2030 with BSPS Case-5) Case-1 A-11 Distribution of School Building C-11 Distribution of Heavily Damaged Building for 2030 with BSPS Case-2 A-12 Distribution of Health Facility Building C-12 Distribution of Heavily Damaged Building for 2030 with BSPS Case-3 A-13 Distribution of Government Building C-13 Distribution of Heavily Damaged Building for 2030 with BSPS Case-4 A-14 Road Network C-14 Distribution of Heavily Damaged Building for 2030 with BSPS Case-5 A-15 Emergency Transportation Road Network (ETRN) C-15 School Building Damage Proposed by JICA RRNE A-16 Distribution of Bridge C-16 Health Facility Building Damage A-17 Distribution of Water Supply Network (Existing) C-17 Government Building Damage A-18 Distribution of Water Supply Network (Planned) C-18 Possible Damage of Road by Liquefaction A-19 Distribution of Sewage Network C-19 Possible Damage of Road by Slope Failure A-20 Estimated Power Pole Distribution C-20 Possible Link Blockage of Road by Building Damage A-21 Distribution of Mobile BTS Tower C-21 Possible Link Blockage of ETRN by Building Damage C-22 Damage of Bridge B. Earthquake Hazard Assessment C-23 Priority Rank of Bridge for Seismic Strengthening B-1 Geomorphological Map C-24 Distribution of Water Supply Network Damage (Existing) B-2 Altitude Distribution Map C-25 Distribution of Water Supply Network Damage (Planned) B-3 Distribution of Collected Borehole Data C-26 Distribution of Sewage Network Damage B-4 Rock Depth Distribution and Location of Borehole C-27 Distribution of Power Pole Damage B-5 Location of Microtremor Measurement C-28 Distribution of Mobile BTS Tower Damage B-6 Geological Cross-Section EW C-29 Distribution of Death in 2016 at Night B-7 Geological Cross-Section NS C-30 Distribution of Death in 2016 at Weekday Noon B-8 Estimated AVS30 from Ground Model C-31 Distribution of Death in 2016 at Weekend Afternoon B-9 Predominant Period of Ground C-32 Distribution of Death Ratio in 2016 at Night B-10 Fault Model of Scenario Earthquake C-33 Distribution of Death Ratio in 2016 at Weekday Noon B-11 Peak Ground Acceleration Distribution C-34 Distribution of Death Ratio in 2016 at Weekend Afternoon B-12 Peak Ground Velocity Distribution C-35 Comparative Vulnerability MAP B-13 Seismic Intensity (MMI) Distribution B-14 Distribution of Liquefaction in Rainy Season B-15 Distribution of Liquefaction in Dry Season B-16 Distribution of Slope Failure B-17 AVS30 Map base on Geomorphological Unit B-18 Liquefaction Susceptibility Map B-19 Earthquake Induced Slope Failure Susceptibility Map Source: JICA Project Team
4-4 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
4.1 Necessity and objectives of Seismic Risk Assessment
According to EM-DAT (Figure 4.1.1), natural disasters in the world has increased sharply from around the 1970s, peaking at around 2005, and remaining at a higher level from then. Natural disasters result in not only life and economic loss, but also as obstructions to the sustainable development of the economy. In order to reduce disaster risk globally, the UN had initiated the International Decade for Natural Disaster Reduction (IDNDR) for the 1990s in 1987 at its 42nd General Assembly. After IDNDR, Hyogo Framework for Action (HFA) was adopted in the 2nd World Conference on Disaster Risk Reduction in 2005. Based on HFA, the world has promoted disaster risk reduction activities and produced steady results with the development of legal frameworks, establishment of DRR organizations, strengthening of collaboration among organizations, capacity development, and structural and non-structural measures to reduce vulnerability. On the other hand, the activities which lead to direct reduction of risk are comparatively limited due to the lack of budget and knowledge. In this circumstance, the Sendai Framework for Disaster Risk Reduction (Sendai Framework) was adopted in the 3rd World Conference on Disaster Risk Reduction in 2015, where sustainable development with mainstreaming DRRM is emphasized through advocating for four priorities for action and assigning seven global targets. Risk assessment is closely related to the Priority for Action 1 of Sendai Framework: Understanding Disaster Risk and Priority for Action 2: Strengthening Disaster Risk Governance to Manage Disaster Risk. Risk assessment is essential to provide basic information for the development of disaster risk reduction policy and plan to achieve the disaster risk reduction goal for human casualties, economy loss and critical infrastructure damage.
Seismic risk assessment of this project is planned to be utilized for the determination of risk reduction goal, the creation of disaster risk reduction and management plan and support for CBDRRM activities of pilot municipalities. It has been used for the development of the Kathmandu Valley Resilient Plan (KVRP) under JICA RRNE project and can be expected to be used for preparation of the disaster risk reduction and management plans of the municipalities in KV other than the three pilot municipalities. It also provides useful information for the formulation of Business Continuity Plans (BCP) for local government as well as utility companies.
4-5 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Drought
Storm
Flood
Earthquake
Source: EM-DAT Figure 4.1.1 Natural disaster trends in the world
A common issue for seismic risk assessment in many countries including Nepal is the lack of basic information such as earthquake damage data and structural inventory. Insufficient earthquake damage data makes it difficult to prepare a specific fragility curve. The lack of inventory will directly affect the precision of risk assessment. This project faced the same issues and some assumptions and presumptions were made to compliment the insufficiency. During the process of risk assessment, the damage data and findings from the Gorkha Earthquake have been used whenever possible. The update of risk assessment may be required in the future in case new findings from the Gorkha Earthquake are obtained and/or the building, infrastructure and lifeline inventories are completed.
4.2 Scenario Earthquakes and Their Occurrence Scene
Three scenario earthquakes were adopted in seismic hazard assessment. Four possible ground motion levels were applied for Central Nepal South Scenario Earthquake because it is located adjacent to the epicenter of the Gorkha Earthquake and the ground motion attenuation characteristics of the Gorkha Earthquake have to be considered. Considering the purpose of the project, among the six ground motion levels of seismic hazard results, four of them are chosen as the target ground motion for risk assessment. Besides, different scenario earthquake occurrence scenes were considered for risk assessment to have alternatives of disaster because the time when an earthquake will happen is not predictable.
4.2.1 Scenario Earthquakes and Ground Motion Used for Risk Assessment
Three scenario earthquakes, i.e. Far-Mid Western Nepal Scenario Earthquake (M = 8.6), Western Nepal Scenario Earthquake (M = 7.8) and Central Nepal South Scenario Earthquake (M = 7.8) were considered in hazard assessment. Ground motion for Far-Mid Western Nepal Scenario Earthquake and Western Nepal Scenario Earthquake was estimated directly from
4-6 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
attenuation formula, while that for Central Nepal South Scenario Earthquake was estimated from the attenuation formula with four modification factors of 1/3, 1/2, 2/3 and 1/1 to consider smaller accelerations recorded in KV from Gorkha earthquake. Considering the possibility of ground motion level and reality of disaster risk reduction countermeasures, the ground motion to be used for risk assessment was discussed in WG1 and 4th JCC meeting, which lead to the conclusion of that (1) Far-Mid Western Nepal Scenario Earthquake would not be forwarded for risk assessment because it is very far from KV and will not cause significant damage due to its weak ground motion, (2) ground motion from Western Nepal Scenario Earthquake (without modification) would be used for risk assessment and, (3) ground motions with modification factor of 1/3, 1/2 and 2/3 from Central Nepal South Scenario Earthquake would be used for risk assessment, and ground motion with modification factor 1/1 (no modification) would not be used for risk assessment since it may largely over-estimate ground motion. As a result, four ground motions noted as WN, CNS-1, CNS-2 and CNS-3 were finally selected for risk assessment as shown in Table 4.2.1. To identify different ground motions from one scenario earthquake, the ground motion used for risk assessment are called scenario ground motion hereinafter.
Table 4.2.1 Scenario ground motion for risk assessment Scenario Earthquake Modification Factor for PGA Notation Western Nepal Scenario Earthquake 1/1 WN 1/3 CNS-1 Central Nepal South Scenario Earthquake 1/2 CNS-2 2/3 CNS-3 Source: JICA Project Team
4.2.2 Earthquake Occurrence Scenes
A difference in occurrence time of earthquakes will not cause a difference in the structural damage of building, infrastructure and lifeline facilities, but do affect human casualties due to the different populations inside building at the time of an earthquake occurrence. In this regard, the earthquake occurrence scenes (occurrence time) for 2016 was considered for night, weekday noon and weekend afternoon based on the survey results of inside building ratio in different day (weekday and weekend) and time. In this project, for the purpose of determination of risk reduction targets, risk assessment was also conducted for 2030, when two scenes were considered: (1) extrapolation of building structure type of 2016, meaning without any strengthening measures to upgrade building seismic performance and (2) assuming total building seismic capacity is somehow improved with five supposed cases to consider possible structure type distribution in the future.
4-7 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
(1) Earthquake Occurrence Scenes for 2016
The different situations, mainly of population distribution, weekday, weekend (holiday), daytime and night were considered for the determination of earthquake occurrence scenes of 2016. As a result, three scenes at night (2:00am), weekday noon (12:00pm) and weekend afternoon (18:00pm) were decided as shown in Figure 4.2.1. Population inside building for each occurrence scene was assigned based on the questionnaire survey. The spreading of fire was not taken as a factor for earthquake occurrence scene because of its low possibility in Nepal. Population distribution for night comes from a census and that for weekday daytime was derived from the origin-destination (OD) survey data, obtained from The Project on Urban Transport Improvement for Kathmandu Valley, JICA, which accounts for the movement of workers and students from their home to office or school. The main features of different occurrence scenes are summarized in Table 4.2.2.
Source: JICA Project Team Figure 4.2.1 Earthquake occurrence scenes for 2016
4-8 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.2.2 Features of different earthquake occurrence scenes Scene Feature More human casualties may occur as compared to the daytime Night (2:00 am) May cause delay on search and rescue due to the difficulty of (Ratio of inside building 100%) personnel mobilization Difficulty in speedy evacuation, especially in winter or rainy season, which may increase human casualties More human casualties may happen in office and commercial facilities, rather than in homes Weekday Noon (12:00 pm) There are a large number of people who have to stay in (Ratio of inside building 90%) offices and commercial facilities due to transportation problems caused by road and bridge damage Less human casualties than the other scenes Weekend Afternoon (18:00 pm) May cause delay on search and rescue due to the difficulty of (Ratio of inside building 70%) personnel mobilization Source: JICA Project Team
(2) Earthquake Occurrence Scenes for 2030
Two scenes for 2030, with and without strengthening of the total building seismic capacity (Figure 4.2.2), were considered to provide information for determination of risk reduction targets in disaster risk reduction and management plan. The scene of without strengthening of the building seismic performance (Extrapolation) means the structure type distribution of 2030 is the same as that of 2016. The scene with strengthening building seismic performance (Seismic Strengthening) changes the structure type distribution of 2016, having different cases with different assumptions. The effect of building seismic performance strengthening could be obtained by comparing the results from the two scenes of with and without. Since the purpose of risk assessment for 2030 is to understand the risk reduction effect, earthquake occurrence time at night is used as a reference. The contents of risk assessment for 2030 are building damage and life loss, and the damages of infrastructure and lifeline facilities were not estimated due to the difficulty and reliability of estimation on future inventories without sufficient information under current incomplete data set. The contents of risk assessment for earthquake occurrence scenes of 2016 and 2030 are described in Table 4.2.3.
Source: JICA project Team Figure 4.2.2 Earthquake occurrence scenes for 2030
4-9 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.2.3 Earthquake occurrence scenes and risk assessment contents Infrastructure Human Economic Year Scene Building / Lifeline Casualty Loss Night 2016 Weekday noon ○ ○ ○ ○ Weekend afternoon Extrapolation 2030 Seismic ○ - ○ - Strengthening Source: JICA project Team
4.3 Approach and Results of Seismic Risk Assessment
Structural damage of general buildings, schools, hospitals, government buildings, roads, bridges, water supply pipeline networks, sewage pipeline networks, power poles and mobile base transceiver stations (BTS) and human casualty and economic loss for each scenario earthquake ground motion was carried out. This approach of risk assessment is basically that commonly used in Japan and modified as necessary and availability of damage data from the Gorkha Earthquake. The approach was discussed in WG meetings as well as individually with counterparts and related organizations, so as to reach common consensus.
4.3.1 Damage to Buildings
Methodology of risk assessment
A flow of risk assessment is shown in Figure 4.3.1. Building inventory data (refer to Section 2.2.1 building), result of hazard assessment (PGA, refer to PR) and proposed damage function of building were used for the damage assessment of each 250m x 250m grid, and collected to estimate the risk of the building. Refer to Appendix 2.4 for the detail.
4-10 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA project Team Figure 4.3.1 Flow of building risk assessment
Proposed damage function for building
Proposed damage function for buildings is shown in Figure 4.3.2. Damage functions for typical (centre) area and perimeter area were used per the predominant period of the ground. This allocation was done based on the response analysis at each grid for average building period of 0.3 to 0.7sec. against predominant period of the ground. Category of the damage function and structural type is shown in Table 4.3.1. The damage function of building shows the damage ratio for general buildings, and the damage probability for public buildings, since the number is limited at each grid.
a) General (centre) area of the Valley b) Perimeter area of the Valley, suffix p (predominant period of the ground, (predominant period of the ground Tg > 1.5sec & Tg ≤0.3sec) 0.3sec< Tg ≤ 1.5sec) Source: JICA project Team Figure 4.3.2 Proposed Damage function for DG 4+ 5 at center area and perimeter area
4-11 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.1 Category of damage function and structural type
Centre (general) area of the Valley: predominant period of the ground, Tg > 1.5sec & Tg ≤0.3sec Perimeter area of the Valley, suffix p: predominant period of the ground 0.3sec< Tg ≤ 1.5sec Source: JICA Project Team
The method and procedure for the development of damage function is described in detail in Appendix 7. The damage function is created based on the building damage data of Gorkha earthquake, soil amplification characteristics of KV, seismic performance capability of buildings as well as the experience of developing damage function in Japan. The damage function of 2002 project is also used as a reference. The damage function is assumed as log-normal distribution. It is created for six structure categories and three damage levels. The parameters of damage function for central area of KV are shown in Table 4.3.2 and Table 4.3.3 is that for perimeter area. Table 4.3.2 Damage function for central area of KV
EMS DL4 & 5 EMS DL3, 4 & 5 EMS DL2, 3, 4 & 5 Standard Standard Standard Mean Mean Mean deviation deviation deviation Masonry 1 5.00 0.66 4.725 0.66 4.525 0.66 Masonry 2 5.40 0.62 5.10 0.62 4.70 0.62 Masonry 3 5.80 0.55 5.55 0.55 5.175 0.55 Masonry 4 6.07 0.55 5.85 0.55 5.55 0.55 RC 1 6.24 0.50 6.00 0.50 5.75 0.50 RC 2 6.40 0.50 6.15 0.50 5.90 0.50 Source: JICA Project Team
Table 4.3.3 Damage function for perimeter area of KV
EMS DL4 & 5 EMS DL3, 4 & 5 EMS DL2, 3, 4 & 5 Standard Standard Standard Mean Mean Mean deviation deviation deviation Masonry 1 4.86 0.66 4.585 0.66 4.385 0.66 Masonry 2 5.26 0.62 4.96 0.62 4.56 0.62 Masonry 3 5.66 0.55 5.41 0.55 5.035 0.55 Masonry 4 5.93 0.55 5.71 0.55 5.41 0.55 RC 1 6.10 0.50 5.86 0.50 5.61 0.50 RC 2 6.26 0.50 6.01 0.50 5.76 0.50 Source: JICA Project Team
4-12 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Damage function of historical building (monument) covers historical buildings (monuments) at Protected Monument Zone, PMZ at three Durbar Squares in KV. Damage ratio is as plotted in Figure 4.3.3 with supposed PGA by the 2015 Gorkha earthquake. Category of “Reconstruction” by damage data of DOA was considered as Damage Grade 4+ 5, and “Retrofitting/ Conservation” was considered as DG 2+3 of EMS 98 respectively. Refer to the damage data of the Appendix 7 of the report.
DG 4+5 DG 3+4+5 DG 2+3+4+5
100%
80%
60%
40%
20%
0% Damage ratio of buildings of ratio Damage 0 200 400 600 800
:DG 4+5 :DG 2+3+4+5 Source: JICA project Team Figure 4.3.3 Proposed damage function for historical buildings (monuments)
Assessment Results
Building damages estimated for each scenario ground motion are summarized in Table 4.3.4 and the details are given after.
4-13 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.4 Estimated building damage (Upper: damage number; Lower: damage ratio)
Scenario earthquake Ground Motion Category Damage Level WN CNS-1 CNS-2 CNS-3 Heavy damage 24,961 65,314 136,060 199,643 (EMS DL4&5) 5.6% 14.7% 30.6% 44.9% Buildings (2016) Moderate damage 21,967 42,940 62,691 67,418 (Total buildings 444,554) (EMS DL3) 4.9% 9.7% 14.1% 15.2% Slight damage 43,564 67,770 77,713 70,462 (EMS DL2) 9.8% 15.2% 17.5% 15.9% 33,763 88,681 185,796 273,269 Case-0 5.6% 14.6% 30.6% 45.1% 28,377 79,075 171,977 258,044 Case-1 4.7% 13.0% 28.4% 42.5% 13,627 56,452 146,361 234,477 Case-2 Buildings (2030, EMS 2.2% 9.3% 24.1% 38.7% DL4&5) 12,162 49,970 131,095 213,481 (Total buildings 606,506) Case-3 2.0% 8.2% 21.6% 35.2% 16,147 52,413 129,904 210,181 Case-4 2.7% 8.6% 21.4% 34.7% 11,138 41,230 111,854 189,357 Case-5 1.8% 6.8% 18.4% 31.2% 237 737 1,654 2,486 Heavy 4.1% 12.9% 28.9% 43.4% Schools 253 539 810 875 Moderate (Total buildings 5,731) 4.4% 9.4% 14.1% 15.3% 568 916 1,057 960 Slight 9.9% 16.0% 18.4% 16.8% 20 64 153 235 Heavy 3.4% 11.0% 26.2% 40.2% Health facilities 24 55 83 94 Moderate (Total buildings 584) 4.1% 9.4% 14.2% 16.1% 51 85 105 97 Slight 8.7% 14.6% 18.0% 16.6% 20 59 126 186 Heavy 4.2% 12.3% 26.4% 38.9% Government buildings 20 44 66 73 Moderate (Total buildings 478) 4.2% 9.2% 13.8% 15.3% 44 71 85 80 Slight 9.2% 14.9% 17.8% 16.7% Source: JICA Project Team
The estimated PGA by the four scenario earthquakes is 150 to 200gal for WN, 150 to 400gal for CNS-1, 250 to 600gal for CNS-2, and 300 to 800gal for CNS-3 generally. The epicenter is close in case of CNS and the PGA changes by the location. Ground motion expressed by MMI is VII for WN, VII to VIII for CNS-1, VIII (partially IX) for CNS-2 and VIII to IX for
4-14 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
CNS-3. PGA of the 2015 Gorkha Earthquake is estimated as 150 to 200gal, and the intensity VII by MMI. Refer to PR for the Hazard Assessment. Design PGA of IS 1893-1 (2002) is 180gal for information.
(1) Damage to general buildings
Results of the assessment by four cases of scenario earthquakes in 2016 and in 2030 are shown below. Damage ratio of building by the scenario earthquakes CNS-1 and CNS-2 have mainly been introduced.
Damage at year 2016
The total number of buildings has been estimated as 444,554 buildings. The ratio of each structural type is 47% for “RC non-engineered”, 28% for “Brick and stone masonry with cement mortar joint”, 14% for “Brick and stone masonry with mud mortar joint”, and 2% for “Others (RC engineered, and others)”. Result of risk assessment for general building in 2016 is shown in Figure 4.3.4 and Table 4.3.4 (Portion of 2016 of the table). Buildings with heavy structural damage (grade 4 and 5 by EMS 98) due to the scenario earthquake CNS-1 has been estimated as 65,314 buildings which is 14.7% of the total number of buildings.
Source: JICA Project Team Figure 4.3.4 Result of risk assessment for general building in year 2016 (CNS-2, Left:Damage number, Right:Damage ratio)
Damage at year 2030
The total number of building has been estimated as 606,506 buildings. The result of risk assessment for general buildings in 2030 is shown in Figure 4.3.5 and Table 4.3.4 (Portion of 2030 and without Building Seismic Performance Strengthening). Buildings with heavy structural damage (grade 4 and 5 by EMS 98) due to the scenario earthquake CNS-1 has been estimated as 88,681 buildings which is 14.6% of the total number of buildings, if no building seismic performance strengthening (Case0) is taken.
4-15 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA Project Team Figure 4.3.5 Result of risk assessment for general building in year 2030 (Without BSPS, CNS-2, Left: Damage number, Right: Damage ratio)
Result of estimated damage of building, heavy structural damage (damage grade 4+5) and moderate structural damage (damage grade 3), at year 2030 is shown in Table 4.3.7 where five cases of BSPS are taken. As far as the number and ratio of structural type, refer to “Section 2.2.1 Building”. Refer to separated “Map book” for the map of damage assessment for Case 1 to Case 5.
Building damage and corresponding economic losses are shown in Table 4.3.5 for the cases of with and without BSPS for CNS-1. Varying according to the cases, the number of heavily damaged buildings is reduced from 11% to 53% and the amount of loss is reduced from about 8% to 30%. The cost for the case of with BSPS, comparing with those of without BSPS, is shown in Table 4.3.6. For Case 1, which means the masonry building with cement mortar are constructed instead of mud mortar for new buildings from 2016 to 2030, the building cost will increase about 9%. In the case of Case 2, where all new and existing mud mortar masonry buildings are replaced by cement mortar buildings, the cost increases about 14%. On the other hand, in the Case 4 and Case 5, which change the majority of masonry buildings from both mud mortar and cement mortar, to engineered RC buildings, leads to a considerable increase of building cost due to the increase of building numbers and high unit construction cost.
4-16 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.5 Building Damage and Loss Amount Reduction in Reduction in Reduction in Buildings with number of Loss amount Case loss amount loss amount heavy damage damaged (mil. NPR) (mil. NPR) (%) buildings Case0 88,681 269,789 Case1 79,075 10.8% 247,974 21,815 8.1% Case2 56,452 36.3% 227,573 42,216 15.6% Case3 49,970 43.7% 187,832 81,956 30.4% Case4 52,413 40.9% 221,349 48,440 18.0% Case5 41,230 53.5% 207,520 62,269 23.1% Source: JICA Project Team
Table 4.3.6 Cost Increase of Cases with BSPS to the Case without BSPS Buildings to be Cost increase to Ratio of Cost Construction Case rebuilt or newly Case 0 (mil. increase cost increase per Cost (mil. NPR) constructed NPR) (%) year Case1 161,953 881,718 79,718 9% 5,694 Case2 243,047 1,072,452 148,391 14% 10,599 Case3 452,543 2,690,179 299,228 11% 21,373 Case4 293,099 2,263,313 894,945 40% 63,925 Case5 389,795 3,009,995 1,233,103 41% 88,079 Note: Building increase from 2016 to 2030 is 161,952; construction cost is 802,000 mil. NPR Cost increase per year is calculated for 14 years, from 2017 to 2030. Source: JICA Project Team
4-17 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.7 Results of building damage assessment (2016, 2030 and 2030 with 5 cases of BSPS)
Building Damage Scinario Eq. Heavily (DG4+5) Partly (DG3) 2016 WN 24,961 5.6% 21,967 4.9% CNS-1 65,314 14.7% 42,940 9.7% CNS-2 136,060 30.6% 62,691 14.1% Total Number: 444,554 CNS-3 199,643 44.9% 67,418 15.2%
Building Damage Scinario Eq. 2030 Heavily (DG4+5) Partly (DG3) without BSPS WN 33,763 5.6% 29,831 4.9% (※1) CNS-1 88,681 14.6% 58,470 9.6% CNS-2 185,796 30.6% 85,520 14.1% Total Number: 606,506 CNS-3 273,269 45.1% 91,892 15.2% Rate of Mitigating Damage Building Damage Scinario Eq. by BSPS (※2) 2030 Heavily (DG4+5) Partly (DG3) Heavily Partly with BSPS WN 28,377 4.7% 26,558 4.4% 16.0% 11.0% Case01 CNS-1 79,075 13.0% 55,103 9.1% 10.8% 5.8% CNS-2 171,977 28.4% 83,859 13.8% 7.4% 1.9% Total Number: 606,506 CNS-3 258,044 42.5% 92,321 15.2% 5.6% -0.5% Rate of Mitigating Damage Building Damage Scinario Eq. by BSPS (※2) 2030 Heavily (DG4+5) Partly (DG3) Heavily Partly with BSPS WN 13,627 2.2% 18,881 3.1% 59.6% 36.7% Case02 CNS-1 56,452 9.3% 50,010 8.2% 36.3% 14.5% CNS-2 146,361 24.1% 83,717 13.8% 21.2% 2.1% Total Number: 606,506 CNS-3 234,477 38.7% 95,133 15.7% 14.2% -3.5% Rate of Mitigating Damage Building Damage Scinario Eq. by BSPS (※2) 2030 Heavily (DG4+5) Partly (DG3) Heavily Partly with BSPS WN 12,162 2.0% 16,590 2.7% 64.0% 44.4% Case03 CNS-1 49,970 8.2% 45,067 7.4% 43.7% 22.9% CNS-2 131,095 21.6% 78,997 13.0% 29.4% 7.6% Total Number: 606,506 CNS-3 213,481 35.2% 93,462 15.4% 21.9% -1.7% Rate of Mitigating Damage Building Damage Scinario Eq. by BSPS (※2) 2030 Heavily (DG4+5) Partly (DG3) Heavily Partly with BSPS WN 16,147 2.7% 17,900 3.0% 52.2% 40.0% Case04 CNS-1 52,413 8.6% 45,293 7.5% 40.9% 22.5% CNS-2 129,904 21.4% 79,695 13.1% 30.1% 6.8% Total Number: 606,506 CNS-3 210,181 34.7% 95,190 15.7% 23.1% -3.6% Rate of Mitigating Damage Building Damage Scinario Eq. by BSPS (※2) 2030 Heavily (DG4+5) Partly (DG3) Heavily Partly with BSPS WN 11,138 1.8% 14,252 2.3% 67.0% 52.2% Case05 CNS-1 41,230 6.8% 40,974 6.8% 53.5% 29.9% CNS-2 111,854 18.4% 77,650 12.8% 39.8% 9.2% Total Number: 606,506 CNS-3 189,357 31.2% 96,203 15.9% 30.7% -4.7% ※1: BSPS: Promotion on Building Seismic Performance Strengthening ※2: Rate of Mitigating Damage by BSPS Compared to 2030 Without BSPS Source: JICA Project Team
4-18 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Damage ratio of buildings with heavy and very heavy structural damage for 2030 without BSPS and for 5 cases of 2030 with BSPS is shown at upper portion of Figure 4.3.6. Ratio of damage reduction against 2030 without BSPS is shown at lower portion of the same figure. The special case that all buildings are “RC engineered” in 2030 is shown as “2030 with BSPS (Reference)” at the right end for reference.
Rate of Damage Building (Heavily: DG4+5)
50% WN CNS-1 CNS-2 CNS-3 45.1% 42.5% Rate(%) = Damage Building Number / Total Building Number 38.7% 40% 35.2% 34.7% 30.6% 31.2% 30% 28.4% 24.1% 25.4% 21.6% 21.4% 20% 18.4% 14.6% 13.0% 13.5% 9.3% 8.2% 8.6% 10% 6.8%
Building Damage Rate (Heavily) 5.6% 4.7% 3.8% 2.2% 2.0% 2.7% 1.8% 0.5% 0% Without BSPS With BSPS Case01 With BSPS Case02 With BSPS Case03 With BSPS Case04 With BSPS Case05 With BSPS (Reference) 2030 with/without BSPS
Rate of Mitigating Damage by BSPS Compared to Without BSPS (Heavily: DG4+5)
100% WN CNS-1 CNS-2 CNS-3 90.1%
Rate(%) = (Without BSPS - With BSPS) / Without BSPS 80% 73.9% 67.0% 64.0% 59.6% 60% 56.0% 52.2% 53.5% 43.7% 43.7% 40.9% 39.8% 40% 36.3% 29.4% 30.1% 30.7% 21.2% 21.9% 23.1% 20% 16.0% 14.2% 10.8% 7.4%5.6%
Rate of Mitigating Damage by BSPS (Heavily) BSPS by Damage Mitigating of Rate 0% With BSPS Case01 With BSPS Case02 With BSPS Case03 With BSPS Case04 With BSPS Case05 With BSPS (Reference) Cases of With BSPS Source: JICA Project Team Figure 4.3.6 Building damage ratio (DG 4+5) in 2030 with BSPS (upper) and ratio of damage reduction (lower) by scenario earthquakes
(2) Damage to school buildings
The total number of school buildings was estimated as 5,731 buildings based on a school building inventory as of 2016 and the result of damage assessment is shown in Figure 4.3.7, Table 4.3.8 and Figure 4.3.8. The ratio of each structural type are, 61% for “Masonry”, 30% for “RC non-engineered”, and 9% for “RC engineered”, and more than half are masonry structure. Buildings with heavy and very heavy structural damage (grade 4 and 5 by EMS 98) due to the scenario earthquake CNS-1 and CNS-2 were estimated as 737 buildings (12.9%) and 1,654 buildings (28.9%) respectively, which is the total of the damage probability of each building. The probability of damage grade 4+5 for each building against the CNS-2 is shown in Figure 4.3.7. Refer to the “Map Book” for more detail.
4-19 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA Project Team Figure 4.3.7 Result of risk assessment for school building (distribution)
Table 4.3.8 Result of risk assessment for school buildings (number and ratio) Scenario Damage Level Total (5,731) Earthquaake DL2 DL3 DL 4 & 5 WN 568 253 237 1,058 18.5% CNS-1 916 539 737 2,192 38.2% CNS-2 1,057 810 1,654 3,521 61.4% CNS-3 960 875 2,486 4,321 75.4% Source: JICA Project Team
Source: JICA Project Team Figure 4.3.8 Result of risk assessment for school buildings (number and ratio)
It was noted the seismic risk assessment (3 damage degrees) for 1,200 public school buildings based on NBC was carried out by MOE and ADB (National Workshop on School Earthquake Safety, 2010) and the results are outlined in Attachment 2.5.1.
4-20 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
(3) Damage to health facility buildings
The total number of health facilities has been estimated as 584 buildings based on a health facility inventory as of 2016 and the result of damage assessment is shown in Figure 4.3.9, Table 4.3.9 and Figure 4.3.10. The ratio of each structural type is, 40% for “Masonry”, 51% for “RC non-engineered”, and 9% for “RC engineered”. Buildings with heavy and very heavy structural damage (grade 4 and 5 by EMS 98) due to the scenario earthquake CNS-1 and CNS-2 is 64 buildings (11.0%) and 153 buildings (26.2%) respectively, which is the total of the damage probability of each building. Ratio of buildings with moderate damage (grade 3) and more are 20.4% and 40.4% respectively. The probability of damage grade 4+5 for each building against the CNS-2 is shown in Figure 4.3.9. Refer to the “Map Book” for details of each building.
Source: JICA Project Team Figure 4.3.9 Result of risk assessment for medical facilities (distribution)
Table 4.3.9 Result of risk assessment for medical facilities (number and ratio) Scenario Damage Level Total (584) Earthquaake DL2 DL3 DL 4 & 5 WN 51 24 20 95 16.3% CNS-1 85 55 64 204 34.9% CNS-2 105 83 153 341 58.4% CNS-3 97 94 235 426 72.9% Source: JICA Project Team
4-21 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA Project Team Figure 4.3.10 Result of risk assessment for medical facilities (number and ratio)
The results of seismic assessment for 28 buildings of fourteen important hospitals based on the four level earthquakes are described in the report by WHO (A Structural Vulnerability Assessment of Hospitals in Kathmandu Valley, 2002). Please refer to Attachment 2.5.1 for the summary.
(4) Damage to governmental buildings
The total number of governmental buildings has been estimated as 478 buildings based on a governmental building inventory as of 2016 and the result of damage assessment is shown in Figure 4.3.11, Table 4.3.10 and Figure 4.3.12. The ratio of each structural type is, 50% for “Masonry”, 2% for “RC non-engineered”, and 48% for “RC engineered” (including buildings not compatible to the latest seismic code). Buildings with heavy and very heavy structural damage (grade 4 and 5 by EMS 98) due to the scenario earthquake CNS-1 and CNS-2 is 59 buildings (12.3%) and 126 buildings (26.4%) respectively, which is the total of the damage probability of each building. Ratio of buildings with moderate damage (grade 3) and more are 20.4% and 40.4% respectively. The probability of damage grade 4+5 for each building against the CNS-2 is shown in Figure 4.3.11. Refer to the “Map Book” for more detail.
4-22 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA Project Team Figure 4.3.11 Result of risk assessment for governmental buildings (distribution)
Table 4.3.10 Result of risk assessment for governmental buildings (number and ratio) Scenario Damage Level Total (478) Earthquaake DL2 DL3 DL 4 & 5 WN 44 20 20 84 17.6% CNS-1 71 44 59 174 36.4% CNS-2 85 66 126 277 57.9% CNS-3 80 73 186 339 70.9% Source: JICA Project Team
Source: JICA Project Team Figure 4.3.12 Result of risk assessment for governmental buildings (number and ratio)
(5) Damage to historical buildings (monuments)
Risk assessment of historical buildings (monuments) at the Protected Monument Zone (PMZ) of 3 Durbar Squares (Hanumandhoka, Patan, and Bhaktapur World Heritage Site
4-23 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
(WHS)) was done. In the 2015 Gorkha Earthquake, 21 monuments (19.4%) were damaged with “Reconstruction” which is equivalent to damage grade 4+5, and 28 monuments (25.95) were damaged with “Retrofitting/Conservation), which is equivalent to damage grade 2+3, and total 49 out of 108 monuments were damaged (45.4%). Damage data at each Durbar Square is shown in Attachment 2.4.1. There is no big difference of PGA at 3 Durbar Squares by each scenario earthquake, and damage number and damage ratio is shown in Figure 4.3.13. Ratio of monument with damage grade 4+5 by scenario earthquake the CNS-1 and the CNS-2 is 50.4% and 72.8% respectively.
Damge Grade 4+5 Dmage Grade 2+3 120 95.9% 91.4% 100 12.3 78.0% 80 18.7
60 52.1% 27.6
40 28.1 20
Damage of Monument (Nos) (Nos) Monument of Damage 23.9 50.4 72.8 83.6 0 WN1234 CNS-1 CNS-2 CNS-3 Scenario Earthquake
Source: JICA Project Team Figure 4.3.13 Result of risk assessment of historical building at 3 Durbar Squares (total 108 monuments)
According to the following reference, PGA 0.266g (g: gravity acceleration) was estimated as a hazard assessment with 475 years return period at the bedrock of the area of Durbar Square. As far as the response at the free surface ground, 0.542g (linear response) and 0.22g (equivalent linear response) at Kathmandu Durbar Square, 0.425g (linear response) and 0.314g (equivalent linear response) at Lalitpur (Patan) Durbar Square, 0.32g (linear response) and 0.208g (equivalent linear response) at Bhaktapur Durbar Square are shown. The result of equivalent linear response is relatively similar to that of the scenario earthquake CNS-1.
Reference: Prem Nam Maskey, 2014 ,“Seismic Hazard at the Monument Zones of the World Heritage Site of the Kathmandu Valley”, CARD, IOE, Tribhuvan University, Nepal
(6) Risk assessment by a detail building survey
Drawings, material strength and weight data for a total of eleven buildings were collected through a detail building survey during the 1st year of the project, and the seismic capacity for both directions (x and y directions) of each building were assessed using a simplified
4-24 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
method. The summary is shown in Figure 4.3.14. Effective PGA (cm/sec2, gal) causing damage grade 4 (heavy damage) and more by EMS-98 was estimated, and the estimated PGA has some range considering the variety of building response by earthquake waves and reliable accuracy of the assessment (shown by blue color dotted line). Data of detail assessment is shown the Attachment 2.5.1 of the report.
Effective PGA of the range of 200~300gal was supposed at the center area in the valley by the scenario earthquake CNS-1, and the area is shown in red color dotted line. Similarly the range of 250~400gal was supposed by the scenario earthquake CNS-2 and the area is shown in red color dotted line for essential public buildings such as school, hospital and governmental building.
It is to be noted at the south side of the valley where PGA of 150~400gal by CNS-1 and PGA of 250~600gal by the CNS-2 have been estimated, but buildings there were supposed to be located at the center area in the valley.
700
Estimated PGA causing
heavy damage of bldg
600
RC eng. (x, y) (x, eng. RC Residential
eng. (x, y)
RC eng. (x, y) (x, eng. RC Hospital
- 500
RC non RC Residential RC eng.. (x, y) (x, eng.. RC Localgovernmental
Brick/cement(x, y) Residential 400 rise resi.
RC eng. (x, y) (x, eng. RC Governmental
Brick/mud (x, y) Residential Brick/cement(x, y) School
PGA by CNS-2 y) (x, eng. RC High - PGA by CNS-2
Brick/mud (x, y) Historical At the center 300 at the center causing heavy structural damage structural heavy causing
Adobe(x, y) Residential PGA by CNS-1 PGA by CNS-1
(gal) (gal) at the center At the center 200 Design PGA by IS 1893-1 (2002) Observed PGA by 2015 Gorkha Earthquake 100 Estimated PGA PGA Estimated
0 Essential pub. bldg Essential pub. bldg 0 1 2 2 4 3 6 7 8 9 10 4 12 5 14 6 16 8 1810 1120 22 Masonry RC Note: Usage, 1~6: Residential, 7 : School, 8 : Hospital, 9 : Historical building, 10 & 11 : Governmental building Structural type, 1: Adobe, 2 & 9 : Brick masonry with mud mortar joint, 3 & 7 : Brick masonry with cement mortar joint, 4 : Non-engineered RC, 5, 6, 8, 10 and 11 : Engineered RC (including building designed by old design code), photos (1~9) by ESS, photos 10, 11 by the project Source: JICA Project Team, photos(1~8)by ESS Figure 4.3.14 Risk assessment of eleven buildings based on detailed building survey
(PGA causing heavy damage was assessed by a simplified method, and the range of PGA by scenario earthquakes at the center in the valley was added.)
4-25 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Observed PGA due to the 2015 Gorkha Earthquake (average 147gal, area of green color dotted line) is shown in the Figure, and design PGA by IS 1893-1 (2002) is shown using yellow color dotted line in the Figure. According to IS 1893-1 (2002), design PGA is 180gal which is the half of 360gal, possible maximum ground acceleration. PGA for design is not shown explicitly in NBC 105. Importance factor 1.5 for essential public buildings such as monument, school, hospital and governmental building is indicated at design codes.
Damage due to the 2015 Gorkha Earthquake
Damage grade 3 of EMS-98 for building No. 1 (Adobe), DG 3 for building No. 2 (Brick masonry with mud mortar joint), DG 3 for building No.4 (RC with soft story), and DG 2 for building No. 6 (High-rise RC) were recorded respectively. No damages for other buildings. Buildings are supposed after the repair work or no damage in the following assessment.
The result of risk assessment of each building which location is supposed at the center of the valley is as follows.
1) Design PGA
A high probability of DG 4 and more has been estimated for building No.1 (Adobe), building No.2 (Brick masonry with mud mortar joint), and building No. 4 (RC with a soft story) for general buildings (against 180gal). A high probability of DG 4 and more has been estimated for building No.7 (Brick masonry with cement mortar, school), building No.9 (Brick masonry with mud mortar joint, historical building) for essential buildings (against 270gal= 180 gal x 1.5).
2) Scenario earthquake CNS-1
A high probability of DG 4 and more has been estimated for building No.3 (Brick masonry with cement mortar joint), and building No. 6 (Engineered hi-rise RC) within general buildings.
3) Scenario earthquake CNS-2
A high probability of DG 4 and more has been estimated for building No.10 (RC, governmental building), and building No.11 (RC, local governmental (municipality) building) for essential buildings. Some probability of DG 4 has been estimated for recently constructed building No. 8 (Engineered RC, hospital) is estimated. Low probability of DG 4 has been estimated for building No. 5 (Engineered RC), and also, it has been estimated that the building was designed with enough allowance, but will suffer DG 3 (brick wall damage).
(7) A study on the result of risk assessment
Estimated PGA by 4 scenario earthquakes is 150 to 200gal for WN, 150 to 400gal for CNS-1,
4-26 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
250 to 600gal for CNS-2, and 300 to 800gal for CNS-3 generally. The epicenter is close in case of CNS and the PGA changes by the location. In the above (6), estimated PGA was supposed as the effective PGA for buildings and the damage risk of each building has been evaluated by supposing that the location is the center of the valley. The damage situation for overall valley has been studied here.
1) Scenario earthquake WN
The size of ground motion for the scenario earthquake WN is similar to that of the 2015 Gorkha Earthquake (PGA is 150gal at the center and 200gal at the perimeter is estimated generally), and the damage degree is also similar.
2) Scenario earthquake CNS-1
The ratio of damage level 4 (heavy structural damage) and more is 15% for general buildings in the case of CNS-1. The ratio of damage level 4 and more is more than 10% for schools, medical facilities and government buildings. PGA at the center of the valley is 200 to 300gal, and the ratio of damage grade 4 and more is approximately 20% for brick masonry with cement mortar, and more than 10% for RC building. Furthermore the damage at the south area of the valley results in almost destruction for brick masonry with mud mortar, damage ratio of heavy structural damage is 40 to 50% for brick masonry with cement mortar, and is 20 to 30% for RC building generally. On the other hand, the damage at the north area of the valley has been estimated to the degree of slightly bigger than that of the 2015 Gorkha Earthquake. The following two points are characteristics.
■ There is a high possibility of damage grade 4 and more at the center and south area of the valley for public buildings; those are the essential facilities for disaster mitigation. ■ There is a big difference of damage situation between the south and north area of the valley.
3) Scenario earthquake CNS-2
The size of ground motion of the CNS-2 is generally 1.5 times of that of the CNS-1. Estimated PGA at the center area is 250 to 400gal, and big damage for masonry mainly including RC building is assessed. Estimated PGA at the south area is 400 to 600gal, and almost destruction for masonry including cement mortar type, and 40 to 60% for RC building. In addition, the damage at the north area has been assessed.
4) Scenario earthquake CNS-3
The ground motion by the CNS-3 might be overestimated, and the assessed damage also might be overestimated.
4-27 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
4.3.2 Damage to transportation Infrastructure
(1) Damage to roads
The purpose of risk assessment on roads is to identify the degree of traffic disturbance induced by an earthquake. Landslides, liquefaction and blockage caused by collapsed building debris were the major reasons of road damage.
1) Damage by slope failure and liquefaction:
Roads located in high potential areas of slope failure and liquefaction are at risk of traffic disturbance by severed roads due to sediment, subsidence of road base and fluctuation of road structures. In this project, possible damaged roads due to slope failure and liquefaction were identified by spatially comparing the 250m mesh-grid potential maps of slope failure and liquefaction based on the scenario seismic ground motion and the map of existing road network for each grid. The road lengths that have been determined to have high potential of slope failure and liquefaction for each scenario earthquake is shown in Table 4.3.11.
4-28 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.11 Summary of lengths of possible damaged road by slope failure and liquefaction
(Km) Slope Failure: High Liquefaction: High Road Class Total Length WN CNS-1 CNS-2 CNS-3 WN CNS-1 CNS-2 CNS-3 NH: National Highway 0.0 0.0 1.3 2.8 0.0 1.8 7.6 11.2 61.6 (DOR)
FRN: Feeder Roads Major 0.0 0.9 8.8 25.1 0.0 5.6 12.8 19.3 290.1 (DOR)
FRO:Feeder Roads Minor 0.0 0.0 1.0 2.4 0.0 0.0 0.1 0.9 84.7 (DOR)
SUR: Strategic Urban Road 0.0 0.0 1.6 4.6 0.0 4.3 13.6 18.3 129.0 (DOR)
District or Village Road 0.0 5.7 85.8 355.6 0.0 64.3 241.0 405.6 5,245.2
Total 0.0 6.6 98.5 390.6 0.0 76.1 274.9 455.3 5,810.6
Remark: There is no Mid-Hill Road (MH) and Postal Road (PR) in the study area based on the annual statistic document preparead by DOR Source: JICA Project Team
2) Road blockage of the narrow streets:
There is a risk of road blockage due to the debris of collapsed buildings by earthquake ground motion in the relatively narrow streets. In this project, the risk of traffic disturbance after the earthquake due to the collapse of adjacent buildings was estimated to calculate the road link blockage rate for each grid targeting the central city where narrow streets are concentrated. First, the road link blockage rate by road segment was calculated based on a typical road width of objective road and a building damage rate of the grid where the road segment is located. Then, the average of road link blockage rates weighted by length of road segments was calculated in each grid.
① Width of representative road: Less than 3.5m Rate of road-link blockage (%)=0.9009×Rate of damaged building(%)+19.845
② Width of representative road: 3.5m to < 5.5m Rate of road-link blockage (%)=0.3514×Rate of damaged building(%)+13.189
③ Width of representative road: 5.5m to < 13m Rate of road-link blockage (%)= 0.2229 ×Rate of damaged building(%)-1.5026 Rate of damaged building = (Rate of total damage + ½×Rate of partial damage) Source: Central Disaster Prevention Council, Japan
The typical road width is set based on the road network data which has the attribute data for road widths, and the building damage rate is the sum of half of the complete destruction rate (DL4+5) and half of partial destruction rate (DL3) for each grid. Figure 4.3.15 shows the result of evaluated road link blockage for the emergency transportation road network proposed by JICA RRNE project.
4-29 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Source: JICA project team Figure 4.3.15 Distribution of road link blockage rates on the ETRN due to the debris of collapsed buildings by scenario earthquake (CNS-2)
From the distribution of potential slope failure based on scenario earthquakes, it can be observed that the potential of slope failure in the mountain area surrounding the Kathmandu Valley is higher than one in the central area, mainly depending on the distance from the epicenter of scenario earthquakes. There are high potential areas for liquefaction along alluvial plain and valley plain formed by tributary stream classified in the geomorphological map.
From distribution of road link blockage rate by grid, in overall, the rates are larger in the south area and mainly depending on the difference in strength of seismic motion due to the distance from the epicenter. In the distribution, when focused only on emergency transportation road network, the rates depend on the width of roads.
(2) Damage to bridges
1) Evaluation Criteria
If the earthquake risk of a bridge is regarded as the functional disturbance of road that can be represented by the event-related collapse of the substructure and the possibility of falling down of the superstructure. In this project, the possibility of collapse of the bridge substructure was evaluated utilizing "Response ductility factor" caused by the earthquake motion, and it was designated as an important index for quantitatively evaluating the risk of the earthquake at the bridge.
The liquefaction of the ground affects the stability of the structure at the time of earthquake, but the liquefying possibility evaluation in this project is based on insufficient soil test of each bridge site. In order to investigate the possibility of liquefaction with sufficient accuracy, it is necessary to carry out at least a boring test including grain size measurement of each layer at bridge site and PS logging, if possible. In this project, no borehole tests were
4-30 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
performed for individual bridge sites, though data were supplemented by referring to the micro topography classification and as such. From the viewpoint of micro topography classification, there is no difference in liquefaction potential for each bridge site, as almost all of the bridges in the Kathmandu Valley are located at the crossing of the river. Considering this point, the influence of liquefaction potential was not considered in the damage assessment. If the influence of liquefaction potential is going to be strongly utilized, the evaluation score of all the bridges may be lowered under investigation by one level.
2) Response Ductility Factor
Collapse of the bridge piers is the most noticeable failure mode for multiple span bridges and this event can be evaluated with "Response ductility factor" as the most quantitative indicator. Deformation of pier may exceed the elastic region when receiving seismic motion designated as the scenario earthquake in this project. "Response ductility factor" (how many times the deformation of elastic limit deformation occurs) is an effective index for evaluating the collapse process of the structure. In order to investigate the ductility factor, it is necessary to analyze the force-deformation relation of the structure from elastic region to plastic region. A general drawing and a reinforcing bar arrangement are necessary to fulfill the purpose of analysis. However, these drawings were not available for almost all the bridges within the survey area. Therefore, a kind of estimation method that is discussed in “Seismic Assessment Tool for Urgent Response and Notification” (the technical note of the National Institute for Land and Infrastructure Management No.71) was utilized. A kind of parameter designated as “yield seismic intensity” mentioned in this document can be utilized. Also, the force-deformation relation of the pier can be estimated using this parameter. As even this estimation method requires the data of external shape and dimensions of each objective, so general drawings of each bridge were created by on-site survey and the drawings were utilized. The response ductility factor was analyzed using the force-deformation relation of the pier and the response spectral value obtained from the result of WG 1. The response spectrum has been assigned for each 250m grid from the result of the hazard assessment within the study area. A regression formula was proposed in the research report of NILIM No.71, based on the investigations on the bridges of Japan, where the response of pier could be estimated from external dimensions of structure without detail information of rebar arrangement. Applying this regression analysis, it is needed to consider the fact that, this regression formula was established for different seismic design standard used in different time period. Regression analysis was made for four periods according to the design concept of the earthquake resistant design code. The coefficients of regression equation: α, β, γ are shown in Table 4.3.13.
4-31 The Project for Assessment of Earthquake Disaster Risk for the Kathmandu Valley in Nepal Final Report (Main Report)
Table 4.3.12 Four period classification of earthquake resistant design code applied for regression analysis Period of earthquake Policy features resistant design code A. before 1980 Linear analysis utilizing response spectra Preliminary study on non-linear analysis was applied considering ductile B. 1980~1990 behavior Large scale earthquake motion (about 400 Years of return period) C. 1990~1995 considering ductile behavior D. after 1995 Special code utilizing the experience of 1995 Hanshin Awaji Earthquake Source: JICA Project Team
The bridges to be surveyed do not conform to the earthquake resistant design code of Japan. Therefore a group of regression formulas were chosen in accordance with a classification closer to the policy of the target bridge design from the above classification with reference to construction year and construction agency etc. However, there are few cases where it can be judged that the bridges targeted this time are provided with reinforcing bars that take into consideration the ductility design concept. Even in the case where the reinforcing bars are designed with consideration of the ductility design concept, it is classified as “A” because it cannot be classified after “B” because there are problems in terms of detailed structure. (for instance, without movement limiting device in bearing)
Technical note of NILIM No.71 set the indicator called the yield seismic intensity Khy to estimate the force-deformation relation of the piers, and the regression equation to obtain this is expressed by the following equations 4.3.1) to 4.3.3) It is defined as follows.