Correlation between Building Damages and Losses with the Microzonation Map of Mataram—Case Study: Lombok Earthquake 2018, Indonesia

: The high intensity of the earthquake on Lombok Island on 5 August 2018, with a magnitude of 7.0 Mb, caused material losses experienced by the affected residential areas. The Indonesian Geological Agency in 2015 published a microzonation map that mapped zones prone to earthquake shocks to mitigate disasters. This study aimed to compare the level of damage and loss in residential areas due to earthquakes in Mataram City with earthquake-prone zones using a microzonation map. The correlation between damage and loss value of residentials with microzonation maps was evaluated using the overlay method. The results showed that the level of damage and the value of the loss of houses in the high disaster-prone zone (red zone) showed the highest loss value. In comparison, the level of losses in the moderate disaster-prone zone (yellow zone) and light disaster-prone zone (blue zone) on the microzonation map shows a low and lower loss value. This study concludes that the microzonation map helps determine the damage zone and the level of disaster vulnerability caused by the earthquake hazard.


Introduction
The location of Indonesia in the collision zone causes frequent natural disasters. Earthquakes are one of such disasters caused by tectonic activities. The dynamic interaction of numerous tectonic plates in eastern Indonesia causes high seismicity rates, which resulted in catastrophic damaging earthquake sequences on Lombok Island in 2018 [1][2][3]. Based on the National Disaster Management Agency (BNPB) records, from the beginning of 2017 to June 2018, there have been 4006 disasters in Indonesia, with 30 tectonic earthquakes on Lombok Island occurring in August 2018 with a magnitude of 5.9-6.9. The earthquake in Lombok claimed at least 481 lives and caused an economic loss of US$ 514,000,000 as of 17 August 2018 [4]. These financial losses are caused mainly by damage to infrastructure and buildings.
Natural disasters are a significant concern in Mataram City. A natural disaster is an unanticipated event or set of circumstances that endanger and disrupts people's lives and livelihoods as a result of natural causes, both by natural factors and/or non-natural factors and human factors, resulting in fatalities, environmental damages, property losses, and psychological impacts [5]. Examples of natural catastrophes include floods, tsunamis, earthquakes, and landslides [6,7]. After the earthquake disaster in Lombok in August 2018, the damage level of the Municipality of Mataram ranked fourth out of six regencies in Lombok with a total loss of US$ 69,008,289. The tremendous damage and loss to the settlement sector reached a value of US$ 34,907,209, damage to the housing subsector reached US$ 32,500,610, and other losses, including cleaning of debris, worth US$ 1,600,470 [4].
Earthquakes are the vibrations or shocks that occur on the surface of the Earth caused by the collisions between the tectonic plates, active fault zones, volcanic eruption activity, and rock debris. Lombok Island is located north of the subduction zone, where the Indo/Australian Plate collides and subducts below the Eurasian Plate with a convergence rate of 44-68 mm/yr [8].
It is crucial to distinguish between risk and vulnerability. Seismic risk is the risk of damage and consequential loss to a particular structure or collection of structures over a specific period. The anticipated losses from all degrees of hazard severity are included in risk, taking their occurrence probability into account, whereas the vulnerability of an element is generally specified for a specific hazard severity level. Seismic vulnerability analysis in urban environments focuses mainly on building structural, geological, and geographical information directly related to potential human and economic losses in the case of a seismic episode [9]. On the other side, Lombok Island is locked between two active faults, the Flores back arc thrust on the north with a convergence rate of ca. 9.9 mm/yr and the normal fault of Sumbawa with a rate of 0.5 mm/yr [10]. Those faults influenced the earthquake events in Lombok Island ( Figure 1).
The earthquake with the greatest strength occurred on 19 August 2018, with a magnitude of 6.9 Mb, and the epicenter at a radius of 30 km NE of East Lombok. At the same time, five other earthquakes occurred with a magnitude of 4.9 Mb-5.7 Mb, whose shaking intensities were scattered in several areas in the island of Lombok, with the distribution shown in Figure 2. The earthquake on the island of Lombok, taking place in July-August 2018, had negative impacts, such as fatalities, damages, and material losses. From January through October 2018, there were 98 earthquakes with a magnitude ranging from 4.1 Mb-6.9 Mb. The depth of the earthquake sources was between 5.7 km-121.9 km [11]; the earthquakes that occurred were dominated in the NW-NE region of the island of Lombok.
The microzonation map of Mataram City from the Geological Agency (2015) is used as an estimated vulnerability zoning model. The map is used as a basis to prepare a zonation map by using the microtremor. Bertelli and Omori were the first to discover microtremors in 1909, while Guttenberg was the first to introduce engineering uses of microtremors [12,13]. Microtremor measurements are valuable for assessing seismic threats in the near-surface substructure [14,15]. Microtremor (or ambient vibration) methods analyze the mechanical qualities of the Earth's subsurface, particularly seismic velocities, by measuring background seismic noise. Seismic noise is the constant vibration of the Earth's surface caused by a combination of the low-frequency (less than 1 Hz) natural phenomena (earthquakes, wind, tides, rivers, rain, variations in atmospheric pressure) and high-frequency (greater than 1 Hz) human activities [16,17]. Idriss and Boulanger [18] showed that damage to building structures due to earthquake intensity of ground shocks during an earthquake was significantly influenced by local geological and soil conditions. The microzonation divides or classifies areas of zones that have relatively similar potential shocks and disasters. Their engineering characteristics have local effects. Hard rock sites will experience a low shock intensity, while soft soil will experience a high shock intensity [19][20][21][22].
The purpose of this research was to determine the accuracy of the microzonation map against actual building losses and damage when the Lombok earthquake in 2018 occurred. The high accuracy of microzonation maps is perfect for development plans such as urban planning and financial risk transfer strategies. The burden of the state budget will be reduced in the future if disaster impact mitigation using the prediction of microzonation maps goes well; in turn, it will also reduce casualties and material losses.

Materials and Methods
The availability of microzonation maps is one of the factors that makes Mataram City a research area, with well-collected data on losses and damages to buildings after the 2018 Lombok earthquake. The city of Mataram has six districts, such as the Ampenan District, Sekarbela District, Selaparang District, Mataram District, Cakranegara District, and Sandubaya District. All districts have different values of loss and deterioration to the buildings and their microzonation areas.

Study Area
Mataram is a land relatively flat, undulating topography, with a slope of 0-2% covering an area of 3.216 hectares and a bumpy area with a slope of more than 2-14% surrounding an area of 2.909 hectares. The altitude of Mataram City is approximately 0-50 m above sea level. These conditions indicate that most of the city is a plain area. The western part of the city is occupied by the flat-sloping area while the eastern has rather high-bumpy relief. The overlay area, with flat physiography, has two conflicting implications. First, flat areas have positive values, namely the construction of infrastructure and facilities which physically have fewer technical obstacles, and the development funding is relatively cheaper. Second, a flat area with an altitude almost parallel to the sea level has terrible implications, including flood-prone areas.
Based on the Geological map of Lombok sheet West Nusa Tenggara [23], three rock formations are covering the research area: Kalikupang Formation (Tqp), Kalibabak Formation (Tqb), and Lekopiko Formation (Qvl) that consists of volcanic rocks, sedimentary rocks, and intrusive rocks whose ages range from Tertiary to Quaternary. The engineering geological map of the Lombok Island, Mataram, lays on the geological unit of the Alluvium deposit (Al), consisting of the river, coastal, and swamp deposits, composed of silty sand, sandy silt, and loose sand [24]. Sedimentary swamp comprises sandy silt, silty clay, black to grey color, fine-medium grained, very soft-firm, medium to high porosity. River sediment consists of sandy silt to silty clay and loose sand, brownish yellow-brown. It contains gravel, fine-medium-grained, very soft-dense, low to moderate plasticity, 3.50-6.50 m thick. Coastal deposits are black to grey sand, fine-grained, high porosity, containing shells. The engineering geological constraints or geological disasters prone are floods, abrasions, and tsunamis.

Damage and Loss Due to the Mataram Earthquake in August 2018
The earthquakes hit Mataram City, impacting six districts, including Ampenan, Sekarbela, Selaparang, Mataram, Cakranegara, and Sandubaya. The first earthquake occurred on 5 August 2018, with a magnitude of 7 Mb at 19 [4].
The impacts of the Lombok earthquakes caused damages and losses in various sectors such as residential, infrastructure, social, economic, and cross-sectoral, which have disrupted public activities and services in the affected areas based on the data from the Post-Natural Disaster Reconstruction and Rehabilitation Action Plan. This natural disaster caused total damage and loss of 69,008,289 US$. The tremendous damage and loss occurred in some sectors comprising the residential sector amounted to 34,905,750 US$; the social sector amounted to 25,335,416 US$; the cross-sector amounted to 7,023,475 US$; the productive economic industries amounted to 1,743,648 US$.
The amount or pattern of building damage caused by the earthquake that impacted Mataram demonstrates the pattern's conformance with the microzonation map. The proportion of building damage by the district is depicted in Figure 3. The Sekarbela district (21.66 percent) and Ampenan (19.76 percent) have suffered the worst damage, while the Mataram district (8.75 percent), Selaparang (8.75 percent), and Sandubaya (7.76 percent) have sustained minor damage. The Cakranegara district in the blue zone is the least affected (1.7 percent).

Microzonation Map
The Potential Earthquake Ground Shaking Hazard Map (Microzonation) Mataram City, West Nusa Tenggara Province, produces three types of zones, namely red, yellow and blue zones (Figure 4). The red zone indicates that the area has a relatively high level of shaking with the morphology of volcanic fluvial plains. The yellow zone means that the area has a fairly medium level of shock in volcanic fluvial plains. The blue one indicates that the area has a relatively low level of shock with the morphology of the undulating volcanic plain of volcanoes. Regional zones on the Microzonation Map of Mataram involve four seismic microzonation parameters, which are dominant periods, soft sediment thickness, average rock/soil amplification, and soil rock site classification. The red zone has the characteristics of a dominant period of more than 0.5 s, a soft sediment thickness of more than 30 m, an average amplification ranging from 1.9 to 2.3, and a site class of soft to medium soil. While the yellow zone has the characteristics of a dominant period ranging from 0.25 to 0.5 s, a thickness of soft sediment ranging from 10 to 30 m, an average amplification ranging from 1.6 to 1.9, and a medium soil site class. Finally, the blue zone has the characteristics of a dominant period of less than 0.25 s, a thickness of soft sediment below 10 m and amplification below 1.6, and a hard soil site class.
After generating three types of zones based on the level of shock on the Earthquake Potential Map (Microzonation) of Mataram City, West Nusa Tenggara Province, these zones can be developed based on disaster mitigation recommendations for safer development.
The red zone has a relatively high level of shock. Therefore, this area is not recommended to construct houses or critical and emergency (IV) facilities. However, once facility (IV) has been built, it needs to be reviewed with special geological studies, site investigations, and specialized foundation designs. The yellow zone with a relatively moderate level of shock is recommended for medium-risk (II) facilities built in this zone. The construction of medium-risk (II) facilities must comply with the provisions of the local building code. The blue zone, which has a relatively low level of shock, is recommended to build facilities in building categories I, II, III, and IV (Table 1). According to the microzonation map, it has been found that each district has a zonation of the potential hazardous ground shocks that differ according to the distribution of soil and rocks characteristics. The largest red zone coverages are respectively in Sekarbela (SK), Ampenan (AM), Mataram MA), Selaparang (SL), Sandubaya (SD) and Cakranegara (CA). The results of microzonation research are directly proportional to the damage and losses that occur. These facts are summarized in Table 1, where it has been found that the Sekarbela (SK) district, with the most significant red zone of 88.52%, experienced enormous damage and loss. Meanwhile, the Cakranegara (CA) district, which has the largest blue zone of 85.37%, experienced the slightest damage and loss.

Worksteps
The Geographic information system is a computer-based information system that focuses on the geography of an area, which can be used and designed to compile, manipulate, process, display, and analyze data that have spatial information. The relation between equal intervals, natural breaks, manual and statistical examination, all of which are GIS-default operations, are the most common mathematical approaches for data classification [25].
The GIS-based data for the land-use suitability analysis has been used in various situations, including geological favorability [26][27][28]. Simple processes, such as laying a road map over a map of local wetlands, or more complex processes, such as multiplying and adding map characteristics of different values to calculate averages and correlations, can be represented by these outputs. This process could be a visual action at its most basic level, but analytical activities require combining one or more data layers [29][30][31][32].
The superimposed research was conducted by correlating the Earthquake Hazard Potential (Microzonation) map of Mataram City, West Nusa Tenggara Province, with data on damage and losses in each district in Mataram Municipality, based on the BNPB report [4]. The comparison of the two data aimed to determine the correlation between both data is either directly or inversely proportional, or even irregular way. The comparisons were also made on the residential area data from each district obtained from the digitization results of Google Earth satellite imagery. It was performed because the variable density of residential could affect the value of losses caused by earthquake disasters.
This research uses the weighted overlay GIS method ( Figure 5). The weighted overlay method analyses spatial data using the overlay technique of several raster maps related to the factors that affect the vulnerability assessment of a problem. Weighted Overlay can combine various inputs in a grid map [33]. This method solves problems with many criteria, such as optimal site selection or suitability modeling. Applying the latter approach to the case of the loss and damage data values and the microzonation map units above (Figure 6), the average element values can be obtained using Equation (1) [26]: where are: i = microzonation units j = loss damage data values z = variable in the source zones ∂ij = overlapping target zone

Results
The level of damage to buildings after the Lombok earthquake (Table 2) The value of losses due to house damage is calculated based on the standard value of the level of damage due to the Lombok earthquake (Table 3), where total loss (100%) is US$ 6320; heavily damaged (80%) worth US$ 5056; moderate damage (50%) worth US$ 3160; and Low (20%) at $1264. Then the value of the loss for each district (Table 4)    The loss rating and dominance of the red zone reflect a pattern of conformity based on the magnitude of losses caused by the earthquake ( Table 5). The ranking of earthquake losses is directly related to the prevalence of the red zone distribution in the districts of Sekarbela (SK) and Ampenan (AM). It implied that areas with red zone domination post-earthquake would also suffer significant losses. Based on the rating of the loss value compared to the level of the vulnerability zone (red), it shows the suitability of the pattern. The value of losses due to the earthquake is directly proportional to the dominance of the red zone distribution in the Sekarbela and Ampenan districts. These show that postearthquake, subdistricts with high red zone dominance will also have high losses.

Discussion
The research area has a thickness of soft rock (volcanic fluvial deposits) which causes earthquake waves to be trapped and amplified, explaining the higher degree of damage. The distribution of damage from the earthquake on Lombok Island was irregular. The study shows that the microzonation map helps determine the damage zone and level of disaster vulnerability due to the earthquake hazards. Earthquake hazard microzonation mapping is an essential tool for land use planning in infrastructure development and mitigation strategies. It provides easy-to-read charts and maps, making it easier for governing bodies to make decisions. It is highly recommended that microzonation studies support urban planning development, residential areas, and areas of people's economic activities in disaster-safe areas.
Microzonation studies over large areas can support urban development plans and designs. This study can mitigate potentially earthquake-prone areas by recording the physical properties of rocks in static conditions so that determining the designation of an area can be more manageable. This research is relatively inexpensive, easy, and suitable for developing and utilizing regions. The blue zone area must be designated for housing, buildings, and areas for community economic activities. While the yellow zone is suitable for developing community activities such as workshops and carpentry, locations included in the red zone allocates as a buffer area for environmental conservation that is reforested by governance programs to reduce carbon emissions.
The level of damage to areas that caused significant losses is dominant in the red zone, although the highest population is not in the areas where the loss rate is high. Local governments should reform land usage and distribution following the disastrous Lombok earthquake.
Based on the comparison between the level of building damage and the value of the loss, it concluded that although the number of damaged houses in the Ampenan district is the highest (4773 units) worth US$ 10,327,512, the loss value is still below that of the Sekarbela district (3774 units), but the loss value is US$ 12,399,840. The number of houses with severe damage dominates the number of damaged in Sekarbela (1540 units), while the value of losses in other districts is proportional to the level of damage. The high level of loss in Sekarbela, caused by the dominance of the number of heavily damaged houses being more than Ampenan.
Future development planning may be more successful in site selection for important infrastructure investment decisions with a scientific understanding of each earthquake hazard level [34,35]. During the 2018-2020 timeframe, 61.000 homes were built for earthquake victims, ranging from mildly damaged to badly damaged homes, while 225.000 buildings were constructed for public amenities and social functions, and millions of dollars were spent to recover the damage [36]. Risk transfer is a strategy to cope with the effects of natural disasters. Natural disasters are unforeseeable, unpredictable phenomena that can have a devastating impact on the population, wreaking havoc on infrastructure and causing enormous human and economic losses [5]. An alternative to risk transfer is disaster insurance, especially earthquake insurance. Earthquake insurance is essential to recover from earthquake disasters and protect from the associated financial losses [37,38].

Conclusions
Its location in the collision zone causes frequent natural disasters in Indonesia. Earthquakes are one of such disasters caused by tectonic activities. The dynamic interaction of numerous tectonic plates in eastern Indonesia causes high seismicity rates, which resulted in catastrophic damaging earthquake sequences on Lombok Island in 2018. The research area in the present study has a thickness of the soft rock (volcanic fluvial deposits), causing earthquake waves to be trapped and amplified, which explains the higher level of damage. The distribution of damage caused by the earthquake on the island of Lombok is irregular. Therefore, it can be concluded that mapping hazard zones based on microzonation are directly proportional to the damage and losses caused by the earthquake. In other words, areas included in the red zone are the most vulnerable to damage and losses. In other words, areas included in the red zone are the most vulnerable to damage and losses.
Microzonation research is a valuable tool for mapping potential earthquake threats in a particular area, making research more detailed and accurate. In addition, it is essential to conduct regional tectonic research and its influence on local geological structures that impact Tertiary bedrock covered by Quaternary volcanic rock from the eruption of Mount Rinjani. It is suggested that local government improve the "building code" rules when applying for building permits.
Based on the findings of this study, it is also suggested that additional models be investigated to predict the amount of seismic activity, such as the Markov Chain model, which has been applied to the Algerian region [39]. As a result, the seismic activity can be observed using multiple models to gain a more comprehensive understanding. Meanwhile, since this region is located in a steep area, land use for forests and farmland does not require considerable consideration due to the earthquakes. As a result, substantial consideration must be given to slope stability, building structures, and available types of settlement materials [40].