Geospatial Analysis of Earthquake Damage Probability of Water Pipelines Due to Multi-Hazard Failure

The main purpose of this study is to develop a Geospatial Information System (GIS) model with the ability to assess the seismic damage to pipelines for two well-known hazards, including ground shaking and ground failure simultaneously. The model that is developed and used in this study includes four main parts of database implementation, seismic hazard analysis, vulnerability assessment and seismic damage assessment to determine the pipeline’s damage probability. This model was implemented for main water distribution pipelines of Iran and tested for two different earthquake scenarios. The final damage probability of pipelines was estimated to be about 74% for water distribution pipelines of Mashhad including 40% and 34% for leak and break, respectively. In the next step, the impact of each earthquake input parameter on this model was extracted, and each of the three parameters had a huge impact on changing the results of pipelines’ damage probability. Finally, the dependency of the model in liquefaction susceptibility, landslide susceptibility, vulnerability functions and segment length was checked out and specified that the model is sensitive just to liquefaction susceptibility and vulnerability functions.


Introduction
Over the past decades, casualties due to natural disasters have increased in the world.One of the most important natural disasters is the earthquake, the neglecting of which leads to irreparable damages [1].Iran's position on the Alpine-Himalayan seismic belt puts this country among one of the world's earthquake-prone countries, and according to geological studies, about 97% of Iran's cities and villages are subjected to earthquake danger [2].In these earthquakes, most damages are related to critical infrastructures; and the study of the damages' history due to the world's destructive earthquakes indicates the importance of lifelines' failure's share of the losses.For example, in the San Fernando earthquake, although there were structural damages, the major damages were related to lifelines' damages, especially pipelines and bridges [3].
Lifelines are crucial lines for the survival of urbanization in the modern world.These lifelines are used to produce and distribute goods and services in urban units, and life possibility in cities depends on the quality and the quantity of these lifelines' function.In the case of cities, supply, transmission and distribution of water play an important role, and the water systems are considered as critical infrastructures of any society, since they are directly related to the vital and basic needs of society's population.Regarding the high prevalence and the wide range of affected people, the damage to the main water distribution pipelines would have great potential for creating massive crises in the society.In other words, small and short-time fluctuations in the main water distribution pipelines' service have wide effects on the performance of the entire urban water distribution network.
The vulnerability of pipelines against disasters is the subject of many studies these days, as the damages to these infrastructures not only is a part of damages, but also can complicate the disasters of other damages and significantly affect the recovery period.Various parameters, such as Peak Ground Acceleration (PGA), Peak Ground Velocity (PGV) and Permanent Ground Displacement (PGD) should be calculated to determine the vulnerability.In lifeline damage analysis, providing spatial information at different levels, such as geological maps, topography maps and lifelines' network of critical infrastructures, is more important than before.Spatial information science according to the Geospatial Information System (GIS) is a method for computerized organization of spatial data that creates a backup system for disaster management by spatial analysis [4].
According to this, the damages to pipelines (especially water pipelines) in the past earthquakes are reviewed, and various studies for pipelines' failure estimation done in the world will be investigated.Afterward, the methodology for damage probability analysis is presented in detail, step by step.It should be noted that the most important limitation of this model is the use of empirical equations to estimate the vulnerability of pipelines that have been used for damages of the previous earthquake.In order to test our model, two different seismic scenarios for the water distribution pipelines of Mashhad city are evaluated.In the following, a sensitivity analysis is performed to find the most effective parameters on pipelines' failure, and finally, the research conclusion is discussed.

Literature Review
The evaluation of the results and negative effects of past earthquakes is a good starting point in order to perform different actions to increase the level of networks' security.The earthquakes occurring in the past years have imposed irreparable damages to pipelines, but there is not enough information and statistics about earthquake damages and the damages to pipelines in the cities.The most important damages to water pipelines from past earthquakes can be summarized as follows [5]:

•
Philippines Mindanao earthquake (1976): a part of the main city water supply pipe was broken due to the collapse of a bridge on the pipelines.

•
El Salvador earthquake (1986): about 80 km of water pipes (20% of total pipelines) and about 65 km of sewage pipelines (22% of total sewage pipelines) were severely damaged.In this earthquake, a part of steel pipe was deformed and tilted.

•
Venezuela Caracas earthquake (1997): drinking water supply pipelines, which had 30 degrees along the fault, were severely damaged.Buried pipelines in soil and sewage treatment installations were damaged, as well.

•
Bam earthquake (2003): in this earthquake, 11 wells were damaged, and six points of the main pipes in central areas of the city were seriously damaged.
After several studies and reviews on the effects of earthquake damage on the pipelines, it can be concluded that damage to the pipelines caused by ground shaking and ground failure needs more attention and research.The ground failure will emerge as three types: faulting, landslide and liquefaction.Furthermore, liquefaction failure is divided into two categories as ground settlement and lateral spreading [6].The ground shaking and ground failure also are known as primary and secondary earthquake damages, respectively.Generally, pipes are under the combined effects of both types of hazards; so that O'Rourke and Deyoe [7] indicated that almost half of the pipe failures in the San Francisco earthquake in 1906 were due to the lateral spread process of liquefaction, and the other half of the damages occurred due to waves crossing and seismic moments in a larger region.Hence, it is very important to provide appropriate vulnerability functions in damage assessment.
A number of pipeline damages from six different earthquakes (four of them in Japan, one in Nicaragua and the last one in the United States) are studied by Katayama et al. [8] to estimate the fragility relations for both Asbestos Cement (AC) and segmented Cast Iron (CI) pipelines in terms of PGA.In this study, three different soil conditions, poor, average and good, were considered in fragility relations analysis.Barenberg [9] published the first PGV-based fragility relation for buried CI pipelines in the late 1980s using the data received from three different earthquakes in the U.S.
A comprehensive study on seismic loss estimation for water resource systems was published by ASCE in 1991 [10].In this study, PGA-based fragility relations were calculated using the collected damage data by Katayama et al. [8] and the Coalinga pipeline damage in 1983.In another study, a new pipeline fragility relation in terms of PGV was presented by O'Rourke and Ayala [11] based on the damage data of Barenberg [9], as well as the available damage information from three other earthquakes in the U.S. It is important to note that the presented fragility relation by O'Rourke and Ayala [11] has been used in the loss assessment methodology HAZUS-MH of FEMA [6].A few years later, a new GIS-based method was introduced by O'Rourke et al. [12] to examine various factors that affect the Los Angeles Department of Water and Power (LADWP) and the Metropolitan Water District (MWD) water supply services after the Northridge earthquake in 1994.In this research, seven seismic parameters including PGA, PGV, PGD, Modified Mercalli Intensity (MMI), Arias Intensity (AI), Spectral Acceleration (SA) and Seismic Intensity (SI) are used to study the relationship between the earthquake repair rate and the seismic intensity.In addition to that, the fragility relations for pipelines also were assessed in terms of MMI, SI, PGA and PGV.O'Rourke et al. [12] concluded that the best-related factor to the pipeline damage is PGV, and hence, the PGV-based fragilities are proposed for steel, CI, Ductile Iron (DI) and AC pipelines.
In another study, a set of algorithms is presented by the American Lifeline Alliance (ALA) to compute the probability of earthquake damage that affects various components of water supply systems [13].The analyses of water resources' pipeline damages related to earthquakes are continued by Pineda-Porras and Ordaz-Schroeder [14] using the data from the Michoacan earthquake (1985) in Mexico City.A detailed PGV map is used to find the best relationship between repair rate and seismic intensity.The results showed that a PGV-based fragility relation is the best option for Mexico City's water system [15].
As can be seen, various studies and relationships were presented along with several experimental vulnerability functions in the name of the equilibrium of repair rate per kilometer for earthquake damage analysis.Table 1 shows the relationships between various vulnerability functions and earthquake data used.In addition to that, strong-motion parameters considered by each study are shown in this table [16].
Moreover, several studies about pipelines' vulnerability assessment have been performed to provide the vulnerability relationships.Rahnama et al. [17] assessed the seismic vulnerability of 11 regional water distribution networks in the city of Tehran.In this study, the uncertainty caused by the earthquake was not considered, and just the damages to the water distribution network were evaluated.The most important results of the study indicated that all 11 regional water distribution network would not face a lot of damages, and just some parts of inflexible pipes with a small diameter would have a problem if the North Tehran Fault activates.
Omidvar et al. [3] presented a model for buried fuel pipelines, and after the determination of the earthquake scenario, the repair rate of the pipeline and the number of fires due to the earthquake around the pipeline were estimated.The model just evaluated the hazards of ground shaking, as the liquefaction probability in this area was extremely low.Eskandari et al. [4] presented a model in order to analyze the seismic damage of critical infrastructures according to a geospatial information system.The model included both types of hazards (ground shaking and the ground failure), and also the uncertainties at a reasonable level are considered.However, the model was used only for hard and heavy components of lifelines (such as refinery, substation, water pumping stations, etc.), and it cannot be used for pipelines' evaluation.
The differences between this study and current literature in this field briefly are: • Although some studies use the geospatial information system in order to analyze pipelines' damage, more evaluation was based on the instrumental view and multi-criteria decision making process than paying attention to simultaneous spatial assessments.

•
In most previous studies in this field, both types of hazards (ground failure and ground shaking) are not considered simultaneously after an earthquake, while here, both types of hazards are considered and combined appropriately with each other.

•
In some studies, the pipeline repair rate just was evaluated on the basis of experimental relationships, and the earthquake position in a pipeline is unknown, while the presented model estimates the position of damage with a high percentage.

•
Uncertainty analysis is another issue in the current literature, while earthquake analysis always comes with uncertainties.This study assesses the uncertainties with many repetitions in software and approximates the assessment results to the reality status.
• By the use of the sensitivity analysis of the model, the study determines the parameters that have many effects on pipelines' damage that have not been considered in the previous studies.

Materials and Methods
According to the complexities in the method of pipelines' damage assessment, the number of parameters and also the information layers that play an important role in damage analysis, modeling using an appropriate system like the geospatial information system is required.The general algorithm and conceptual model of this software system are shown in Figure 1.It is important to note that the visual studio with C# programming language along with .NET Frameworks and Arc GIS Engine for reading and writing files and processing on spatial data are used in this study.In the following, the required steps for the implementation of the system in the seismic damage assessment of pipelines are explained.

Case Study and Collecting Data in the Database
The local database filed by the ESRI platform is used in this study because the database should be able to support the spatial data (both vector and raster) and also be able to make a connection with the programming area.In order to implement the existing system and display outputs, a case study is required, so the data of Mashhad city is used in this research.
Mashhad is the second largest city in Iran located in an area between the Kopet Dagh Mountains from north and northeast and the Binalud Mountains from the southwest.Earthquake evaluation of this area shows intensive activities in past centuries, especially in the eighteenth century and then a relatively calm period [18].The existence of capable and active faults on both sides of the plain near this city show the high potential of earthquake hazard in this area.From east and southeast, this city has about a 20-km distance from a fault and from south and southwest less than 2 km from another fault.In 2007, in the urban complex of Mashhad, 275 micro earthquakes and earthquakes occurred, three of which were above 4.5 on the Richter scale, and the most important was 6.6 on the Richter scale [19].Therefore, there are capable and active faults around Mashhad, as is shown in Figure 2. The water demand of this city is supplied by 270 deep wells, 5 lime wells and 3 dams (Kardeh, Targh and Doosti).In this study, the damages to the main water distribution pipelines in Zones 7 and 8 of this city due to the earthquake are analyzed.The main water distribution pipelines in Zones 7 and 8 of this city include 9.2 km with asbestos cement, 16.4 km with ductile iron and 6.8 km with welded steel, which provide the water requirements of these two zones.The main layers of the database in this model include vector layers of water distribution pipelines, the descriptive layer of each pipeline's specifications, such as pipe material, pipe diameter, The main layers of the database in this model include vector layers of water distribution pipelines, the descriptive layer of each pipeline's specifications, such as pipe material, pipe diameter, etc., the vector layer of faults' position and historical and instrumental earthquakes around the studied area [20,21], the geology vector layer of the area, such as the vector layers of groundwater depth, deposit type, geological age, slope angle and the groundwater condition, and determining the map of liquefaction susceptibility.

Seismic Hazard Analysis
The seismic hazard analysis process is an important step in the seismic damage assessment on buildings and lifelines.In this study, the damages caused by both ground shaking and ground failure are analyzed [6], and the following steps are taken:

Selection of the Earthquake Scenario
Earthquake scenario selection has an important role in the analysis of outputs, and hence, the values of the earthquake scenario should be considered in a probable and pessimistic manner [22].Two scenarios are selected for this analysis in which the first one is on an active fault of the area (Figure 2) and the other one is selected randomly.Scenario 1 (on Shandiz fault): In this step, the maximum magnitude expected for Shandiz fault, which is a strike-slip fault with an 85-km failure length, is determined after preparing the seismic database of the area, according to presented relationships in Table 2.In addition to that, the focal depth is considered according to the average of the historical and instrumental earthquakes around each fault, and the earthquake position is also selected randomly on the Shandiz fault line in each analysis.Scenario 2 (selected randomly): According to 10 historical earthquakes and 2 instrumental earthquakes between 6 and 7.5 Richter in less than 100 km of the city [23], the random probable magnitude is selected in a range between 6 and 7.5 Richter on the basis of the uniform distribution function.The earthquake focal depth also is chosen in a range between 3 to 25 km because more than 85% of previous earthquakes occurred in this range [23].In order to determine the coordinates of an earthquake location, an area of the map is selected with a radius of 50 km around the city.The summaries of the two earthquake scenarios are shown in Table 3. Randomly in every time of analysis, a position is selected within the radius of 50 km from the city.

Pipelines Separation into Equal Length to Increase the Accuracy of Analysis
First of all, the pipelines are separated into small segments (1000 m) in order to increase the accuracy of analysis.Then, after recognizing the attenuation relationship and the appropriate vulnerability functions for each probable scenario, the hazard analysis for each segment is performed, and the output parameter values of hazard analysis for both ground shaking and ground failure in each segment of pipe are calculated separately.The main reasons for pipelines' segmentation are to accurately calculate the distance of each part of the pipe with the earthquake position, and the geological characteristics, as shown in Figure 3, in each small region varies continuously; so if the pipelines have a long length, the extraction of the geological parameters are faced with many uncertainties.Randomly in every time of analysis, a position is selected within the radius of 50 km from the city.

Pipelines Separation into Equal Length to Increase the Accuracy of Analysis
First of all, the pipelines are separated into small segments (1000 m) in order to increase the accuracy of analysis.Then, after recognizing the attenuation relationship and the appropriate vulnerability functions for each probable scenario, the hazard analysis for each segment is performed, and the output parameter values of hazard analysis for both ground shaking and ground failure in each segment of pipe are calculated separately.The main reasons for pipelines' segmentation are to accurately calculate the distance of each part of the pipe with the earthquake position, and the geological characteristics, as shown in Figure 3, in each small region varies continuously; so if the pipelines have a long length, the extraction of the geological parameters are faced with many uncertainties.

Hazard Analysis of Ground Shaking Based on Attenuation Relationships
The earthquake output parameter can be calculated as a function of earthquake magnitude, the distance between the earthquake source to the site, the site structure conditions and the fault type of the area for each segment of the pipeline, according to the attenuation relationship for each earthquake scenario.In this study, in order to reduce the uncertainties, 3 attenuation relationships including Zare et al. [29], Ghodrati et al. [30], and Campbell and Bozorgnia [31] are used, and the relevant relationships have been coded in the prepared system.

Hazard Analysis of Ground Failure
The ground failure is divided into the three main following categories: liquefaction, landslide, and faulting:

•
Liquefaction is the most important hazard due to ground failure that always threatens the pipelines.In this model, in order to consider the failure caused by soil liquefaction, at the beginning, the geologic map is drawn, and the liquefaction susceptibility of each segment of pipeline is determined.The liquefaction susceptibility of various types of soil deposits is estimated by assigning a qualitative susceptibility rating based on the general depositional environment, the geologic age of the deposit and the material type [32].Then, the liquefaction

Hazard Analysis of Ground Shaking Based on Attenuation Relationships
The earthquake output parameter can be calculated as a function of earthquake magnitude, the distance between the earthquake source to the site, the site structure conditions and the fault type of the area for each segment of the pipeline, according to the attenuation relationship for each earthquake scenario.In this study, in order to reduce the uncertainties, 3 attenuation relationships including Zare et al. [29], Ghodrati et al. [30], and Campbell and Bozorgnia [31] are used, and the relevant relationships have been coded in the prepared system.

Hazard Analysis of Ground Failure
The ground failure is divided into the three main following categories: liquefaction, landslide, and faulting:

•
Liquefaction is the most important hazard due to ground failure that always threatens the pipelines.In this model, in order to consider the failure caused by soil liquefaction, at the beginning, the geologic map is drawn, and the liquefaction susceptibility of each segment of pipeline is determined.The liquefaction susceptibility of various types of soil deposits is estimated by assigning a qualitative susceptibility rating based on the general depositional environment, the geologic age of the deposit and the material type [32].Then, the liquefaction probability for each segment of pipeline is calculated.The likelihood of experiencing liquefaction at a specific location is primarily influenced by the susceptibility of the soil, the amplitude and duration of ground shaking and the depth of groundwater.The probability of liquefaction for a given susceptibility category can be determined using the following relationship [6]: where P [Liquefaction SC | PGA = a] is the conditional liquefaction probability for a given susceptibility category at a specified level of peak ground acceleration, KM is the moment magnitude correction factor, Kw is the ground water correction factor [33,34] and P ml is the proportion of the map unit susceptible [35].
At the end, the permanent ground displacement caused by the two liquefaction hazards, lateral spreading and ground settlement is determined [36].

•
The landslide hazard evaluation requires the characterization of the landslide susceptibility of the soil/geologic conditions of a region or sub-region.Susceptibility is characterized by the geologic group, slope angle and critical acceleration [37].Landslide susceptibility is measured on a scale of 1 to 10.The site condition is analyzed using three geologic groups and groundwater level.The groundwater condition is divided into either dry condition (groundwater below level of the sliding) or wet condition (groundwater level at ground surface).The critical acceleration is then estimated for the respective geologic and groundwater conditions and the slope angle [37].Furthermore, the percentage of map area having a landslide-susceptible deposit is estimated [38]; and finally, the permanent ground displacements are determined based on the expected displacement correction factor [39], the induced acceleration and the number of cycles [33].

•
The median maximum displacement is estimated based on the earthquake moment magnitude [27].
Finally, the output parameters of ground shaking (PGV and PGD) and the output parameters of ground failure (liquefaction probability, the permanent ground displacement caused by lateral spreading, ground settlement, landslide and faulting) are evaluated for each part of the pipeline.

Seismic Vulnerability Functions of Pipelines
The vulnerability function is a relationship that expresses the expected failure for a structure as a function of strong ground movement or any other seismic parameter and also can be expressed in different states such as vulnerability relationships and failure curves.Generally, the damage to pipelines is expressed as a function of Peak Ground Velocity (PGV) or Permanent Ground Displacement (PGD) or in the parameter of the Repair Rate (RR) [6].The repair rate is presented by different models, such as ALA (American Lifeline Alliance) [13], LESSLOSS (Risk Mitigation for Earthquakes and Landslides Integrated Project) [40] and the Federal Emergency Management Agency [6].In this study, the Federal Emergency Management Agency [6] is used (Table 4).The repair rate relationship due to PGV was formed according to experimental data of O'Rourke and Ayala [11] and on the basis of four observed earthquakes in America and two earthquakes in Mexico, as well as the repair rate relationship due to Permanent Ground Displacement (PGD) was extracted from Honegger and Eguchi [41] results for San Diego.The RR variable is the repair rate and represents the number of repairs in each kilometer of a pipeline [6].The pipelines' damage probability is calculated based on the following experimental relationship (after determination of the repair rate assuming that the number of pipe failures follows Poisson's distribution) [42]: In this relationship, L is the pipeline length in km, RR is the pipeline repair rate and P f is the pipeline failure probability.

Seismic Damage Assessment of Pipelines
Damages to pipelines are either leak mode or break mode in each segment of the pipeline.In order to analyze the damage to network pipelines, the following steps should be taken:

•
First step: the output values' extraction of seismic the damage assessment parameter (PGA, PGV, and PGD due to liquefaction or landslide) for each segment of pipelines.

•
Second step: calculation of the repair rate and failure probability for each segment of pipelines.

•
Third step: failure probability combination due to ground shaking or ground failure.
After determination of the failure probability due to ground shaking and ground failure for each segment, the failure total probability is combined based on union and intersection probability rules of 2 events as the following relationship: • P COMB : the final damage probability • P GF : failure probability due to ground shaking • P GS : failure probability due to ground failure After determination of the total failure probability for each segment, the following operation is repeated:

•
First, a random value is chosen between zero and one.If the final damage probability for the segment is more than the selected random value, the segment fails and goes to the next step.

•
Each pipe segment for which its final damage probability is higher than the selected value has failed.Now, in this step, it should be determined that the segment breaks or leaks.Based on the studies of O'Rourke and Ayala [11], the hazard effect due to ground shaking (PGV) in pipelines occurs as 80% leak and 20% break, and the hazard effect due to ground failure and permanent displacements (PGD) in pipelines occurs as 80% break and 20% leak.Therefore, if we assume that the failure probability due to ground shaking is PGS and the failure probability due to ground failure is PGF, then: If P GS > P GF , a random value is chosen between 0 and 100.If the chosen value for the pipe segment was between 0 and 80, the segment has leaked; otherwise, the segment has been broken.

•
If P GS < P GF , a random value is chosen between 0 and 100.If the chosen value for the pipe segment was between 0 and 80, the segment has been broken; otherwise, the segment has leaked.
Therefore, in each earthquake scenario, the pipeline part that is subjected to break or leak can be determined.In this study, all operations are repeated ten thousand-times based on Monte Carlo simulation, and then, the achieved results for each pipeline segment are averaged.Finally, the numbers of leaks, breaks, leak rate, break rate and the final failure probability are determined for each pipeline segment based on ten thousand-times repetition.Based on the leak rate and break rate outputs for each segment, the damage states can be determined (Table 5).All of these steps for all pipeline segments are repeated 10,000 times based on Monte Carlo simulation to consider the uncertainty of the system, as well (Figure 4).Monte Carlo methods are especially useful in modeling systems with many random variables with known probability distributions.Monte Carlo simulation, which is based on the iteration and generation of random variables from a specific range, is one of the most widely-used numerical methods [44].
ISPRS Int.J. Geo-Inf.2017, 6, 169 11 of 19 simulation, and then, the achieved results for each pipeline segment are averaged.Finally, the numbers of leaks, breaks, leak rate, break rate and the final failure probability are determined for each pipeline segment based on ten thousand-times repetition.Based on the leak rate and break rate outputs for each segment, the damage states can be determined (Table 5).All of these steps for all pipeline segments are repeated 10,000 times based on Monte Carlo simulation to consider the uncertainty of the system, as well (Figure 4).Monte Carlo methods are especially useful in modeling systems with many random variables with known probability distributions.Monte Carlo simulation, which is based on the iteration and generation of random variables from a specific range, is one of the most widely-used numerical methods [44].

Results
In this study, the seismic hazard analysis and the seismic damage analysis of water distribution pipelines are discussed.The desired analyses are performed for two scenarios of Table 3.As it is not possible to show all of the outcomes, the results of hazard analysis are not shown here, and only the city water distribution pipelines' damage is presented.In the next step, the effects of various earthquake parameters on the pipelines, as well as the sensitivity of input parameters are evaluated.
The zoning map of the pipe leak rate and pipe break rate for Scenarios 1 and 2 is shown in Figures 5 and 6, respectively.As can be seen in these figures, the amount of damage in Scenario 1 is more than the expected damage from Scenario 2.Moreover, it can be concluded that a segment cannot achieve a leak rate over 50% and a break rate over 50% simultaneously.The results obviously showed that the leak rate is higher than that of the break in various segments.However, the pipeline break rate is higher than the leak rate as the displacement due to liquefaction or landslide is greater than PGV, and the dominant hazard is the ground failure hazard regarding the high susceptibility of landslide or liquefaction.

Results
In this study, the seismic hazard analysis and the seismic damage analysis of water distribution pipelines are discussed.The desired analyses are performed for two scenarios of Table 3.As it is not possible to show all of the outcomes, the results of hazard analysis are not shown here, and only the city water distribution pipelines' damage is presented.In the next step, the effects of various earthquake parameters on the pipelines, as well as the sensitivity of input parameters are evaluated.
The zoning map of the pipe leak rate and pipe break rate for Scenarios 1 and 2 is shown in Figures 5  and 6, respectively.As can be seen in these figures, the amount of damage in Scenario 1 is more than the expected damage from Scenario 2.Moreover, it can be concluded that a segment cannot achieve a leak rate over 50% and a break rate over 50% simultaneously.The results obviously showed that the leak rate is higher than that of the break in various segments.However, the pipeline break rate is higher than the leak rate as the displacement due to liquefaction or landslide is greater than PGV, and the dominant hazard is the ground failure hazard regarding the high susceptibility of landslide or liquefaction.The zoning map of pipeline damage probability based on Scenario 2 is shown in Figure 7.As can be seen in this figure, the pipeline damage probability in Scenario 1 is higher than Scenario 2, as the leak and break rates in Scenario 2 are less than Scenario 1.The damage in Scenario 2 also is less than Scenario 1 since the distance of earthquake location to the pipeline is not constant, and it varies in each analysis.This is the reason that the fault distance to the pipeline for Shandiz fault is less than the case we have for Scenario 2 as its location is selected randomly in a range of 50 km around the city.
The results of damage analysis for the main water distribution of Mashhad city are presented in the form of a chart in Figure 8.As can be seen from the figure, the leak percentage in Scenarios 1 and 2 is about 13% and 21% higher than the break percentage of pipelines, respectively.Therefore, it can be concluded that the break percentage is nearly equal to the leak percentage in the analyses, and this is due to moderate to high susceptibility of liquefaction in the area.On the other hand, the pipes' break rate is very important for pipeline damage analysis because about 45% of pipelines' damage is dependent on the pipes' break rate.The zoning map of pipeline damage probability based on Scenario 2 is shown in Figure 7.As can be seen in this figure, the pipeline damage probability in Scenario 1 is higher than Scenario 2, as the leak and break rates in Scenario 2 are less than Scenario 1.The damage in Scenario 2 also is less than Scenario 1 since the distance of earthquake location to the pipeline is not constant, and it varies in each analysis.This is the reason that the fault distance to the pipeline for Shandiz fault is less than the case we have for Scenario 2 as its location is selected randomly in a range of 50 km around the city.
The results of damage analysis for the main water distribution of Mashhad city are presented in the form of a chart in Figure 8.As can be seen from the figure, the leak percentage in Scenarios 1 and 2 is about 13% and 21% higher than the break percentage of pipelines, respectively.Therefore, it can be concluded that the break percentage is nearly equal to the leak percentage in the analyses, and this is due to moderate to high susceptibility of liquefaction in the area.On the other hand, the pipes' break rate is very important for pipeline damage analysis because about 45% of pipelines' damage is dependent on the pipes' break rate.The zoning map of pipeline damage probability based on Scenario 2 is shown in Figure 7.As can be seen in this figure, the pipeline damage probability in Scenario 1 is higher than Scenario 2, as the leak and break rates in Scenario 2 are less than Scenario 1.The damage in Scenario 2 also is less than Scenario 1 since the distance of earthquake location to the pipeline is not constant, and it varies in each analysis.This is the reason that the fault distance to the pipeline for Shandiz fault is less than the case we have for Scenario 2 as its location is selected randomly in a range of 50 km around the city.The zoning map of pipeline damage probability based on Scenario 2 is shown in Figure 7.As can be seen in this figure, the pipeline damage probability in Scenario 1 is higher than Scenario 2, as the leak and break rates in Scenario 2 are less than Scenario 1.The damage in Scenario 2 also is less than Scenario 1 since the distance of earthquake location to the pipeline is not constant, and it varies in each analysis.This is the reason that the fault distance to the pipeline for Shandiz fault is less than the case we have for Scenario 2 as its location is selected randomly in a range of 50 km around the city.
The results of damage analysis for the main water distribution of Mashhad city are presented in the form of a chart in Figure 8.As can be seen from the figure, the leak percentage in Scenarios 1 and 2 is about 13% and 21% higher than the break percentage of pipelines, respectively.Therefore, it can be concluded that the break percentage is nearly equal to the leak percentage in the analyses, and this is due to moderate to high susceptibility of liquefaction in the area.On the other hand, the pipes' break rate is very important for pipeline damage analysis because about 45% of pipelines' damage is dependent on the pipes' break rate.The results of damage analysis for the main water distribution of Mashhad city are presented in the form of a chart in Figure 8.As can be seen from the figure, the leak percentage in Scenarios 1 and 2 is about 13% and 21% higher than the break percentage of pipelines, respectively.Therefore, it can be concluded that the break percentage is nearly equal to the leak percentage in the analyses, and this is due to moderate to high susceptibility of liquefaction in the area.On the other hand, the pipes' break rate is very important for pipeline damage analysis because about 45% of pipelines' damage is dependent on the pipes' break rate.

Discussion
The main parameters and relationships that affect the pipeline damage analysis are discussed in this section.

Determine the Effect of Earthquake Input Parameters
In this step, the effect of earthquake input parameters on the probability of damage to pipelines is determined.In other words, by changing the magnitude of leak percentage, the break percentage and the final failure, the probability of failure is measured (Figure 9).The probability of failure does not change too much for magnitudes less than 5 Richter, while in the range of 4 to 6 Richter, the break rate is higher than the leak rate.It is important to note that the failure probability of the leak rate increases and would be higher than the break rate for magnitudes more than 6.30.The reason is that the ground failure hazard is more dependent on ground conditions, and so, with the increasing magnitude, the ground condition is not changing.The ground shaking that effects the leak values is highly dependent on the earthquake input parameters; and by increasing magnitude, the PGV value becomes more than the PGD, and the repair rate becomes more than the break rate.The leak percentage, the break percentage and the pipelines final damage probability are proportional to the change in focal depth, focal distance, PGV and PGD, respectively (Figure 10A to Figure 10D).Based on Figure 10A, if the focal depth is more than 50 km, the damage probability is very small.According to Figure 10B, the damage gradient at a focal distance less than 30 km is small, and the results are very similar to each other; while with the increasing distance, the gradient is increased.Figure 10C,D shows the changes in the output parameters of the seismic hazard analysis

Discussion
The main parameters and relationships that affect the pipeline damage analysis are discussed in this section.

Determine the Effect of Earthquake Input Parameters
In this step, the effect of earthquake input parameters on the probability of damage to pipelines is determined.In other words, by changing the magnitude of leak percentage, the break percentage and the final failure, the probability of failure is measured (Figure 9).The probability of failure does not change too much for magnitudes less than 5 Richter, while in the range of 4 to 6 Richter, the break rate is higher than the leak rate.It is important to note that the failure probability of the leak rate increases and would be higher than the break rate for magnitudes more than 6.30.The reason is that the ground failure hazard is more dependent on ground conditions, and so, with the increasing magnitude, the ground condition is not changing.The ground shaking that effects the leak values is highly dependent on the earthquake input parameters; and by increasing magnitude, the PGV value becomes more than the PGD, and the repair rate becomes more than the break rate.

Discussion
The main parameters and relationships that affect the pipeline damage analysis are discussed in this section.

Determine the Effect of Earthquake Input Parameters
In this step, the effect of earthquake input parameters on the probability of damage to pipelines is determined.In other words, by changing the magnitude of leak percentage, the break percentage and the final failure, the probability of failure is measured (Figure 9).The probability of failure does not change too much for magnitudes less than 5 Richter, while in the range of 4 to 6 Richter, the break rate is higher than the leak rate.It is important to note that the failure probability of the leak rate increases and would be higher than the break rate for magnitudes more than 6.30.The reason is that the ground failure hazard is more dependent on ground conditions, and so, with the increasing magnitude, the ground condition is not changing.The ground shaking that effects the leak values is highly dependent on the earthquake input parameters; and by increasing magnitude, the PGV value becomes more than the PGD, and the repair rate becomes more than the break rate.The leak percentage, the break percentage and the pipelines final damage probability are proportional to the change in focal depth, focal distance, PGV and PGD, respectively (Figure 10A to Figure 10D).Based on Figure 10A, if the focal depth is more than 50 km, the damage probability is very small.According to Figure 10B, the damage gradient at a focal distance less than 30 km is small, and the results are very similar to each other; while with the increasing distance, the gradient is increased.Figure 10C,D shows the changes in the output parameters of the seismic hazard analysis The leak percentage, the break percentage and the pipelines final damage probability are proportional to the change in focal depth, focal distance, PGV and PGD, respectively (Figure 10A to Figure 10D).Based on Figure 10A, if the focal depth is more than 50 km, the damage probability is very small.According to Figure 10B, the damage gradient at a focal distance less than 30 km is small, and the results are very similar to each other; while with the increasing distance, the gradient is increased.Figure 10C,D shows the changes in the output parameters of the seismic hazard analysis on the damage of the pipes.As it be seen in Figure 10C, the dominant damage is the leak percentage associated with the ground shaking hazard and the PGV; while in Figure 10D, the dominant damage is the break that is dependent on the ground failure hazard and the PGD.It should be noted that the pipeline damage probability is associated with the earthquake magnitude, focal depth, focal distance, PGV and PGD according to Figures 9 and 10.
ISPRS Int.J. Geo-Inf.2017, 6, 169 14 of 19 on the damage of the pipes.As it be seen in Figure 10C, the dominant damage is the leak percentage associated with the ground shaking hazard and the PGV; while in Figure 10D, the dominant damage is the break that is dependent on the ground failure hazard and the PGD.It should be noted that the pipeline damage probability is associated with the earthquake magnitude, focal depth, focal distance, PGV and PGD according to Figures 9 and 10.

Sensitivity Analysis
The sensitivity analysis is performed to assess the dependency level of the developed model in this study on the input parameters including liquefaction susceptibility, landslide susceptibility, pipelines vulnerability and pipeline length.
Based on the sensitivity analysis, by changing the liquefaction susceptibility and landslide between three different states (without susceptibility, moderate susceptibility and high susceptibility) and calculating values for the leak percentage, the break percentage and the pipelines final damage probability, it can be concluded that the liquefaction susceptibility is able to effect the break rate and the final damage probability, while it does not change the leak rate.The reason is that the leak rate is dependent on the ground shaking hazard while changing the ground conditions, only increasing the ground failure hazard.Figure 11 also shows that the outputs of the model are not sensitive to landslide susceptibility as the hazard due to the landslide is not considered in the relationships used in this study for pipelines' vulnerability.
the leak rate is dependent on the ground shaking hazard while changing the ground conditions, only increasing the ground failure hazard.Figure 11 also shows that the outputs of the model are not sensitive to landslide susceptibility as the hazard due to the landslide is not considered in the relationships used in this study for pipelines' vulnerability.
Figure 12A shows the sensitivity analysis of the parameters of the vulnerability function and pipeline length in the damage probability analysis.This model uses the Federal Emergency Management Agency [6] relationship for which the anticipated damage value is higher than Eidinger [45] and ALA (American Lifeline Alliance) [13] relationships and is somewhat pessimistic.Based on Figure 12A, the final damage probability in the FEMA [6] relationship is greater than the other and in the ALA [13] relationship is less than the other.Furthermore, the break rate value of the Eidinger [45] and ALA [13] relationships has been estimated with greater values in comparison to the FEMA [6] relationship.It can be concluded that the failure probability based on the FEMA [6] relationship with respect to the Eidinger [45] and ALA [13] relationships has been estimated to be about 19% and 47%, respectively.
In seismic hazard analysis, in order to consider the distance effect and the geological parameters in each location, the pipelines are divided into a number of the same-sized segments.Then, for each segment, the output parameters are calculated.In this study, each segment length has been considered as 1000 m as the default in the analyses.Therefore, in order to determine what would be the effect of the change in length of segments on the outputs, three segments with a length of 100 m, 500 m and 5000 m are re-calculated.After comparing the results, it was observed that changing the pipeline length has a small impact on the probability of damage.Figure 12A shows the sensitivity analysis of the parameters of the vulnerability function and pipeline length in the damage probability analysis.This model uses the Federal Emergency Management Agency [6] relationship for which the anticipated damage value is higher than Eidinger [45] and ALA (American Lifeline Alliance) [13] relationships and is somewhat pessimistic.Based on Figure 12A, the final damage probability in the FEMA [6] relationship is greater than the other and in the ALA [13] relationship is less than the other.Furthermore, the break rate value of the Eidinger [45] and ALA [13] relationships has been estimated with greater values in comparison to the FEMA [6] relationship.It can be concluded that the failure probability based on the FEMA [6] relationship with respect to the Eidinger [45] and ALA [13] relationships has been estimated to be about 19% and 47%, respectively.
In seismic hazard analysis, in order to consider the distance effect and the geological parameters in each location, the pipelines are divided into a number of the same-sized segments.Then, for each segment, the output parameters are calculated.In this study, each segment length has been considered as 1000 m as the default in the analyses.Therefore, in order to determine what would be the effect of the change in length of segments on the outputs, three segments with a length of 100 m, 500 m and 5000 m are re-calculated.After comparing the results, it was observed that changing the pipeline length has a small impact on the probability of damage.

Conclusions
The past earthquakes experiences have shown that if the damage to lifelines leads to stopping public services after disasters' occurrence, especially an earthquake in a city, the damages and injuries would be multiplied.Undoubtedly, the vulnerability analysis and investigation of pipelines after earthquakes can be effective in planning for disasters and for establishing an appropriate approach to prevent these events in major cities.In this study, a model is developed to perform numerous spatial analyses in the geospatial information system.The four main parts of the desired model include: the collecting data, seismic hazard analysis, vulnerability assessment and damages' assessment of pipelines.The Monte Carlo simulation also is applied to consider uncertainties in damage analysis.
The results of this study briefly are: • Although the model presented in this paper was used for water pipelines as a case study, it can easily be evaluated for other pipelines just by having the vulnerability function in accordance with the relevant pipeline.

•
The experimental vulnerability assessment models are only able to estimate the leak and break values in the pipeline system due to an earthquake.Therefore, in order to solve this problem, few relationships are presented in the methodology for damage assessment by using the geospatial information system, which can specify the failure locations.

•
The relationships suggested in this model are capable of analyzing the damage of pipelines for two hazards (ground shaking and ground failure) due to earthquake simultaneously, unlike the similar studies in the literature.Some studies have not considered the ground failure hazard, while the high effect of this hazard was investigated in this study.

•
In addition, most of the studies are only based on the experimental relationships of pipelines' vulnerability without considering the uncertainty in the results, while our model is considering the uncertainty based on Monte Carlo simulation.Hence, it is expected that the results are much closer to reality.

Conclusions
The past earthquakes experiences have shown that if the damage to lifelines leads to stopping public services after disasters' occurrence, especially an earthquake in a city, the damages and injuries would be multiplied.Undoubtedly, the vulnerability analysis and investigation of pipelines after earthquakes can be effective in planning for disasters and for establishing an appropriate approach to prevent these events in major cities.In this study, a model is developed to perform numerous spatial analyses in the geospatial information system.The four main parts of the desired model include: the collecting data, seismic hazard analysis, vulnerability assessment and damages' assessment of pipelines.The Monte Carlo simulation also is applied to consider uncertainties in damage analysis.
The results of this study briefly are: • Although the model presented in this paper was used for water pipelines as a case study, it can easily be evaluated for other pipelines just by having the vulnerability function in accordance with the relevant pipeline.

•
The experimental vulnerability assessment models are only able to estimate the leak and break values in the pipeline system due to an earthquake.Therefore, in order to solve this problem, few relationships are presented in the methodology for damage assessment by using the geospatial information system, which can specify the failure locations.

•
The relationships suggested in this model are capable of analyzing the damage of pipelines for two hazards (ground shaking and ground failure) due to earthquake simultaneously, unlike the similar studies in the literature.Some studies have not considered the ground failure hazard, while the high effect of this hazard was investigated in this study.

•
In addition, most of the studies are only based on the experimental relationships of pipelines' vulnerability without considering the uncertainty in the results, while our model is considering the uncertainty based on Monte Carlo simulation.Hence, it is expected that the results are much closer to reality.

19 Figure 1 .
Figure 1.The conceptual model for seismic damage assessment of pipelines using geospatial analysis.

Figure 1 .
Figure 1.The conceptual model for seismic damage assessment of pipelines using geospatial analysis.

Figure 1 .
Figure 1.The conceptual model for seismic damage assessment of pipelines using geospatial analysis.

Figure 2 .
Figure 2. Map of faults and water distribution pipelines.

Figure 2 .
Figure 2. Map of faults and water distribution pipelines.

Figure 4 .
Figure 4. Developed computer program for the seismic damage analysis of pipelines.

Figure 4 .
Figure 4. Developed computer program for the seismic damage analysis of pipelines.

Figure 5 .Figure 6 .
Figure 5. Zoning map of the leak rate in two probable scenarios.

(A) Scenario 1 (B) Scenario 2 Figure 7 .
Figure 7. Zoning map of the damage probability in two probable scenarios.

Figure 5 .Figure 5 .Figure 6 .
Figure 5. Zoning map of the leak rate in two probable scenarios.

(A) Scenario 1 (B) Scenario 2 Figure 7 .
Figure 7. Zoning map of the damage probability in two probable scenarios.

Figure 6 .
Figure 6.Zoning map of the break rate in two probable scenarios.

Figure 5 .Figure 6 .
Figure 5. Zoning map of the leak rate in two probable scenarios.

Figure 7 .
Figure 7. Zoning map of the damage probability in two probable scenarios.

Figure 7 .
Figure 7. Zoning map of the damage probability in two probable scenarios.

19 Figure 8 .
Figure 8. Summary of the damage results of water distribution pipelines in different scenarios.

Figure 9 .
Figure 9. Determine the impact of magnitude on the failure probability of pipelines.

2 Figure 8 .
Figure 8. Summary of the damage results of water distribution pipelines in different scenarios.

19 Figure 8 .
Figure 8. Summary of the damage results of water distribution pipelines in different scenarios.

Figure 9 .
Figure 9. Determine the impact of magnitude on the failure probability of pipelines.

Figure 9 .
Figure 9. Determine the impact of magnitude on the failure probability of pipelines.

Figure 10 .
Figure 10.The impact of various earthquake parameters on the failure probability of pipelines.

Figure 10 .
Figure 10.The impact of various earthquake parameters on the failure probability of pipelines.

Figure 11 .
Figure 11.The sensitivity analysis of the parameters of the liquefaction and landslide susceptibility in the failure probability of pipelines.

Figure 11 .
Figure 11.The sensitivity analysis of the parameters of the liquefaction and landslide susceptibility in the failure probability of pipelines.

Figure 12 .
Figure 12.The sensitivity analysis of the parameters of the vulnerability function and the pipeline length in the failure probability of pipelines.

Figure 12 .
Figure 12.The sensitivity analysis of the parameters of the vulnerability function and the pipeline length in the failure probability of pipelines.

Table 2 .
Determine the maximum magnitude earthquake due to activating the Shandiz fault based on different relationships.

Table 3 .
Summary of the two earthquake scenario parameters.

Table 3 .
Summary of the two earthquake scenario parameters.

Table 5 .
Classification of damage levels for water pipelines.

Table 5 .
Classification of damage levels for water pipelines.