Shallow crustal earthquake models, damage, and loss predictions in Banda Aceh, Indonesia

The city of Banda Aceh stands on Holocene fluvial basin sediment, surrounded by the Aceh and Seulimeum fault segments where large magnitude earthquakes can occur at any time. Such earthquakes could cause extensive physical infrastructure damages, injuries, and economic loss. This research aims to produce several earthquake scenarios, to determine the damage ratio of the buildings and its distribution in these different scenarios and to estimate the number of potential casualties and economic loss. Data analysis in this research includes modelling an earthquake by applying a ground motion model for shallow crustal earthquakes to yield an intensity map after the correction for site effects. The damage ratios for different types of buildings were calculated using the fragility curves of buildings that were developed by other researchers. Building occupants at different times, building damage ratios, and injury ratios were used to determine the number of injured residents in the earthquake scenario. The results show that Banda Aceh could be potentially experience earthquake ground motion within the intensity range from MMI VII to IX due to earthquakes of Mw 6.5 to Mw 7.0. The Mw 7.0 earthquake on the Aceh segment is the worst case scenario, causing building damage throughout Banda Aceh with casualty rates of 3.5–20% of the population in buildings in general if it happened during daytime, with economic loss about 3320 million USD. Mitigation and preparedness programs for the Banda Aceh community are therefore very important in order to reduce the potential level of damage, personal injury, and economic loss due to future earthquakes.


Introduction
The ability to reliably estimate earthquake damage is critical in evaluating building fragility, and forecasting damage and casualties due to a potential earthquake that could happen at any time and at any magnitude. Related to this, the vulnerability assessment of urban centres has been introduced by Baiocchi et al. (2012), which includes the vulnerability assessment of public infrastructure, private investment, and residential construction. The investigation of the building vulnerability, whether public infrastructure or residential construction, should also consider the mapping of the groundwater and high pore pressure sediment below the surface for the possibility of the sediment liquefaction and potential site amplification (Wood and Steedman 1992).
Assessment of vulnerability with respect to seismic activity for Indonesian cities like Banda Aceh is a concern not only because Indonesia lies on the ring of fire where seismic activity is high, but also because Banda Aceh stands on Holocene fluvial basin sediment and is surrounded by the Aceh and Seulimeum fault segments. The Aceh and Seulimum fault segments are part of the Great Sumatra fault (GSF). According to Ito et al. (2012), the northern part of the GSF contains a 200 kmlong segment that shows no seismic activity in at least 170 years. However, the northern portion of this fault segment, aseismic creep occurs from a 7.3 ± 4.8 km depth at the rate of 2.0 ± 0.6 cm/year. In the southwest part of the segment, Ito et al. (2012) estimated a locking depth of 14.8 ± 3.4 km, with a down-dip slip of 1.6 ± 0.6 cm/year. Those segments of the faults are capable of producing a Mw 7.0 earthquake that could put Banda Aceh in the high risk zone. Hence, the study of the site characteristics, which involves the generation of a seismic hazard map, and simulating earthquake loss prediction, are crucial for Banda Aceh. The Earthquake loss prediction is a method used to estimated damages, including infrastructure, human injuries, and economic loss under different earthquake scenarios (ground shaking and liquefaction) as well as the types of building which cause human casualties (Miura et al. 2008;Hashemi and Alesheikh 2011;Karimzadeh et al. 2014;Rusydy et al. 2017a;Guettiche et al. 2017;Saputra et al. 2017;Karimzadeh et al. 2017;Rusydy et al. 2018a;Rusydy et al. 2018b;Liu et al. 2019).
This paper aims to estimate the distribution and ratio of building damage as the results of several earthquake scenarios with different intensities, magnitudes and accounting for site amplification, and liquefaction effects. This paper also aims to predict the casualty numbers and economic loss in each earthquake scenario. In this study, the modelling of the ground shaking intensity resulting from different earthquakes scenarios were generated, building vulnerability depend on the fragility curve and liquefaction effect was computed, and the number of casualties in Banda Aceh estimated. The times of occurrence of earthquakes modelled in this research were 2 AM for night-time and 2 PM for day-time scenarios. The result of this study will contribute to enhancing earthquake disaster mitigation in the city. This paper will also highlight the parts of Banda Aceh that could suffer the most damage and economic loss due to ground shaking and associated liquefaction from potential future earthquakes.

Tectonics and seismicity of the study area
The great Aceh-Andaman Mw 9.1 earthquake of December 26, 2004 and the associated tragedy in Aceh province brought the active tectonics of this area into global focus (McCaffrey 2009). In the western part of Sumatra, a subduction zone is present between the Indo-Australian Plate and the Eurasian Plate. The Indo-Australian plate moves northward and is subducted beneath the Eurasia plate at approximately 5 cm/year (McCaffrey 2009). The oblique convergence system between the Indo-Australia plates and the Eurasian plate establish the Great Sumatra Fault (GSF), a major rightlateral trench fault system 1900 km in length (McCaffrey 2009;Sieh and Natawidjaja 2000).
The GSF completely traverses Sumatra Island, controlling much of the island's tectonics and has great potential to produce a large magnitude earthquakes (Sieh and Natawidjaja 2000;Chlieh et al. 2008;Muksin et al. 2018). Sieh and Natawidjaja (2000) divided the GSF into 19 major segments. In details, the 17 segments of the GSF are less than 200 km in length, with the other 2 being greater than 200 km. The two longest segments are the Aceh segment (200 km) and Renun segment (220 km).
Furthermore, Sieh and Natawidjaja (2000) defined the GSF in Aceh province as four major segments including: the Tripa segment (180 km-long), Aceh segment (200 km-long), Seulimeum segment (120 km-long) and Batee segment (150 km-long). The Aceh and Seulimeum segments are the closest segments to Banda Aceh. The Seulimeum segment experienced a magnitude M 7.2 earthquake in 1936 (Muksin et al. 2018). In 1965 another M 6.5 earthquake destroyed Krueng Raya, which is located closer to the Seulimeum segment (Sieh and Natawidjaja 2000). There were no major events recorded in the last 170 years on the Aceh segment. This earthquake history with few events in recent years suggests the potential for a large seismic event to happen in this region in the future. Ito et al. (2012) estimated that the Aceh segment could produce a Mw 7 earthquake, and endanger more than 260,000 residents of Banda Aceh.
Banda Aceh is bounded to the south-west by the Aceh segment and to the northeast by the Seulimeum segment, as shown in Fig. 1. From the west of coastal range moving south, the surface geology bounding the Aceh segment consists of pre-tertiary limestone. According to Barber and Crow (2005), this limestone formed in late Jurassic to early Cretaceous from a portion of Gondwanaland. Volcanic foothills, of Pleistocene to Holocene ages, are located on the northern coast of both sides of the Seulimeum segment. The alluvial deposit closest to the coastal area of Banda Aceh extends to a depth of more than 200 m and becomes thinner upstream (Bennet et al. 1981).
Banda Aceh is located on the Krueng Aceh basin, which was formed in response to fluvial processes. According to Bennet et al. (1981), the area below Banda Aceh city was formed during the Pleistocene and Holocene from undifferentiated alluvium and coastal sediments which consist of gravels, sands, and muds. Culshaw et al. (1979) noted that ancient sand-bars run sub-parallel to the coastline of Banda Aceh city. These sand-bars consist of brown to grey moderately weathered fine to medium sized grains with some silty sands. The lithology of the sand-bars are very similar to new beach sand. Having sand-bar formations along with Holocene sediments makes Banda Aceh a likely liquefaction-prone zone; furthermore, it will amplify the ground shaking of the earthquake. Banda Aceh city covers 9 sub-districts and 90 villages having different areas and populations; according to local authorities, 261,034 people live in Banda Aceh city as of September 2017.

Methodology
The tectonic and geological data obtained from Bennet et al. (1981), Sieh and Natawidjaja (2000) were digitized using ArcGIS software and overlaid with a SRTM map from NASA. First, the active fault location and magnitude estimation had to be determined from earlier studies prior to developing earthquake scenarios. The earthquake model applied to the various scenarios was well documented and based on previous studies of local and regional tectonics from Bennet et al. (1981), Sieh and Natawidjaja (2000), and Ito et al. (2012). The ground shaking map was developed using Intensity Prediction Equations (IPEs) as designed by Allen et al. (2012). After site coefficient correction, the IPEs equation was used to produce an intensity map at the levels of the bedrock and surface. Damage estimates of several types of buildings were determined using building fragility curves. As there were no published building fragility curves specific to Indonesia, this study used fragility curves developed for the Philippines by Tingatinga et al. (2013) and Naguit et al. (2017). Using the Philippines's fragility curve was appropriate because the building construction types in Indonesia are quite similar to those of the Philippines. The detailed methodology including data preparation, earthquake scenario development, building damage assessment analysis, casualty and economy loss estimation, are shown in Fig. 2.

Data preparation
GIS analysis for earthquake damage prediction requires shape-file data to be input into the GIS software. The tectonic and geological data were obtained from Bennet et al. (1981), Sieh and Natawidjaja (2000), and Barber and Crow (2005). The 2011 Banda Aceh building data were received from the Development Plan Department of Banda Aceh city. The 2016 Banda Aceh city building data using the GIS tool were then updated and digitized.

Earthquake scenario
The understanding of seismology has developed significantly over the past decades, especially the study and modelling of earthquake-generated ground motion. Empirical ground motion models have been developed that, give an earthquake's magnitude and distance, provide numerical estimates of ground motion intensity, in terms of both physical indicators like Peak Ground Acceleration (PGA), as well as macroseismic intensities, or Modified Mercalli Intensity (MMI) (Rusydy et al. 2017a;Rusydy et al. 2018a;Rusydy et al. 2018b). The latter are referred to as IPEs (Intensity Prediction Equations), and are particularly useful for rapidly and simply communicating earthquake hazard and risk to the media and public (Allen et al. 2012;Rusydy et al. 2017a;Rusydy et al. 2018a).
Different tectonic settings have different applicable ground motion models, but the IPEs developed by Allen et al. (2012) were specifically designed for a wide range of tectonic environments. Also, the Allen et al. (2012) IPEs were designed to be applicable for distance, depth and magnitude ranges that match well the earthquake scenarios of our study: distance 2.5-15 km, depth < 20 km, and 5.0 < Mw < 7.9. For these reasons we chose the Allen et al. (2012) IPEs as the most suitable for our study. These IPEs were generated using 13,077 earthquakes data around the world since 1973; 1613 of the earthquakes are from Asia. The IPEs (Eq. 1) were used to estimate the intensity of the earthquake at each distance to the rupture zone and Eq.
(2) was proposed by Borcherdt (1997) and Cinicioglu et al. (2007) to compute the site amplification effect and it's utility as applied to shallow crustal earthquakes in Banda Aceh city: Where I stands for the intensity (MMI), M is the earthquake magnitude (Mw), R rup is the closest distance to the rupture zone (km), c 0 , c 1 , c 2 and c 3 are coefficients with the values 3.950, 0.913, − 1.107, and 0.813, respectively, and S is the site correction where F represents the amplification capacity of the local soil (Fv for long periods, Fa for short periods).

Site coefficient and liquefaction effect
The site coefficient for Banda Aceh city was determined using shear-wave velocity travel-time averaged over a depth of 30-m (Vs30) data from Muzli et al. (2014) and Jamaluddin et al. (2018) in order to calculate the amplification factor. Muzli et al. (2014) and Jamaluddin et al. (2018) determined Vs30 using the MASW (Multichannel Analysis of Surface Wave) method from 53 locations, which were placed approximately 1 km apart. The result showed that the Banda Aceh area has stiff and soft soil layers. To determine the amplification factor from the Vs30 data, this study used Eq. (3) from Borcherdt and Eeri (1994) to consider how a mid-short period amplification of earthquake ground motion could affect Banda Aceh: Where Fv is the mid-short period of the amplification factor, 1050 is the shear wave velocity in firm to hard rock and mw is the value for the peak ground acceleration (PGA). According to the Tim Revisi Peta Gempa Indonesia (2010), the PGA value for 10% exceedance in 50 years in Banda Aceh is 0.4 g, and is dominated by the GSF. According to a prior study of Borcherdt and Eeri (1994), the mw value at 0.4 g is 0.45.
Liquefaction is one of the side effects of the seismic ground motion that can damage buildings and infrastructure. Several factors influence liquefaction following an earthquake, including magnitude, intensity and the water table level (Holzer et al. 2011). The water table levels of the groundwater in Banda Aceh city were determined using a 1D Geo-electrical method in this research. This method was applied to 20 locations and combined with the liquefaction vulnerability map of Banda Aceh City that was developed by Tohari et al. (2015) as shown in Fig. 3. The water table map and liquefaction vulnerability map were then utilised to produce the new liquefaction map. Groundwater depth is a critical parameter in studying the liquefaction potential zone (Chung and Rogers 2013), and the effect of varying the depth of groundwater influenced liquefaction by 20% (Tonkin and Taylor Ltd 2013). The detailed weights used in the GIS analysis in this study are shown in the flowchart (see Fig. 2).

Earthquake loss prediction analysis
Banda Aceh city lies 2.5 km far from the Aceh segment and 15 km from the Seulimeum segment. Earthquakes occurring on these faults are only a short distance from the city, and given that these are shallow crustal earthquakes large events are likely to result in high levels of building damage and casualties among residents (Allen et al. 2012). Different building typologies have different resistances to ground shaking, which also varies from one area to another and can be estimated from the fragility curves developed by Tingatinga et al. (2013) and Naguit et al. (2017). The types of building studied by them were similar to the typical buildings proposed by FEMA (2010)  to 1 as a function of ground shaking/intensity in MMI, with the value 1 corresponding to 100% damage. The data of 79,421 buildings in Banda Aceh city were inventoried and divided into three typologies; (1) Low-Rise Concrete Moment Frame (C1-L), 78,006 buildings; (2) Low-Rise Light Wood Frame (W1-L), 802 buildings; (3) Medium Rise Concrete Moment Frame (C1-M), 612 buildings. The distribution of the types of buildings in Banda Aceh city is shown in Fig. 4a. To calculate the damage ratio for the different construction types, the fragility curves in Fig. 4b were used, which express damage ratio in terms of the intensity experienced by each building. C1-L and C1-M buildings have nearly similar trend of fragility curves, but the W1-L curve has lower damage ratios. The lighter wood material makes the W1-L building more resistant to the earthquake ground motion compared to C1-L and C1-M, as has been shown by field study conducted by Naguit et al. (2017).
The damage ratio associated with each building yields varying numbers of casualties among people trapped in a collapsed or heavily damaged building (Hashemi and Alesheikh 2011;Rusydy et al. 2017a;Rusydy et al. 2018a). Building collapse, the primary cause of human injury in earthquakes, results in a fatality rate of 75% (Coburn and Spence 2002). The casualty rate for each building is more difficult to estimate as it differs from one earthquake to another, and the existing data on casualty rates yields an inadequately wide range of numbers of people injured. The number of injuries depends on each building's population. It is difficult to record the population for each building to solve this equation; therefore, this study calculated the predicted population data of Banda Aceh city and the total area of all houses in Banda Aceh city. This calculation is applicable only to residential buildings. The equations of house population are shown in Eqs. (4) and (5). This paper used the Eq. (5) that developed by Hashemi and Alesheikh (2011) to estimate the number of injuries. The number of injured in that equation was based on the occupancy at the time of the earthquake, damage ratio, population of each building, and injury ratio. The worst case scenario for the injury ratio is as shown in Table 1. Table 2 shows the population distribution in all sub-districts in Banda Aceh based on the 2017 population census.
Injury ¼ Occupancy ratio x population x damage ratio x injury ratio Based on building damage ratio models, this study was able to estimate the economic loss for each buildings as well as each sub-district. The economic loss estimation model are an important information for decision makers and policy makers (Guettiche et al. 2017). Financial loss estimation study due to earthquake event, has been introduce by Sun et al. (2012); Jaiswal and Wald (2013); Guettiche et al. (2017). This study, applying the Eq. (6) to estimate the economic loss each sub-districts in Banda Aceh city.
Economic Loss ¼ damage ratio x building area x stories x price of building The damage ratio of a building was derived from the model of damage for different earthquake scenarios on the Aceh and Seulimeum segments. The building information that have been collected previously by GIS analysis and building inventory, include area and number of Fig. 4 a The building typology distribution in Banda Aceh city measured in 2016. b Heuristic fragility curve of 3 variants of buildings in Banda Aceh (Tingatinga et al. 2013;Naguit et al. 2017) Table 1 Injury ratio for worst-case damage scenario in different type of building according to Hashemi and Alesheikh (2011) No.
Type of buildings Injury ratio   (Indonesia Rupiah) and USD (United States Dollar) at Banda Aceh was divided among 5 different types of buildings and depends on the different type of occupancy and type of construction provided in Table 3. Public buildings, including schools and governmental offices, are typical of building type B, while special public buildings such as hospitals, laboratories, and hotels were referred as building type A. One and three storeys reinforced concrete houses (C1-L) were referred as houses type B, while timber construction houses were referred to as houses of type C.

Validation of the models
The accuracy of the earthquake intensity and damage prediction models is of concern. As a first step, the models of IPEs (Intensity Prediction Equations) and building damage for each building typology should be validated. Such a validation was performed for previous recorded earthquakes in the region. Earthquake models are very sensitive to magnitude and shaking attenuation. The damage models rely on the fragility curve, and the injury model depends on the injury ratio. Indonesia has numerous large magnitude shallow crustal earthquakes. The IPEs equation proposed by Allen et al. (2012) was used in this study to develop the earthquake models. The IPEs was validated by a review of recorded historical shallow crustal earthquakes in Aceh province. The most memorable and well documented shallow, crustal earthquake in Aceh province was the Tanoh Gayo earthquake, which occurred on July 2, 2013 with magnitude Mw 6.1, at a depth of 10 km (Rusydy et al. 2015). The Tanoh Gayo earthquake was determined to be a shallow crustal earthquake. Rusydy et al. (2015) and Daryono and Tohari (2016) noted that this earthquake source had not been recognized as part of the GSF. The earthquake source belonged to another segment of faults that had not been recognized in the earthquake map previously. Using macro seismic surveys and surface ruptures, Rusydy et al. (2015) and Daryono and Tohari (2016) found the segment is part of a previously unknown fault that crosses Takengon City striking to the north-west. Based on their findings, this study used the rupture zone as the source of the Tanoh Gayo earthquake as shown in Fig. 5. The Vs30 used to estimate the amplification factor for the Aceh Tengah and Bener Meriah districts, which were affected by the Tanoh Gayo earthquake, was determined by Rusydy et al. (2017b) using multichannel analysis of surface waves (MASW) measurements. The actual intensity of the 2013 Tanoh Gayo earthquake was derived from Rusydy et al. (2015) in a survey performed a few days after the Tanoh Gayo earthquake.
The intensity validations showed some differences between the IPEs and the actual intensity (see Table 4). For the Timang Gajah area, the actual intensity was in the range of 7-9 but the IPEs appear model the maximum intensity in this range. For other areas, the IPEs values are shown to be higher than the actual intensity. Overall, the IPEs intensity are overestimates by about 0.1-0.5 MMI. These overestimates could be due to the IPEs effectively assuming uniform slip over the rupture area. In fact, each rupture involves along-strike variations in slip which results in different spatial distributions of intensity. On other hand, the local geology can also cause dissimilar amplification values; accordingly, it varies from one site to others and it depends on local soil classifications as studied previously by Karimzadeh et al. (2017).
The damage model validation was conducted using a visual assessment method. A few days after the 2013 Tanoh Gayo earthquake, photographs of damaged buildings at selected validation points were taken. Those photos were used to determine the damage ratio. Hence, the visual classification of damage ratio referred to the work of Okada and Takai (2000). The result of the validation process is shown in Table 5.
In Timang Gajah subdistrict, the Low-Rise Concrete Moment Frame (C1-L) buildings experienced moderate (0.2) to heavy damage (0.6), as did those in Bebesan and Silih Nara due to the Tanoh Gayo earthquake. In the Ketol subdistrict, C1-L buildings experienced heavy damage (0.4-0.6). Accordingly, the damage model for C1-L buildings produced estimates for damage in the Tanoh Gayo earthquake that are in the range of the actual damage. Low-Rise Light Wood-Frame (W1-L) buildings suffered moderate (0.1-0.2) to heavy damage (0.4-0.6) in Timang Gajah and Ketol but none to moderate actual damage in Silih Nara. The actual damage of W1-L buildings was similar to the model in Timang Gajah, but overestimated in Ketol.

Earthquake intensity models
The earthquake intensity model was based on an IPEs designed for global applications (see Eq. (1)). This equation was used to predict the intensity of the earthquake in the bedrock and at the surface. The effect of site amplification on the surface intensity are calculated using Eq.
(2) with appropriate site coefficient. At large scale (see Fig. 6), the Vs30 data for calculating the site coefficient using Eq. 3 was taken from the global Vs30 model developed by Allen and Wald (2009), yet for small scale (see Fig. 7a), the Vs30 data were derived from Muzli et al. (2014) and Jamaluddin et al. (2018). These data then were applied in GIS analysis to produce the earthquake scenario ground motions shown in Fig. 6. Earthquake scenarios correspond to Mw 7 and Mw 6.5 earthquakes for each of the Aceh and Seulimeum fault segments. A previous study noted that the highest potential magnitude earthquake that could be produced in the Aceh segment is in the range Mw 6.5-7.0 (Ito et al. 2012).
The 1964 earthquake at the Seulimeum segment had a magnitude of 6.5 (Sieh and Natawidjaja 2000) and the maximum magnitude event could be Mw 7. The models showed that the degree of intensity at a particular point depends on the distance from the earthquake and the local site coefficient. The USGS value of Vs30 used to calculate the site coefficient is very suitable for regional earthquake models, but is less applicable to local earthquake models. To determine the details of spatial variation of the earthquake scenario intensities in Banda Aceh itself, this study used values of Vs30 taken form the MASW studies of Muzli et al. (2014) and Jamaluddin et al. (2018).

Site coefficient analysis
Detailed Vs30 studies in Banda Aceh city were conducted by Muzli et al. (2014) and Jamaluddin et al. (2018), and these data were used here to calculate the site coefficient and amplification factors. The locations of the Vs30 data and the interpolated results are shown in Fig. 7a. Based on these results, this study calculated the amplification factor for Banda Aceh. Those results are illustrated in Fig. 7b. The amplification factor (Fv) Fig. 5 Rupture zone for the Tanoh Gayo earthquake 2013 and intensity map produced using the IPEs and the location of validation points. The rupture zone is set according to Rusydy et al. (2015) and Daryono and Tohari (2016). Digital elevation model derived from NASA 2007 was substituted into Eq.
(2) to compute the site coefficient (S). The results of the site coefficient data used in developing the earthquake intensity models associated with the nearest fault system to Banda Aceh are shown in Fig. 7c, d, e, f. Soft soils (low Vs30) trigger a higher amplification factor, and stiff soils (high Vs30) yield a lower amplification factor (Rusydy et al. 2017b).
The highest amplification factors for the central and western parts of Banda Aceh city ranged from 2.3 to 2.7. The lower amplification factors were found for the northern and southern areas of the city. The estimated ground motion intensity produced in Banda Aceh due to the Aceh segment scenario earthquakes are 8.4-9.3 MMI for magnitude Mw 7 and 8.2-9.2 MMI for magnitude Mw 6.5. The high intensity located in the southwest of Banda Aceh city is near to the source of the earthquake, as shown in Fig. 7c, d. The earthquake intensities estimated for the Seulimeum segment scenario earthquake of magnitude Mw 7 are in the range of 7.6 to 8.6 MMI, with the high intensities occurring in eastern Banda Aceh, as shown in Fig. 7e. Table 5 The damage comparison between our damage model and actual damage from the 2013 Tanoh Gayo earthquake   Sub-district  Intensity

Liquefaction effect
The groundwater table level is one parameter that triggers the liquefaction (Chung and Rogers 2013). It influences the liquefaction effect by about +/− 20%, compared to other parameters (Tonkin and Taylor Ltd 2013). Figure 8a shows the depth groundwater table in Banda Aceh city. These data are based on 20 locations of 1D Geo-electrical measurement. The map shows the interpolated data Fig. 7 Site coefficient map and intensity in Banda Aceh city. a Vs30 map around Banda Aceh city modified from Muzli et al. (2014) and Jamaluddin et al. (2018). b Amplification factor map and c Intensity map of Mw 7 magnitude earthquake from Aceh segment. d Mw 6.5 from the Aceh segment, e Mw 7 the Seulimeum segment, and f Mw 6.5 the Seulimeum segment among the sites. The calculated liquefaction map of Banda Aceh based on the groundwater table and vulnerability liquefaction map is shown in Fig. 8b. This study considered additional building damage caused by liquefaction, depending on the level of liquefaction vulnerability indicated in Fig. 8b. Liquefaction is assumed to result in 2.5, 5, 7.5 and 10% additional damage to buildings in light, medium, high, and very high vulnerability zones, respectively.

Damage distribution in several earthquake models
In this study we generated four earthquake scenarios, two produced by the Aceh segment (Mw 7 and Mw 6.5) and two by the Seulimeum segment. Using GIS analysis, the earthquake intensity maps (Fig. 7c, d, e, and f) are clipped by the building polygon data (Fig. 4a). Building damage associated with various types of buildings was calculated using the fragility curve shown in Fig. 4b. This study calculated the damage ratios for each type of building in Banda Aceh city due to earthquake intensity and up to 10% additional damage cause by liquefaction. The results of the distribution of the building damage for each earthquake scenario are shown in Fig. 9.

Human injury estimation
As part of the injury estimation analysis, population size, damage ratio, and intensity had to be precisely determined. This study attempted to build the intensity map, the building population sizes, and the damage ratio as precisely as possible. The injury estimations in different earthquake scenarios were examined at two times during the day (2 PM and 2 AM), even though FEMA (2010) has suggested the use of three times per day (2 PM, 5 PM and 2 AM). This study considered two times per day because we assumed that people have similar patterns of activity at 2 PM and 5 PM.

Injury prediction during the day-time (2 PM)
The time of occurrence at 2 PM is a busy time when many people in Banda Aceh city are away from their homes and working in their offices. In a daytime event, the occupancy ratio used was 1 for offices, 0.3 for houses, and 0.5 for dormitories and hotels. Prediction of the number of injuries to people was calculated using Eq. (5) and the percentage of the population. Figure 10 denotes the distribution of injury numbers around Banda Aceh city for several earthquake scenarios in percentages. This scenario considers that the people of Banda Aceh have no quick response plan or adequate access to evacuation. The scenario I in Fig. 10a is the worst-case scenario associated with an earthquake produced by the Aceh segment for Banda Aceh city. This scenario will injure 3.5-20% of the population buildings in almost all sub-districts. Scenario II also will harm 3.5-20% of the population in buildings. Scenario III is the worst-case scenario for an earthquake event from the Seulimeum segment of magnitude Mw 7. This earthquake will injure 2.3-15% of population in buildings. The IV scenario, with Mw 6.5 will injure 2-15% of the population in buildings.

Injury prediction at night-time (2 AM)
In a night-time earthquake, all the residents of Banda Aceh city are assumed to be in their houses, and the offices have zero occupation. In this scenario, this study could compute the number of injuries based on the population data of residences as shown in Table 2. The predicted population was in turn calculated with GIS analysis using population data from each village, the total area of the houses in the village and the area of each building, as determined in Eq. (4). The occupancy ratio at this time is 1.0 for all homes and residences, and 0 for offices, schools, and shops. The potential number of injured people in Banda Aceh city during a night-time event is shown in Fig. 11.
In a worst-case scenario for an earthquake on the Aceh segment, the number of injured will range from 1 to 2 people in each building. There was little difference in the number injured among scenarios I, II, III, and IV. In the total, approximately 20,264 people, or 7.8% of the population of 261,034 of Banda Aceh, would be injured in scenario I. The most casualties would be in the subdistricts of Kuta Alam (3032 casualties), Baiturrahman (2455 casualties), and Syiah Kuala (2427 casualties).
Scenario III is the worst-case earthquake from the Seulimeum segment. This segment is located in the eastern part of Banda Aceh, opposite to the Aceh segment. This event will injure 1-4 people in each building. The total number injured in this event would be approximately 14, 857 or 5.7% of Banda Aceh residents. The greatest numbers of casualties are in the sub-districts of Kuta Alam (2282 injuries), Syiah Kuala (2081 injuries), and Baiturrahman (1819 injuries). The detailed injury numbers for each district is shown in Table 6.

Economic loss estimation
The analysis of the damage distribution from the highest magnitude earthquake scenarios associated with the Aceh and Seulimeum segments are provided in Fig. 9. Figure 12 shows that the most affected districts in term of economic loss are sequentially Kuta Alam, Syiah Kuala, and Baiturrahman for the Mw 7 Aceh segment and Seulimeum segment scenarios. Generally, the Aceh segment scenario is predicted to have higher damage and associated cost than the Seulimeum segment. This is because Kuta alam district, the most populated district in Banda Aceh city with many government offices, is closer to the Aceh compared to the Seulimeum segment. In total, economic loss in Banda Aceh would be 46,554, 706 million IDR or 3320 million USD if the Mw 7 earthquake occur on the Aceh segment and 34,710,843 million IDR or 2475 million USD if it occurred on the Seulimeum segment. More detail of economic loss estimation for each district can be seen in Fig. 12.

Conclusions
This study successfully developed several earthquake scenarios and determined the worst case scenario affecting the residents of Banda Aceh during daytime and nighttime earthquake events. The earthquake scenario ground motions in MMI were developed using globally applicable, macroseismic IPEs designed by Allen et al. (2012). Because these IPEs were developed for shallow crustal earthquakes, they suit the context of Banda Aceh situated along the Great Sumatran Fault (GSF). Two segments, i.e. the Aceh and the Seulimeum segments which are part of the GSF, are characterized as shallow crustal earthquake sources that have potential to generate severe impacts once major earthquake occurs. This study estimated the human casualties potentially caused in Banda Aceh by earthquakes generated from these segments. The intensity validation model used with respect to Aceh province, along with the effects of local geology, work well to correctly estimate an intensity model. The fragility curve of C1-L developed by Tingatinga et al. (2013) and Naguit et al. (2017) are compatible with parameters found in Aceh province.
Various types of buildings and ground shaking will give rise to different damage ratios. In the worst scenario from the Aceh segment, 20,264 residents of Banda Aceh will be injured, which are about 7.8% of the total population of the city. The worst-case scenario associated with the Seulimeum segment would injure 14,857 people. Both of the scenarios were assumed to occur at 2 AM, Fig. 10 Injury percentages distribution for the daytime scenario of an earthquake. a Intensity map of Mw 7 from Aceh segment, b Mw 6.5 from Aceh segment, c Mw 7 from Seulimeum segment, and d Mw 6.5 from the Seulimeum segment when most people are at home. As explained earlier, the quality control on private and non-government building construction in Indonesia is weak, this could also contribute to the large number of potential human casualties due to the earthquakes. The emergency situation potentially generated by such an earthquake could put disaster emergency response agencies under huge pressure. We suggest the Disaster Management Agency of the city to increase their capacity in order to give punctual and effective measures to respond to such emergency situations. Implementation of building codes and quality control on private and non-government buildings should be imposed.  Fig. 12 The infrastructures economic loss each districts in Banda Aceh city due to Mw 7 earthquake from Aceh segment and Seulimeum segment